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ABSTRACT 


Under  contract  to  the  Philadelphia  Corps  of  Engineers,  Camp  Dresser  &  McKee 
developed  a  Daily  Flow  Model  for  the  Delaware  River  Basin  from  its  headwaters 
to  the  Delaware  Memorial  Bridge.  A  task  committee  was  set  up  to  direct  the 
study.  Members  of  the  committee  included  persons  from  the  Philadelphia  District 
Office  of  the  U.S.  Amiy  Corps  of  Engineers,  the  Pennsylvania  Department  of 
Environmental  Resources,  the  New  York  Department  of  Environmental  Conservation, 
the  Delaware  River  Basin  Commission  and  the  United  States  Geological  Survey. 

All  decisions  concerning  model  operations  were  a  direct  result  of  committee 
participation.  As  a  preface  to  the  work  contained  in  this  report,  the  assump¬ 
tions  and  limitations  of  this  project  must  be  set  forth.  The  essence  of  the 
results  of  this  project  are  of  a  comparative  basis,  not  to  be  interpreted  as 
absolute  yields  of  the  systems  analyzed  as  they  are  presently  operated.  The 
following  paragraphs  define  the  report  structure  and  give  a  short  synopsis  of 
the  capabilities  of  the  Daily  Flow  Model. 

The  results  of  this  work  which  were  conducted  in  two  phases,  are  presented 
in  two  separate  volumes  in  this  report  and  five  appendices.  The  first  volume 
presents  Phase  I  of  the  work  and  contains  the  preliminary  work  of  developing 
50  years  of  historical  natural  inflows  for  selected  reaches  within  the  Delaware 
River  Basin.  It  also  includes  an  analysis  of  two  50  year  model  simulations 
of  the  flows  in  the  entire  river  basin.  One  excludes  all  reservoirs  and  simply 
routes  flows  through  the  basin  while  the  other  includes  the  operation  of  the 
three  New  York  City  reservoirs,  Pepacton,  Cannonsville  and  Neversink,  to  main¬ 
tain  a  low  flow  objective  at  Montaque,  N.J.,  while  at  the  same  time  allowing 
the  diversion  of  water  to  the  New  York  City  water  supply  system.  The  results 
of  these  two  simulations  are  presented  in  Appendices  A  and  B  respectively. 

The  Appendices  contain  statistical  low  flow  frequency  and  duration  analyses 
at  various  locations  within  the  basin. 

The  second  phase  of  this  work  is  described  in  Volume  II  of  this  report 
and  contains  an  analysis  which  simulates  the  use  of  both  existing  and  pro¬ 
posed  reservoirs  to  maintain  various  flow  conditions  at  Trenton,  N.J.  The 
results  include  a  matrix  of  the  maximum  maintainable  low  flow  values  of  vari¬ 
ous  return  frequencies  which  correspond  to  17  different  combinations  of  these 


reservoirs.  The  combinations  range  from  a  single  reservoir  online  to  up  to 
six  reservoirs  online  to  maintain  the  Trenton  flow  objectives  under  different 
hydrologic  regimes.  Three  of  the  17  combinations  were  also  analyzed  using  the 
low- flow  frequency  and  duration  statistics  as  in  Phase  I.  The  results  of  the 
three  simulations  are  presented  in  Appendices  C,  D  and  E. 

The  Daily  Flow  Model  is  capable  of  handling  various  basin  operations. 

First  and  foremost  the  model  routes  flows  from  upstream  to  downstream  using 
lag  functions  derived  from  actual  maximum  and  minimum  flow  times  through 
reaches  in  the  watershed.  Fifty  (50)  years  of  historical  daily  flow  data 

were  used  to  develop  incremental  inflows  at  61  key  locations  in  the  basin. 

In  order  to  develop  the  incremental  inflows  it  was  necessary  to  develop 

complete  fifty  year  flow  records  at  all  USGS  gaging  stations  in  the  basin. 

Incomplete  or  regulated  records  were  filled-in  and  extended  by  correlation 
to  a  nearby,  long-term  recording  station. 

In  addition  to  routing  flows,  the  model  can  operate  up  to  ten  reservoirs 
to  augment  low  flow  conditions  at  two  locations  by  establishing  minimum  flow 
targets.  In  this  report,  Montaque  and  Trenton,  New  Jersey  were  used  as  these 
locations.  Three  reservoirs  in  New  York,  Pepacton,  Cannonsville  and  Neversink, 
were  used  to  augment  Montaque's  flow  according  to  the  1954  Supreme  Court  Decree. 
All  other  reservoirs,  whether  online  or  proposed,  were  used  to  augment  the  flow 
at  Trenton. 

Certain  simulation  procedures,  particularly  concerning  reservoir  operation, 
were  adapted  in  the  model  that  tend  to  cause  the  simulation  to  overestimate  the 
real  world  obtainable  system  yields.  Most  of  these  simplifications  of  the  real 
operating  decisions  were  made  in  order  to  limit  the  amount  of  computer  time 
required  to  execute  the  model.  This  constraint  is  a  very  practical  concern  when 
a  model  is  developed  to  handle  fifty  years  of  daily  flow  data  at  numerous  loca¬ 
tions  in  such  a  large  watershed.  Perhaps,  the  most  important  model  operation 
decision  concerns  the  manner  in  which  the  simulation  procedure  makes  decisions 
on  reservoir  releases.  For  example,  when  a  release  is  needed  to  augment  any 
target  flow,  the  release  is  taken  from  the  fullest  reservoir  using  percentage 
as  the  criteria.  This  release  scheme  does  not  take  into  account  hydropower 
optimization  or  water  quality  effects  as  is  done  be  the  controlling  parties  in  the 
Delaware  basin  in  the  real-life  situation.  Other  adapted  procedures  that 


affect  the  simulation  results  are  the  following: 

•  only  50  years  of  exact  historical  hydrology  are  used  for 
simulation  of  future  conditions; 

•  "perfect  foresight"  of  the  next  days  flow  conditions  is 
used  in  operating  the  NYC  reservoirs; 

•  reservoirs  are  operated  with  changing  releases  on  a  daily 
basis; 

•  releases  are  always  made  from  the  percent  fullest  reservoir 
resulting  in  an  exactly  even  draw  down  of  all  facilities;  and 

•  no  considerations  are  made  in  the  model  of  hydropower  optimi¬ 
zation,  water  quality  conditions  or  existing  physical  limita¬ 
tions  such  as  pipe  capacity  when  releasing  from  reservoirs. 

Even  with  these  limitations,  the  model  has  been  proven  to  be  a  powerful 
analysis  tool  for  Delaware  basin  planners.  It  should  be  used  as  a  comparative 
analysis  tool,  not  necessarily  reflecting  the  exact  ultimate  yields  of  the 
system  as  it  is  presently  operated.  However,  the  model  can  be  made  to  more 
closely  simulate  with  relatively  minor  modifications,  which  may  have  harsh 
affects  on  simulation  execution  times,  the  present  operations  of  the  system. 

Phase  II  of  this  report  goes  into  much  more  detail  concerning  some  of  the 
useful  applications  listed  below: 

•  analysis  of  bringing  different  reservoirs  on  and  off  line 
to  augment  low  flow  conditions; 

•  the  effects  of  using  different  operations  for  the  NYC 
reservoirs; 

•  the  hydrologic  possibilities  and  consequences  of  maintaining 
different  low  flow  objectives  at  various  locations;  and 

•  determining  the  effects  of  constructing  peak  flow  skimming 
facilities  on  the  main  stem  or  tributaries  to  the  Delaware. 

The  Daily  Flow  Model  is  still  relatively  new  and  in  the  development  phase. 

It  has  been  coded  for  the  Delaware  River  Basin  and  would  require  modification 
to  be  used  on  any  other  watershed.  The  Fortran  coding  is  not  throughly  doc¬ 
umented  nor  is  the  code  itself  straight  forward  and  clean.  The  model  was 
developed  and  used  as  an  analysis  tool  to  solve  a  particular  problem  and  that 
it  did  exceedingly  well.  The  purpose  of  this  project  was  not  the  development 
of  a  generalized  computer  program,  although  with  minor  changes  the  model  could 
fairly  easily  be  adapted  to  simulate  operations  in  other  river  basins  where 
daily  information  has  a  utility.  For  the  Delaware  basin  itself  the  model  can 
be  improved  upon  to  make  it  more  reflective  of  actual  field  operating  conditions. 


For  example,  the  model's  reservoirs  can  be  changed  to  operate  more  closely 
to  the  real-life  operations.  Hydropower  optimization  and  outlet  works 
capacities  can  be  included.  Other  special  operating  techniques  can  be  added. 
However,  these  future  considerations  do  not  in  any  way  detract  from  the  mean¬ 
ingful  and  highly  useful  results  given  in  this  report. 
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I.  INTRODUCTION 


AUTHORITY 

This  study  is  authorized  by  Section  2?  of  the  Hater  Resources 
Developne.;t  Act  (R.L.  93-251)  and  Section  214  of  the  Rivers  and  Harbors 
Act  of  1965  (R.L.  89-290).  The  study  was  requested  by  the  Pennsylvania 
Department  of  Environnental  Resources  (PER)  and  the  State  of  Mew  York 
Department  of  Environnental  Conservation  (DEC).  The  Philadelphia  District 
office  of  the  Corps  of  Engineers  is  responsible  for  the  conduct  of  the  work 

and  the  study  is  being  directed  by  a  group  consisting  of  staff  from  the 
Corps  of  Engineers,  the  Pennsylvania  HER,  the  New  vork  DEC,  the  Delaware 
River  Basin  Commission  (DPPC),  and  the  United  States  Ceological  curvey 
(USGS).  The  Philadelphia  District  Office  contracted  with  Water 
Resources  Engineers  (EIRE)  an  operatina  unit  of  Cam p  Dresser  &  McKee  Inc.  and 
the  firm  is  performing  work  under  contract  number  DACW61-78-C-0127. 

PURPOSE 

This  study  is  being  performed  in  two  najor  phases.  The  purpose 
of  phase  one,  which  is  documented  in  this  report,  is  to  investigate  the 
flow  characteristics  of  the  Delaware  River  Basin  under  naturalized  flow 
conditions  for  a  50-year  base  period  and  to  determine  the  effects  on 
these  flow  characteristics  caused  by  the  operation  of  the  Cannons vi lie, 
Pepacton  and  Me vers ink  Reservoirs.  These  reservoirs  are  a  source  of  water 
supply  to  the  City  of  New  York  and  are 'required  to  be  operated  to  meet 
specific  flow  objectives  at  Montague,  New  Jersey. 

The  purpose  of  phase  two,  is  to  determine  the  capabilities  of 
various  combinations  of  existing  and  proposed  reservoir  projects  to 
meet  maximum  sustained  low  flow  objectives  at  certain  specified  key 
locations. 


DESCRIPTION  OF  STUDY  AREA 
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The  study  area  encompasses  the  Delaware  River  Basin  from  its 
headwaters  to  the  Delaware  Memorial  Bridge  at  Wilmington,  Delaware.  The 
basin  drains  parts  of  southern  New  York  State,  eastern  Pennsylvania, 
western  New  Jersey  and  northern  Delaware.  The  drainage  area  of  the  basin 
above  Wilmington  is  11,030  square  miles  and  the  basin  averages  about  55 
miles  in  width  and  200  miles  in  length.  Elevations  range  from  near  2,600 
feet  above  mean  sea  level  in  the  headwaters  to  10  feet  at  Trenton,  New 
Jersey.  The  basin  is  divided  by  the  Blue-Kittatinny  Mountain  range  which 
runs  northeasterly  through  the  middle  of  the  basin.  The  upper  basin 
consists  of  the  eastern  slopes  of  the  Poconos  and  the  western  slopes  of 
the  Catskills  and  adjacent  mountains.  In  the  lower  portion  of  the  basin 
rolling  hills  gradually  diminish  to  the  low  areas  in  the  Atlantic  Coastal 
Plain.  Figure  1-1  outlines  the  total  Delaware  River  basin  and  shows  the 
major  tributaries  and  existing  and  proposed  reservoirs.  The  daily  flow 
model  nodes  and  the  key  locations  for  duration  and  frequency  analysis  are 
also  shown  on  the  figure. 

STUDY  PROCEDURES 

The  major  study  effort  involves  the  development  of  a  historical 
daily  flow  simulation  model  of  the  Delaware  River  Basin.  This  model  is 
used  to  produce  naturalized  daily  flows  and  to  analyze  the  effects  of  the 
operation  of  the  three  New  York  City  reservoirs. 

Natural  and  regulated  daily  flows  are  simulated  at  key  locations 
(model  nodes)  for  a  50-year  historical  period  from  October  1927  to  September 
1977.  The  first  step  in  the  modelling  effort  is  to  generate  50  years  of 
natural  inflows  to  each  model  node.  Model  nodes  are  put  at  existing 
USGS  gages,  points  of  proposed  or  existing  diversion,  places  of  proposed 
or  existing  impoundments,  or  other  points  of  interest.  Most  model  nodes 
are  USGS  stream  gaging  stations.  Natural  inflows  are  determined  by 
subtracting  historical  USGS  daily  flow  records  of  adjacent  stations. 

When  records  are  not  available  for  the  entire  50-year  period,  the  daily 
flows  of  a  station  are  generated  with  a  correlation  and  extension 
procedure.  Historical  records  are  extended  using  linear  regression 
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equations  to  long  term  recording  stations.  Natural  inflows  to  some 
headwater  stations  are  generated  with  a  correlation  and  extension  pro¬ 
cedure  even  though  records  are  available  because  the  record  had  become 
regulated.  The  correlation,  regression  analysis  and  record  extensions 
are  discussed  in  Chapter  II. 

With  the  use  of  the  naturalized  inflows  and  travel  times  which 
are  developed  for  various  reaches  throughout  the  basin,  the  simulation 
model  produces  50  years  of  natural  flows  at  60  key  locations  in  the 
basin.  The  daily  flow  model  and  its  development  is  described  in  Chapter 
III.  In  order  to  analyze  the  natural  flow  characteristics  in  the  basin, 
daily  flow  duration  curves  and  low  flow  frequency  curves  are  developed 
for  44  key  lcoations.  These  results  are  also  discussed  in  Chapter  III, 
and  are  presented  in  Appendix  A,  Naturalized  Daily  Flows. 

A  daily  regulated  flow  model  is  developed  for  the  Delaware  River 
Basin  to  determine  the  capability  of  the  three  New  Yrok  City  reservoirs 
to  meet  the  daily  flow  requirement  of  1,750  cfs  at  Montague,  New  Jersey. 
This  model  uses  the  natural  daily  inflows  developed  for  the  natural  flow 
model  and  the  release  criteria  specified  in  the  DRBC  Docket  D77-201  for 
the  reservoirs.  The  analysis  assumes  that  the  reservoirs  are  initially 
full  and  that  the  reservoirs  are  in  operation  for  the  entire  50-year 
period.  The  New  York  City  operating  procedures,  the  model  simulation  of 
the  reservoir  systems,  and  the  results  are  discussed  in  Chapter  IV. 

In  order  to  analyze  the  flow  characteristics  of  the  basin  simu¬ 
lated  by  the  operation  of  the  reservoirs  to  meet  the  flow  objectives  at 
Montague,  duration  and  low  flow  frequency  curves  are  developed  for  18 
key  locations.  These  results  are  presented  in  Appendix  B,  Regulated 
Daily  Flows. 


^Task  Group  Report  DRBC  Docket  No.  D-77-20 

Appraisal  of  Upper  Basin  Reservoir  Systems,  Drought 
Emergency  Criteria  and  Conservation  Measures,  Delaware 
River  Basin  Commission,  March  1979. 
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II.  NATURALIZATION  OF  MEAN  DAILY  FLOW  DATA 


INTRODUCTION 

Natural  daily  inflows  for  the  Delaware  River  Basin  above  the 
Delaware  Memorial  Bridge,  located  at  Wilmington,  Delaware,  are  developed 
at  key  locations  for  a  base  period  of  50  years.  The  base  period, 

October  1927  to  September  1977,  covers  water  year  1923  through  water 
year  1977. 

Existing  USGS  surface  water  stations  are  used  to  develop  the 
natural  inflows  for  the  base  period.  Some  of  these  stations  have 
complete  records  for  the  50  year  period;  others  have  only  partial 
records.  Some  stations  are  affected  by  regulation  from  upstream 
reservoirs  for  all  or  part  of  the  base  period  of  record.  Therefore,  the 
development  of  natural  inflows  for  the  entire  base  period  involves  an 
analysis  of  the  existing  records  and  periods  during  which  the  records 
are  affected  by  regulation. 

Basin  Network 

A  schematic  line  diagram  of  the  Delaware  River  Basin  from  its 
headwaters  to  the  Delaware  Memorial  Bridge  is  shown  in  Figure  1 1 - 1 . 

This  figure  incorporates  all  the  key  locations  for  analysis  required  in 
the  project.  The  diagram  shows  the  location,  number  and  name  of  the 
USDS  gaging  stations  used  in  the  development  of  the  natural  flows  for 
the  basin.  Reservoir  sites  and  names  are  also  shown.  The  USGS  gaging 
stations  shown  on  Figure  II-l  are  listed  in  Table  il-l  from  furthermost 
upstream  to  downstream.  The  table  includes  the  USGS  number,  location, 
period  of  record,  earliest  date  of  regulation  by  upstream  reservoir(s) 
and  the  drainage  area.  The  list  of  stations  also  includes  locations 
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needed  to  develop  natural  inflows  to  prospective  projects  on  the  Musconetcong, 
Lehigh  and  Schuylkill  Rivers.  The  period  of  record  indicated  in  the  table 
only  includes  the  base  period  from  1  October  1927  to  30  September  1977.  The 
periods  of  record  given  on  Table  1 1 - 1  are  those  for  which  daily  flow  data 
were  available  from  National  Water  Data  Storage  and  Retrieval  System 
(WATSTORE)  of  the  U.S.  Geological  Survey,  Reston,  Virginia.  These  periods 
of  records  may  vary  slightly  from  the  periods  of  record  given  in  the 
USGS  Surface  Water  Records  publications. 

The  reservoirs  shown  on  Figure  I I - 1  are  grouped  into  two  major 
categories:  1)  those  reservoirs  whose  effects  of  regulation  are  to  be 
removed  from  the  river  system  for  which  natural  inflows  are  developed 
and  2)  those  whose  releases  are  considered  to  be  a  part  of  the  natural 
system  for  the  purposes  of  this  study.  For  the  second  case  the  recorded 
flows  at  the  immediate  downstream  gaging  stations  are  used  without  any 
adjustments.  The  effects  of  the  Mongaup  Reservoir  System  and  Wal lenpaupack 
Reservoir  were  incorporated  into  the  natural  flow  model.  The  reservoirs 
are  operated  for  hydroelectric  power,  not  for  water  supply  or  flood  control, 
and  are  therefore  assumed  to  be  part  of  the  natural  system.  Table  1 1 -2 
is  a  list  of  the  existing  reservoirs  for  which  natural  inflows  for  the  50- 
year  base  period  are  developed.  The  date  each  reservoir  became  operational 
is  also  given.  These  dates  are  also  shown  on  the  line  diagram  of  Figure 
1 1 - 1  and  are  located  within  the  reservoir  triangular  symbol.  Natural 
reservoir  systems  are  indicated  by  an  "N." 

GENERAL  PROCEDURES  FOR  CORRELATION  AND  REGRESSION  ANALYSIS 

To  develop  the  natural  flows  at  USGS  gaging  station  sites, 
specific  procedures  for  correlation  and  extension  of  flow  records  are 
carried  out  for  several  difference  cases.  They  include: 

•  Extension  of  flows  for  all  short  term  stations  during  the 
natural  non-regulated  period  of  1  October  1927  to  2  June 
1953  based  on  long  term  non-regulated  stations, 

•  Extension  of  flows  for  all  non-regulated  short  term  stations 
for  the  remaining  portion  of  the  50-year  base  period  based 
on  non-regulated  long  term  stations. 


EXISTING  RESERVOIRS  IN  THE 
DELAWARE  RIVER  BASIN 
FOR  WHICH  NATURAL  INFLOWS  ARE  DEVELOPED 


•  Development  of  natural  inflows  at  existing  reservoir  sites 
for  50-year  base  period,  and 

•  Extension  of  several  short  tern  regulated  stations  based  on 
adjacent  long  term  regulated  stations  to  obtain  natural 
incremental  flow  data. 

This  section  summarizes  the  general  procedures  used  for  data  management, 
correlation  and  extension  of  records  for  all  cases.  The  development  of 
natural  inflows  to  the  existing  reservoirs  and  other  special  consid¬ 
erations  are  discussed  in  detail  in  subsequent  sections. 

Data  Acquisition  and  Data  Management 

The  mean  daily  stream  flow  data  from  1  October  1927  through  30 
September  1977  for  the  stations  listed  in  Table  1 1 - 1  were  obtained  on 
magnetic  tape  from  the  WATSTORE  System  of  the  USGS  in  Reston,  Virginia. 

A  computer  program,  written  by  WRE  called  FIND  is  used  to  produce  an 
output  file  of  daily  flows  which  are  needed  as  input  for  the  correla¬ 
tion  and  regression  analysis. 

In  general,  daily  flow  values  for  a  short  term  station  whose 
record  is  to  be  extended  are  retrieved  from  the  daily  flow  tape  with  the 
daily  flov/s  for  nearby  long  term  stations.  The  period  of  record  re¬ 
trieved  for  each  station  is  identical.  This  permits  a  direct  cor¬ 
relation  of  the  short  term  station  to  each  of  the  long  term  stations  for 
concurrent  periods  of  daily  flow.  The  six  major  long  term  unregulated 
stations  used  for  correlation  analysis  are: 

t  01469500,  Little  Schuylkill  River  at  Tamaqua,  Pennsyl vania , 

•  0147200,  Schuylkill  River  at  Pottstown,  Pennsylvania, 

•  01743000,  Perk i omen  Creek  at  Oraterford,  Pennsyl vania, 

t  0144R000,  Lehigh  River  at  Tannery,  Pennsylvania, 

•  01430500,  (old  #23^),  Lackawaxen  River  at  West  Hawley, 
Pennsyl vania,  and 

•  01431500,  Lackawaxen  River  at  Hawley,  Pennsylvania. 
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Other  long  term  stations  which  are  regulated  for  only  part  of 
their  period  of  record  are  also  used  for  specific  cases  (such  as 
determining  natural  inflow  at  reservoir  sites)  where  these  stations  are 
more  appropriate  because  of  the  correlation  period  needed  and  the 
proximity  of  these  long  term  stations  to  the  stations  whose  flow  records 
must  be  extended. 

Correlation  and  Regression  Analysis 

Correlations  and  the  determinations  of  linear  regression  equations 
are  performed  by  the  stepwise  regression  program  (BMD02R)  of  the  Bio¬ 
medical  Computer  Programs  developed  at  the  University  of  California,  Los 
Angeles.*  In  general,  the  daily  flow  data  input  to  the  program  consists 
of  one  short  term  station,  assigned  as  the  dependent  variable,  and  selected 
long  term  stations,  assigned  as  the  independent  variables.  Prior  to 
performing  the  regression  analysis  the  program  computes  the  means,  standard 
deviations,  and  the  correlation  matrix.  The  program  then  computes  in  a 
stepwise  manner  a  series  of  multiple  linear  regression  equations  for 
each  independent  variable  added  to  the  analysis.  The  first  independent 
variable  (i.e.  long  term  stations)  chosen  to  compute  the  regression  equa¬ 
tion  is  the  one  which  has  the  highest  correlation  coefficient.  Subsequent 
long  term  stations  are  added  to  the  regression  equation  one  at  a  time  in 
order  to  determine  if  the  predictability  of  the  correlation  equation  can 
be  improved.  At  each  step  the  program  computes,  among  other  statistical 
parameters,  the  multiple  correlation  coefficient,  standard  error  of  the 
estimate,  and  the  linear  regression  equation's  constant  and  coefficients. 

In  this  study  the  addition  of  more  than  one  independent  variable 
(long  term  station)  to  the  stepwise  regression  analysis  of  daily  flows 
did  not  produce  a  significant  increase  in  the  correlation  coefficient  or 
a  significant  reduction  in  the  standard  error  of  the  estimate  from  those 
computed  by  using  only  the  one  long-term  station  which  initially  had  the 

*BMD,  Biomedical  Computer  Programs,  Health  Sciences  Computing  Facility 
Department  of  Biomathematics,  School  of  Medicine,  University  of  Cali¬ 
fornia,  Los  Angeles,  University  of  California  Press,  January  1,  1973. 


best  correlation.  Therefore,  all  the  equations  developed  to  extend  the 
short-term  stations  are  based  on  the  simple  linear  regression  to  one 
station  only. 

At  the  beginning  of  the  project,  regression  analyses  were  per¬ 
formed  on  the  actual  flow  data  and  the  log  transform  of  the  flow  data. 

In  all  cases,  the  log  transform,  which  reduce  the  flow  extremities  to  a 
more  compact  scale,  produced  better  correlations  than  the  actual  data. 
Thus,  in  all  analyses,  the  log  transform  of  the  data  is  used  to  develop 
the  correlations  and  regression  equations. 

The  final  linear  regression  equation  for  the  extension  of  short¬ 
term  stations  as  produced  by  the  BMD  program’s  computation  of  the  con¬ 
stant  and  coefficient  is  in  the  form: 

L0G1QY  =  A  +  BL0G1QX 

where 

Y  =  flow  of  short-term  station, 

A  =  constant, 

B  =  coefficient,  and 
X  =  flow  of  long-term  station. 

This  equation  represents  the  relationship  of  the  short-term  station  to  a 
long-term  station  for  a  concurrent  period  of  flow  records.  It  is  then 
used  to  extend  or  fill  in  missing  periods  of  flow  record  for  short-term 
stations  and  to  develop  natural  flows  for  other  stations  during  periods 
of  regulation. 

Flow  Record  Extension 

The  extension  of  flow  records  of  the  short-term  stations  require 
the  manipulation  of  the  daily  flow  data  tape  and  the  use  of  the  regression 
equations  to  compute  the  extended  flows.  These  procedures  are  done  within 
a  program  written  by  WRE  called  FILLIN.  The  linear  regression  equations 
are  used  to  compute  the  daily  flows  of  missing  periods  of  record  for  the 
short-term  stations  or  the  natural  inflows  to  reservoirs.  The  extended 
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records  or  natural  inflows  to  reservoirs  are  combined  with  the  actual 
records  to  produce  the  complete  50-year  base  period  of  daily  flows. 

Some  stations  are  correlated  to  more  than  one  station;  the  different 
correlations  are  used  to  extend  different  parts  of  the  missing  or 
regulated  period  of  record. 

Table  1 1 -3  presents  a  summary  of  the  correlation  and  regression 
analyses  and  flow  record  extension  in  the  Delaware  River  Basin.  For 
each  station  listed  the  USGS  number  and  location  are  given.  The  correla¬ 
tion  or  correlations  which  are  finally  used  to  extend  a  station's  record 
are  given.  The  table  indicates  if  a  station's  record  is  not  extended. 

Some  station's  records  are  not  extended  because  these  stations  are  not 
used  as  model  nodes.  Other  station's  records  are  not  extended  because 
the  period  of  record  of  these  stations  exists  only  for  a  recent  regulated 
period,  making  it  impossible  to  extend  the  record  into  the  unregulated 
period  prior  to  1953.  For  each  correlated  station,  the  table  gives  the  cor¬ 
relation  period,  correlation  coefficient,  standard  error  of  the  estimate, 
and  the  constant  and  coefficient  of  the  linear  regression  equation  for  the 
log  transformation  of  the  observed  flow  data.  The  last  column  of  the 
table  indicates  the  period  of  record  which  is  extended  for  each  station. 

DEVELOPMENT  OF  NATURAL  INFLOWS  TO  EXISTING  RESERVOIRS 

The  development  of  the  natural  inflows  to  the  existing  reser¬ 
voirs  listed  in  Table  1 1 -2  are  done  by  correlating  the  reservoir  station 
to  specific  USGS  stations  for  unregulated  periods.  The  reservoir  station 
records  are  extended  for  any  missing  unregulated  periods,  and  natural 
inflows  are  developed  for  any  regulated  periods.  Descriptions  of  the 
correlation  and  regression  analysis  procedures  for  each  reservoir  are 
given  below. 

Pepacton  Reservoir 

The  Pepacton  Reservoir,  located  on  the  East  Branch  Delaware 
River  near  Downsville,  began  operating  in  September  1954.  The  USGS 
station  at  Downsville  (No.  01417000)  is  used  to  develop  the  natural 
inflows  to  the  reservoir.  For  the  50-year  base  period,  this  station's 
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TABLE  1 1-3.  CORRELATION  AND  FLOW  RECORD  EXTENSION  SUMMARY 

(continued) 
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Also  correlated  to  all  other  appropriate  predesignated  long-term  stations. 

^Used  to  determine  inflow  to  F.E.  Walter  Reservoir  for  period  2/1/61-9/30/77  with  drainage  area  adjustments. 


TABLE  I 1-3.  CORRELATION  AND  FLOW  RECORD  EXTENSION  SUMMARY 

(continued) 
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Also  correlated  to  all  other  appropriate  predesignated  long-term  stations. 


TABLE  11-3.  CORRELATION  AND  FLOW  RECORD  EXTENSION  SUMMARY 

(continued) 
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Also  correlated  to  all  other  appropriate  predesignated  long-term  stations. 
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Also  correlated  to  all  other  appropriate  predesignated  long-term  stations. 
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Flows  of  two  stations  added  together. 


TABLE  I 1-3 .  CORRELATION  AND  FLOW  RECORD  EXTENSION  SUMMARY 

(continued) 
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period  of  record  is  from  October  1 P4 1  to  September  1977.  For  the 
nonregulated  period  from  1927  to  1941,  the  Downsville  records  are 
extended  by  correlation  of  the  Downsville  station  to  the  station  at 
Fishs  Fddy  (No.  01421000)  on  the  East  Branch  Delaware  River.  Natural 
inflov/s  for  the  regulated  period  of  record  from  1954  to  1977  are  created 
by  correlation  of  the  Downsville  station  to  the  station  at 
Margerettevi 1 le  (No.  01413500),  which  is  located  above  the  reservoir  on 
the  East  Branch  Delaware  River. 

Cannonsville  Reservoir 

The  Cannonsville  Reservoir,  which  began  operation  in  October 
1963  is  located  on  the  West  Branch  Delaware  River  near  Stilesville.  The 
USDS  station  at  Stilesville  (No.  01425000)  is  used  to  develop  the 
natural  inflov/s  to  the  reservoir.  Its  period  of  record  is  from  October 
1952  to  September  1977.  For  the  non-regulated  period  from  1°27  to  lq52, 
the  Stilesville  records  are  extended  by  correlation  to  the  Hale  Fddy 
Station  (No.  01426500),  located  downstream  on  the  West  Branch  Pel  av/are 
River.  Natural  inflov/s  for  the  regulated  period  of  record  at 
Stilesville  are  determined  by  correlation  to  the  station  at  Walton  (No. 
D1423000),  located  above  the  reservoir  on  the  West  Branch. 

Prompton  Reservoir 

The  Prompton  Reservoir,  which  is  located  on  the  West  Branch 
Lackawaxen  River,  became  operational  in  July  196D.  The  USDS  station  at 
Prompton  (No.  D142900D)  is  used  to  develop  the  natural  inflov/s.  The 
period  of  record  for  this  station  is  from  October  1940  to  September 
1977.  All  long-term  stations  are  used  to  find  the  best  correlation  with 
Prompton.  The  downstream  station  at  Hawley  (No.  91431500)  is  selected 
to  extend  the  non-regulated  period  of  record  from  1927  to  1940  because 
it  has  the  best  correlation.  There  does  not  exist  a  gaging  station 
above  Prompton  Reservoir  which  could  be  used  to  develop  the  natural 
inflows  to  the  reservoir  during  the  regulated  period  by  correlation  with 
the  Prompton  station,  as  is  done  for  Pepact.on  and  Cannonsville 
Reservoirs.  Non-regulated  stations  in  the  upper  portion  of  the  basin 
are  chosen  for  correlation  analysis.  The  station  at  Margerettevil le 
(Mo.  0141 35DOO)  is  selected  for  correlation  and  is  used  to  extend  the 
flow  record  at  the  Prompton  site  during  the  regulated  period  of  record. 


The  extended  flows  at  Prompton  based  on  Margeretteville  are 
compared  for  verification  to  the  mean  daily  inflows  of  the  Monthly 
Reservoir  Operation  reports  prepared  by  the  Philadelphia  District  office 
of  the  Corps  of  Engineers.  The  reports  give  the  daily  inflow,  outflow, 
and  pool  elevation.  Records  are  available  from  October  1965  with  some 
missing  periods.  Non-regulated  flows  at  Prompton  based  on  Margeretteville 
compare  favorably,  with  some  scatter,  to  the  inflows  of  the  operation 
reports.  Therefore,  the  method  of  correlation  and  extension  with  Margeret¬ 
teville  is  considered  to  produce  an  adequate  representation  of  the  natural 
inflow  to  the  reservoir. 


Jadwin  Reservoir 

Jadwin  Reservoir,  which  became  operational  in  October  1°59,  is 
located  on  Pyberry  Creek,  a  tributary  to  the  Lackawaxen  River.  The  USDS 
station  located  on  Dyberry  Creek  near  Honesdale  (No.  R14?ogoo)  is  used 
to  develop  the  natural  inflows  to  the  reservoir.  Its  period  of  record 
is  from  Octoher  1943  to  Septenher  1977.  Flow  records  are  extended  by 
correlation  to  the  Lackawaxen  River  Station  at  Hawley  (No.  01431500)  for 
the  non-regul ated  period  and  by  correlation  to  Margerettevi 1 le  (No. 
01413500)  for  the  determination  of  the  natural  inflows  to  Jadwin  during 
the  regulated  period  of  record. 

Neversink  Reservoir 

fevers  ink  Reservoir,  which  is  located  on  Neversink  River, 
became  operational  in  June  1953.  The  USDS  station  at  Neversink  (No. 
01436000)  is  used  to  develop  the  natural  inflows  to  the  reservoir.  The 
period  or  record  during  the  50-year  base  period  is  from  October  1941  to 
Septenher  1977.  The  Neversink  flow  record  is  extended  by  correlation  to 
the  Oakland  Valley  Station  (No.  01437000)  on  the  Neversink  River  for  the 
non-regulated  period  from  1927  to  1941.  The  creation  of  the  natural 
inflows  to  the  reservoir  during  the  regulated  period  of  the  Neversink 
station  from  1953  to  1977,  is  done  by  correlation  to  the  Claryville 
Station  (No.  01435000)  which  is  located  above  the  Neversink  Reservoir. 
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F.  E.  Walter  Reservoir 


F.  E.  Walter  Reservoir,  which  became  operational  in  February  1961, 
is  located  on  the  Lehigh  River  near  White  Haven.  The  US6S  station  near 
White  Haven  (No.  014478000)  is  used  to  develop  the  natural  inflows  to  the 
reservoir  during  the  non-regulated  period.  Its  period  of  record  is  from 
October  1957  to  September  1977.  The  White  Haven  flow  record  is  extended 
by  correlation  to  the  station  on  the  Lehigh  River  at  Tannery  (No.  01448000) 
for  the  period  from  1927  to  1957.  The  natural  inflows  to  the  reservoir  from 
1961  to  1977  are  developed  by  prorating  the  sum  of  the  flows  at  two  upstream 
stations,  Stoddartsvi lie  on  the  Lehigh  River  and  Blakeslee  on  Tobyhanna 
Creek,  to  the  damsite  by  a  drainage  area  proportion  ratio. 

The  natural  inflow  of  the  reservoir  based  on  a  drainage  area  ratio 
and  the  two  upstream  stations  are  compared  to  the  mean  daily  inflows  of  the 
Monthly  Reservoir  Operation  reports  prepared  by  the  Philadelphia  District 
office  of  the  Corps  of  Engineers.  Records  are  available  from  February  1973 
through  September  1977  with  some  missing  periods.  The  flows  developed 
from  the  two  upstream  stations  compare  very  well  to  the  flows  of  the 
operation  reports  and  this  method  is  used  to  determine  the  natural  inflows 
to  the  reservoir  during  the  regulated  period  of  the  White  Haven  station. 

Beltzville  Lake 

Beltzville  Lake,  which  is  located  on  Pohopoco  Creek,  a  tributary 
to  the  Lehigh  River,  became  operational  in  February  1971.  The  USGS 
station  on  Pohopoco  Creek  below  Beltzville  dam  site  (No.  01449800)  is 
used  to  develop  the  natural  inflows  to  the  lake.  The  period  of  record 
at  this  station  is  from  1967  to  1977  for  the  50-year  base  period.  For 
the  non-regulated  period,  the  flow  record  at  tne  Dam  Site  is  extended, 
based  on  correlation  to  the  station  at  Parryville  (No.  01450000)  whose 
records  are  extended  based  on  correlation  to  the  Lehigh  River  station 
at  Tannery  (No.  01448000).  The  creation  of  the  natural  inflows  to  the 
lake  at  the  Dam  Site  for  the  regulated  period  of  record  is  based  on 
the  correlation  of  the  Dam  Site  station  to  the  station  at  Kresgeville 
(No.  01449360)  which  is  located  upstream  of  the  reservoir  and  has  a 
period  of  record  from  1966  to  1977. 


The  extended  flows  at  the  Dam  Site  based  on  Kresgeville  are 
compared  for  verification  to  the  mean  daily  inflows  of  the  Monthly 
Reservoir  Operation  reports  prepared  by  the  Philadelphia  District 
Office  of  the  Corps  of  Engineers.  Records  are  available  from  February 
1973  through  September  1977  with  some  missing  periods.  The  natural  inflow 
at  the  Dam  Site  based  on  correlation  with  Kresgeville,  compare  favorably 
to  the  inflows  of  the  operation  reports.  Therefore,  the  method  of  corre¬ 
lation  and  extension  with  Kresgeville  is  considered  to  produce  an  adequate 
representation  of  the  natural  inflow  to  the  Lake  for  the  regulated 
period  of  the  Dam  Site  station. 

Nockamixon  Reservoir 

The  Nockamixon  Reservoir,  which  became  operational  in  December 
1973,  is  located  on  Tohickon  Creek  near  Pipersville.  Natural  flows  are 
developed  at  the  USGS  station  at  Pipersville,  whose  period  of  record  is 
from  October  1935  to  September  1977.  The  non-regulated  period  of  record 
at  Pipersville  is  correlated  to  the  long-term  stations.  The  long-term 
station  which  produces  the  best  correlation  is  at  Graterford  on 
Perkiomen  Creek  (No.  01473000).  This  station  is  used  to  extend  the  flow 
record  of  Pipersville  from  1927  to  1935  for  the  non-regulated  period  and 
also  from  1973  to  1977  for  the  regulated  period  of  the  Pipersville 
station. 

SPECIAL  CONSIDERATIONS 

In  addition  to  the  general  procedures  described  for  correlation, 
regression  analysis,  extension  of  flow  records,  and  the  development  of 
the  natural  inflows  to  the  existing  reservoirs,  special  considerations 
are  required  for  the  development  of  natural  flows  at  the  Philadelphia 
station  on  the  Schuylkill  River  and  for  the  Delaware  and  Raritan  Canal 
diversions. 
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Schuylkill  River  Station  at  Philadelphia 


The  recorded  daily  flows  at  Philadelphia  do  not  include  the 
diversions  above  the  station  by  the  City  of  Philadelphia  for  municipal 
water  supply.  The  diversion  flow  is  accounted  for  in  order  to  produce 
the  natural  flow  at  the  Philadelphia  station  for  its  period  of  record 
from  October  1931  to  September  1977,  and  to  correlate  the  paged  flows 
plus  the  diversion  flows  to  a  long-term  station  in  order  to  extend  the 
flow  record  from  1927  to  1931. 

Pumping  records  were  ohtained  from  the  City  of  Philadelphia  for 
the  Belmont,  'lueen  Lane,  and  Shawriont  Schuylkill  Piver  stations.  These 
records  and  the  gaging  station  records  are  used  to  create  the  natural 
flow  of  the  Schuylkill  River  at  Philadelphia  from  1931  to  1977.  The 
extension  of  the  natural  flow  record  from  1927  to  1931  is  accomplished 
by  correlating  the  Philadelphia  station  flows  plus  the  diversion  flows 
to  the  sun  of  the  flows  recorded  at  the  Pottstown  station  on  the 
Schuylkill  River  (No.  01472(100)  and  *-he  Oraterford  station  on  Perkionen 
Creek  (No.  01473000). 

Oelaware  and  Raritan  Canal 

The  Delaware  and  Raritan  (Pf.R)  Canal  diverts  water  from  the 
Delaware  River  at  Raven  Rock,  which  is  located  north  of  Trenton,  New 
Jersey.  Flow  records  are  available  beginning  March  1947  to  the  present 
from  the  HSOS  gaging  station  located  at  Kingston,  New  Jersey. 

In  order  to  determine  the  diversion  flow  prior  to  1947  the  New 
Jersey  Bureau  of  Water  Facility  Operations,  were  contacted  for 
assistance.  No  flow  records  for  this  period  are  available,  and  the 
canal  went  through  several  operational  and  rehabi  1  itat.ion  changes 
throughout  this  period.  The  canal,  which  was  oriqinally  open  for 
navigation  in  1B34,  was  closed  to  navigation  in  1°33  by  the  railroad 
company  which  owned  the  canal  at  that  tine.  In  1934  ,  the  State  of 
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Jersey  became  the  new  ov/ner,  and  for  the  next  ten  years,  studied  the 
possible  uses  of  the  canal.  In  1944,  the  DAR  Canal  was  rehabil  itated 
for  use  as  an  industrial  water  supply,  and  in  1W9  it  was  rebabi  1  i tated 
for  potable  water  supply. 

An  estimate  of  the  canal  diversions  before  the  establishment  of 
the  gage  at  Kingston  was  based  on  an  interview  with  operational 
personnel.  For  the  ungaged  period  from  1927  to  1947,  the  daily 
diversion  flow  for  the  DRR  Canal  is  estimated  to  he  50  cfs. 

CORRELATION  ACCURACY  CHECK 

Additional  checks  on  the  correlation  results  are  conducted  for 
eight  selected  stations.  They  include: 

0  01417000,  East  Branch  Delaware  River  at  Downsville,  NY, 

•  01425000,  West  Branch  Delaware  River  at  Stilesville,  NY, 

•  01428500,  Delaware  River  near  Barryville,  NY, 

0  01436000,  Neversink  River  at  Neversink,  NY, 

•  01437000,  Neversink  River  at  Oakland  Valley,  Nv, 

0  01438500,  Delaware  River  at  Montague,  NY, 

0  01451000,  Lehigh  River  at  Walnutport,  PA,  and 

0  01474500,  Schuylkill  River  at  Philadelphia,  PA. 

Concurrent  years  of  flow  for  these  stations  and  the  stations  used  in  the 
correlation  analysis  are  used  with  the  correlation  equations,  presented 
in  Table  II-3,  to  re-qenerate  the  actual  flows  for  the  eight  stations. 

The  uownsville  and  Stilesville  gages  actually  use  two  correlation  equations 
to  complete  the  record.  However,  for  the  correlation  check  only  the  equa¬ 
tions  developed  by  using  Fishs  Eddy  for  Downsville  and  Hale  Eddy  for 
Stilesville  are  used. 

The  accuracy  of  the  re-generated  flows  are  analyzed  by  correlating 
the  re-generated  flows  to  the  actual  flows  using  the  regression  program 
BMD02R  discussed  and  referenced  in  the  section  Correlation  and  Regression 
Analysis.  The  results  of  the  correlation  of  observed  and  confuted  flows 
for  the  eight  stations  are  presented  in  Table  1 1 -4 .  The  table  gives  the 
station  number  and  location,  the  correlation  coefficient,  the  observed 
and  computed  mean  discharge  for  the  concurrent  years  of  flow  used,  and 
the  constants  A  and  B  for  the  linear  regression  equations.  The  equation 
is  in  the  form: 


Y  =  A  +  B  ( X ) 
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where 


Y  =  Observed  flow  (cfs), 

X  =  Computed  flow  (cfs)  and 
A,B  =  constants 


Six  of  the  eight  stations  have  correlation  coefficients  greater  than  0.95. 
Lower  correlation  coefficients  are  calculated  for  the  Neversink  River, 

0.93  at  Neversink  and  0.87  at  Oakland  Valley.  The  differences  in  the 
observed  and  computed  mean  discharges  for  the  period  of  record  analyzed 
are  less  than  two  percent  for  all  stations  except  the  two  on  the  Neversink 
River.  For  the  station  at  Neversink  the  computed  mean  is  6  percent  low 
and  for  the  station  at  Oakland  Valley  the  computed  mean  is  16  percent  high 

The  equations  presented  in  Table  1 1-4  demonstrate  how  well  the 
re-generated  (or  computed)  flows  compare  to  the  observed  flows.  If  the 
computed  flows  were  re-generated  to  be  equal  to  the  observed  flows  for 
each  and  every  day  of  the  period  of  record  analyzed  the  constant  A  would 
equal  zero  and  the  constant  B  would  equal  one.  Graphically,  this  perfect 
fit  would  produce  a  line  with  the  y-intercept  at  zero  and  a  slope  of  one 
(45°  angle)  when  the  observed  versus  computed  values  are  plotted.  Figures 
1 1 - 2  through  1 1 -9  present  plots  of  the  regression  equation  developed  for 
the  concurrent  years  of  discharges  for  the  eight  stations.  In  each  figure 
the  observed  is  given  on  the  ordinate  and  the  computed  is  given  on  the 
abscissa.  The  regression  equation  is  plotted  with  triangular  symbols  to 
a  point  approximately  equal  to  twice  the  mean  flow.  The  figures  also 
include  a  straight  line  at  a  45°  angle  without  a  symbol  to  represent  a 
perfect  fit.  These  figures  demonstrate  that  the  correlation  equations 
used  to  fill-in  missing  or  regulated  periods  of  record  are  highly  accurate 
with  the  exceptions  of  the  Neversink  River  at  Neversink,  Figure  II-5, 
the  Neversink  River  at  Oakland  Valley,  Figure  II-6,  and  the  Lehigh  River 
at  Walnutport,  Figure  II-8. 

The  Neversink  River  at  Neversink  extended  flows  are  based  on  the 
Neversink  River  gage  at  Oakland  Valley  as  shown  in  Table  1 1 -3 .  For  the 
concurrent  period  of  record  analyzed,  the  re-generated  flows  for  Neversink 
from  Oakland  Valley  show  that  the  filled-in  flows  are  over  estimated  for 
discharges  less  than  150  cfs  and  are  underestimated  for  discharges  greater 
than  150  cfs.  Therefore,  during  low  flow  periods  of  the  actual  extended 


period  from  October  1927  through  September  1941  the  flows  at  Neversink 
are  slightly  overestimated.  For  example  the  7Q10  for  Neversink  based  on 
the  50  years  of  naturalized  flows,  which  includes  filled-in  and  recorded 
flows,  is  15  cfs.  The  7Q10  at  Neversink  based  on  the  regression  equation 
of  Table  I I -3  and  the  7Q10  at  Oakland  Valley  is  20  cfs. 

The  Neversink  River  gage  at  Oakland  Valley  is  extended  for  one 
year  only  (water  year  1928)  and  is  based  on  the  flows  at  West  Hawley  on 
the  Lackawaxen  River.  Figure  1 1 -6  shows  that  for  the  concurrent  period 
of  record  analyzed  the  re-generated  flows  are  always  overestimated.  Thus, 
the  filled-in  flows  for  the  water  year  1928  are  probably  greater  than  those 
that  actually  occurred  for  that  one  year. 

The  Lehigh  River  at  Walnutport  extended  flows  are  based  on  the 
Lehigh  River  flows  at  Bethlehem  as  shown  in  Table  1 1  - 3 .  For  the  concurrent 
period  of  record  analyzed  the  re-generated  flows  are  underestimated  for 
discharges  less  than  1700  cfs  and  are  over  estimated  for  discharges  greater 
than  1700  cfs  as  shown  in  Figure  II-8.  This  relationship  is  produced  for 
the  concurrent  period  of  record  from  water  year  1953  through  water  year 
1960.  However,  analyzing  the  50-year  of  naturalized  flows  the  7Q10  for 
Walnutport  is  205  cfs,  and  the  7Q1 0  at  Walnutport  based  on  the  7Q10  at 
3ethlehem  and  the  regression  equation  of  Table  1 1-3  is  222  cfs  which  is  an 
eight  percent  increase.  Therefore,  in  this  case,  the  analysis  of  the 
concurrent  period  of  record  may  not  be  indicative  of  the  relationship 
produced  by  the  regression  equation  developed  from  the  correlation  of 
Walnutport  to  Bethlehem.  Thus  the  regression  equation  is  maintained  to 
predict  the  missing  period  of  record  for  Walnutport. 

In  general  the  correlation  analysis  and  fill-in  procedures  are 
adequate  to  produce  50  years  of  daily  naturalized  flows  for  modeling  of 
the  Delaware  River  basin.  The  total  flows  developed  for  these  stations 
and  the  other  stations  in  the  basin  are  used  to  produce  incremental  inflows 
by  lagging  and  subtracting  the  filled-in  flow  records  of  adjacent  stations 
as  described  in  the  following  Chapter. 
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TABLE  I I -4 


CORRELATION  OF  OBSERVED 
AND  COMPUTED  (RE-GENERATED)  FLOWS 


Station 

Number  &  Location 

Correlation 

Coefficient 

Mean  Discharge1 
(cfs) 

Observed  Computed 

Equations 
Y  =  A 
A 

;  Constants2 
+  B  ( X ) 

B 

01417000 

E.  Br.  Del.  R.  at 
Downsville  NY 

0.979 

698 

693 

24 

0.972 

01425000 

W.  Br.  Del .  R.  at 
Stilesville  NY 

0.996 

811 

796 

-4 

1.024 

01428500 

Delaware  River  nr 
Barryville  NY 

0.972 

3783 

3721 

48 

1 .004 

01436000 

Neversink  R.  at 
Neversink  NY 

0.933 

246 

232 

-27 

1 .179 

01437000 

Neversink  R.  at 
Oakland  Valley  NY 

0.870 

434 

503 

-37 

0.936 

01438500 

Delaware  R.  at 
Montague  NY 

0.997 

6451 

6472 

117 

0.979 

01451000 

Lehigh  R.  at 
Walnutport  PA 

0.987 

1738 

1748 

224 

0.867 

01451000 

Schuylkill  R.  at 
Philadelphia  PA 

0.973 

2924 

2910 

33 

0.993 

Jfnean  discharge  for  correlation  period 
2y  =  observed  flow,  X  =  computed  flow 
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c?m<?omcocoo 


DOWNSVILLE,  N.Y. 


FIGURE  1 1-2.  Correlation  Accuracy  Check:  01417000 

East  Branch  Delaware  River  at  Downsville,  NY 
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oncxmwoo 


STILESVILLE,  N.Y. 


-200  200  600  1000  t400  1800 


COMPUTED 

FIGURE  1 1-3.  Correlation  Accuracy  Check:  01425000 

West  Branch  Delaware  River  at  Stilesville,  NY 
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BARRYVILLE,  N.Y. 


FIGURE  II- 


■1—  I  -  i  ■  1  I  L. 

000  4000  6000  8 

COMPUTED 


4.  Correlation  Accuracy  Check:  01428500 
Delaware  River  near  Barryville,  NY 


cmcximcoooo 


orncsomcoooo 


OAKLAND  VALLEY,  N.Y. 


1200 


-200  L 

0  200  400  600  800  1000  1200 

COMPUTED 


FIGURE  1 1 - 6 .  Correlation  Accuracy  Check:  01437000 
Neversink  River  at  Oakland  Valley,  NY 
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OOQCOUJQ£>UJO 


MONTAGUE,  N.J. 


FIGURE  1 1-7.  Correlation  Accuracy  Check:  01438500 
Delaware  River  at  Montague,  NO 
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OCDCOLUQ£>UJO 


PHILADELPHIA,  PA. 


COMPUTED 


FIGURE  1 1-9.  Correlation  Accuracy  Check:  01474500 
Schuylkill  River  at  Philadelphia,  PA 
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III.  NATURAL  FLOW  MODEL 


PURPOSE  or  THE  NATURAL  FLOW  MODEL 

The  purpose  of  the  natural  flow  model  is  to  negate  the  effect  of 
all  existing  reservoirs  (except  the  Wal lenpaupack  Reservoir  and  the  Mongaun 
Reservoir  system)  in  order  to  develop  natural  daily  flows  at  key  stations  in 
the  basin  for  the  50-year  base  period  of  record.  Natural  daily  flow  duration 
and  low  flow  frequency  analyses  are  performed  at  those  key  stations. 

The  natural  flov;  model  consists  of  a  branching  network  of  nodes 
as  shown  on  Figure  1 1 - 1 .  Most  nodes  are  also  USDS  surface  water 
stations.  The  model  is  driven  by  a  set  of  natural  incremental  inflows 
to  each  node.  These  incremental  inflows  are  conhined  with  a  simple 
linear  lag  function  frori  node  to  node  to  produce  natural  daily  flows. 

The  model  beoins  at  the  rrost  upstream  nodes  on  each  day.  The  natural 
inflows  to  thrse  headwater  nodes  is  lagged  to  downstream  nodes  while  the 
inflow  between  the  nodes  is  added,  producing  a  daily  natural  flow  at  the 
lower  nodes.  It  is  very  often  the  case  that  the  headwater  nodes  are 
existing  reservoirs. 

The  natural  inflows  to  the  reservoirs  are  taken  from  the  analyses 
described  in  Chapter  II.  The  incremental  inflows  at  the  interior  nodes 
are  determined  by  subtracting  the  routed  recorded  or  extended  flows  of 
an  upstream  node  from  a  downstream  node.  This  procedure  assumes  that 
the  natural  inflow  between  two  adjacent  nodes  is  independent  of  regula¬ 
tion.  These  natural  incremental  inflows  drive  the  natural  and  regulated 
flow  models.  The  natural  and  regulated  flow  models  assume  identical 
rainfal 1-runoff  as  that  of  the  50-year  base  period. 
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MODEL  OVERVIEW 


The  basic  architecture  of  the  natural  flow  model  is  complicated 
considerably  by  the  requirements  of  this  particular  project.  The  concepts 
of  the  model  are  easy  to  comprehend  but  the  extremely  large  data  set 
which  is  used  as  input  data  has  made  it  necessary  to  disagrega  the  various 
components  of  the  model  to  minimize  redundant  calculations.  To  this  end 
the  model  is  divided  into  several  completely  independent  programs.  Many 
of  these  programs  are  simply  processors  which  manipulate  the  data  to 
insure  efficient  and  less  costly  simulation.  Table  1 1 1  - 1  highlights  these 
processors  and  indicates  the  numerous  tapes  used  in  this  project.  The 
various  input  files,  programs  and  output  files  are  executed  in  the  order 
shown  in  the  table.  The  following  describes  the  various  procedures  fol¬ 
lowed  in  each  of  these  eight  programs. 

Program  SORT 

The  first  step  of  the  simulation  is  to  preprocess  the  input  data 
into  a  format  which  is  compatible  with  the  other  programs.  The  raw  data 
in  its  most  basic  form  is  daily  discharge  as  recorded  by  the  USGS.  This 
data  is  provided  on  a  card  image  magnetic  tape  using  the  format  depicted 
in  Figure  1 1 1  - 1 .  Four  cards  are  required  for  each  month  of  daily  flows. 

Each  card  has  the  station  number,  year  and  month.  Eight  days  of  flow  data 

are  given  on  each  of  the  four  cards  except  the  last  which  has  the  remaining 

25-31  day  flows.  This  tape  is  ordered  "by  station,  by  date,"  and  an  example 

is  shown  in  Figure  1 1 1 - 2 .  The  first  eight  digit  number  after  the  card  type 

identifier  "3"  is  the  station  number  (01479000)  which  is  followed  by  the 
year,  month,  card  number  (1  through  4),  and  the  daily  flows.  Card  type 
"2,"  shown  in  the  middle  of  the  figure  is  contained  in  the  USGS  data  at 
the  beginning  of  each  water  year.  It  contains  the  station  number  (01479000), 
the  parameter  code  (00060  for  stream  flow  in  cfs),  the  statistics  code 
(00003  for  mean  value  for  each  day)  and  "ENT." 

The  shear  volume  of  data  (93,218  records)  required  for  this  project 
makes  this  format  unacceptable.  For  this  reason  the  SORT  program  reorders 
the  data  set  into  a  "by  date,  by  station"  format,  as  shown  in  Figure  1 1 1 -3 . 


1 1 1 -2 


available  for  that  particular  day. 
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FIGURE  1 1 1  -  2 .  Example  of  Recorded  U.S.G.S.  Data  Ordered  by  Station  by  Date 
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The  first  eight  digit  numbers  after  the  card  type  "3"  are  the  station  numbers. 
This  example  shows  the  flows  of  days  1-8  (Card  No.  1)  for  February  1928 
(192802  in  the  figure)  at  several  stations.  They  are  followed  by  the  flows 
for  days  9-16  (Card  No.  2)  for  February  1928  for  several  stations.  This 
reformatted  tape  eliminated  the  necessity  of  rereading  the  entire  data  set 
for  each  fill-in  operation. 

The  heart  of  this  program  is  a  system  subroutine  called  SORTMERGE. 

The  advantage  of  the  system  subroutine  is  extreme  efficiency  relative  to 
procedures  using  conventional  FORTRAN  sort  routines.  Unfortunately  system 
subroutines  are  also  system  specific.  Thus  it  will  be  necessary  to  make 
some  slight  modifications  before  this  program  could  be  executed  on  a  system 
other  than  the  Control  Data  Corporation  ( CDC )  CYBER  176  system  used  in  this 
study . 

The  output  tape  generated  in  SORT  is  used  as  input  data  for 

FILL  IN. 

Program  FILLIN 

The  correlation  and  regression  analysis  used  to  both  extend  short 
term  stations  and  eliminate  the  effect  of  regulation  on  long  term  stations 
is  described  in  detail  in  Chapter  II.  FILLIN  is  the  program  which  applies 
the  resultant  correlation  equations  to  extend  the  existinq  USGS  record 
thereby  obtaining  a  complete  50-year  record  for  the  gaging  stations. 

Because  of  problems  with  data  acquisition  and  stations  with  special  consider¬ 
ations,  FILLIN  requires  some  special  input  tapes  and  files  to  complete  the 
period  of  record  of  all  stations.  The  New  York  City  Reservoir  daily  inflows 
are  also  written  onto  the  output  tape  to  consolidate  all  of  the  input  data. 
Figure  1 1 1 -4  is  an  example  of  the  filled  in  data  set.  The  outout  tape 
generated  by  this  program  is  of  the  same  general  format  as  the  original  USGS 
tape  except  that  it  is  ordered  by  date  by  station.  This  tape,  however, 
includes  50  years  of  records  for  all  stations.  The  records  of  those  stations 
which  are  not  required  for  further  analysis  are  filled  with  999999  where 
recorded  data  are  not  available.  The  output  tape  from  FILLIN  is  used  as 
input  data  for  INFLOW. 
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FIGURE  II 1-4.  Example  of  the  Filled-In  Data  Set 


Program  INFLOW 


This  program  is  perhaps  the  most  sensitive  and,  therefore,  most 
important  facet  of  the  entire  simulation  process.  The  output  tape  from 
this  processor  drives  both  the  natural  flow  model  and  the  regulated  flow 
model . 


The  actual  mechanics  of  this  program  are  discussed  in  the  subsection 
entitled  Incremental  Runoff  and  Routing  Procedures.  Incremental  flows  are 
calculated  by  lagging  and  subtracting  filled-in  flow  records  of  adjacent 
stations.  When  the  incremental  inflows  at  all  model  nodes  which  are  USGS 
stations  are  realistically  represented,  the  inflows  at  nodes  which  are  not 
USGS  stations  are  calculated  using  specific  prorating  procedures.  Inflows 
to  three  USGS  stations  whose  records  could  not  be  extended  (Callicoon, 

Tocks  Island,  and  Glendon)  are  also  calculated  using  specific  prorating 
procedures.  The  prorating  procedures  used  to  calculate  inflows  at  nodes 
which  are  not  USGS  stations  or  are  USGS  stations  whose  record  could  not 
be  extended  are  presented  on  Figure  1 1 1-5 .  For  example,  the  incremental 
inflow  to  Callicoon  is  calculated  by  subtracting  the  flows  at  Fishs  Eddy 
and  Hale  Eddy  from  the  flow  at  Barryville  and  multiplying  this  value  by  a 
drainage  area  ratio.  The  drainage  area  ratio  is  determined  by  subtracting 
the  drainage  areas  of  Fishs  Eddy  and  Hale  Eddy  from  the  drainage  area  of 
Callicoon  and  dividing  this  result  by  the  value  obtained  when  the  drainage 
areas  of  Fishs  Eddy  and  Hale  Eddy  are  subtracted  from  the  drainage  area 
of  Barryville. 

The  output  tape  supplies  the  input  data  to  both  the  natural  and 
regulated  models.  An  example  of  the  incremental  inflow  file  format  is 
shown  in  Figure  1 1 1-6. 

Program  NATFLO 

To  this  point  in  the  simulation  sequence,  all  of  the  programs 
have  been  data  preprocessors.  NATFLO,  on  the  other  hand,  is  the  natural 
flow  model.  It  takes  the  inflow  generated  in  INFLOW,  routes  them  through 
the  system  using  a  linear  lag  formulation  and  creates  the  output  in 
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several  formats.  These  various  output  formats  are  designed  to  facilitate 
any  further  analysis  which  may  be  required. 


The  daily  flows  for  each  month  are  printed  t  in  yearly  tables. 
The  tables  for  each  year  also  include  for  each  month,  the  minimum,  maximum, 
mean  and  total  daily  flows  and  the  standard  deviation.  Figure  1 1 1 -7  is  an 
example  printout  showing  the  1953  daily  discharges  for  the  Delaware  River 
at  Barryville.  As  a  backup,  an  unformatted  file  is  saved  on  magnetic  tape. 
This  tape  is  an  efficient  means  of  recalling  the  stations,  dates  and  flows 
to  be  reformatted  in  the  remaining  post-processors.  It  can  also  be  used  to 
create  multiple  copies  of  the  results  without  rerunning  the  entire  program. 

Programs  REFORM1/SORT/REFORM2 

These  remaining  three  programs  are  post-processors  which  reformat 
the  model  output  into  a  form  compatible  with  the  USGS  WATSTORE  Daily  Values 
statistics  program  A969  which  produces  duration  and  frequency  tables.  The 
first  program  converts  the  binary  file  to  standard  USGS  card  image  format 
ordered  "by  date,  by  station."  The  second  program  reorders  the  file  to  "by 
station,  by  date"  and  the  last  program  simply  inserts  the  annual  header 
cards  required  by  the  A969  program.  At  this  point,  the  model  output  is 
ready  for  processing  by  the  USGS  A969  program. 

INCREMENTAL  RUNOFF  AND  ROUTING  PROCEDURES 

As  previously  stated,  historical  incremental  runoff  (or  inflow) 
drives  the  natural  and  regulated  flow  models. 


The  entire  simulation  process  is  not  intended  to  he  a 
rainfal 1/rimoff  processor.  The  purpose  of  this  modelling  effort  is  to 
*ocus  on  the  low  flow  events  in  the  basin.  Painfal 1-runoff  processors 
by  their  formulation,  tend  to  misrepresent  small  rainfall  events  and 
groundwater  accretion.  This  being  the  case,  rather  than  simulating 
runoff  from  rainfall,  the  recorded  and  extended  USGS  record  is  used  to 
determine  the  historical  runoff  between  adjacent  gaging  stations.  This 


II 1-12 


r 


i 

i 


UJ 

o  a. 

o  CJ 


_3s 

U.  O 
H  U4 

-J  cc 
«*  <  «* 
a  -4  3 

d  r>< 

w.  a;  _j 

<  HU 
2  MO 


1  < 

:  j 

Q.  -J 

©  ! 

IJ 

1 

O 

Ui 

« 

©  I 

/v 

.  « 

<  j 

i 

1 

• 

1  • 

u>  j 

i 

1A  i 

«  « 

« 

1 

i  i 

°  i 

| 

p  ! 

-j  * 

1 

i 

* 

i 

a  ; 

; 

-/  i 

: 

<  i 

l 

j 

o 

j 

j 

1 

$ 

» 

c* 

i 

j 

# 

1 

f 

1 

1 

© 

9 

w 

j 

1 

u 

1 

7 

i 

I 

►*4 

i 

1 

© 

1 

z 

l 

1 

U' 

1 

1 

p 

1 

1 

©  *- 

i 

Ui  Z 

I  j 

Ci  T- 

i 

1 

or  o 

{ 

1 

3  « 

i 

1 

O  *4 

; 

9 

©  > 

4  i 

UJ 

1 

©  • 

m 

i 

1  j 

a 

p 

i 

I  | 

©  -j 

o 

i 

i 

| 

Ui  UJ 

CM 

^  M 

*  l 

t 

P  u 

' 

1 

3  © 

!  1 

« 

Z 

1 

1 

P  •- 

ui 

1 

r  © 

0 

1  | 

a  Cl. 

O 

i 

1 

u  © 

Z 

i 

i 

!  | 

i 

■J 

1  u  oi 
•  UJ  43 
.  O  CD 


*  O  (M  ! 

w  m  in  ( 


o>  o  o  o 
43  CO  p  43 
CD  ©  *0  -« 
f\l  CM  CM  fvi 

rj  cm  p  m 

^  o  «h  p 


O'  in  in 

m  cJ 

©  rnj 

CM  CJ 


*o  O' 

43  c- 


mt*o  i 
pjm  i 


.  o 
.  z 


Jofl'P' 
cd  p  p 
ok  <r  m 


©  P  i 
4)  cO  « 
•H  ©  < 


i/>  ro  p  cc 

p|  m  »-•  in 


■  PJ  If)  *0  CO 
I  CO  4)  P  CM 
»  cV  cm  cj  co 


*"1 
a 


*-  eo 
OJ  *n 
Cj  co 


4)  k  m  in  ■ 

♦|0  D  M 


co;  O'  »o,  m  i 
«f "O  m  (v« 


cof  co  *n  p 
pi  p  cop 


cc  «-•  f 
tc  cc  » 


cm.  -<  ©  »o 
co;  rv 

*Ol  CO  C0*0 


Z5  . 

:->  « 

I  « 


lo 


in  45  • 

*»#  *4  I 

*o  co  i 


©  cv  p  p1  p  ©1 
•n  lie  m  o  w 


n  m  * 
4)  o  r— 
CC  P'P 


O  CO  v© 

c-  Off 


fro  I 
cc  p  . 

VO  4)  ’ 


p  in  ' 
—i  co  < 


«•  ©’O' 
p  roitt 
op 
CJ  tVj-t 

I 

^  ©Ip 

-  e'|<r 


©  © 
O  4) 


p  un  i 
tv.cr 


o  cpi© 

O  P'10 

<v  «•>*► 

P  CJl© 


cc  ©  i 

P'  CO  l 


I 


©Im  < 

i  4)  for  ' 


♦  PIP 
rv  ©ivO 
•  •! 
o  ©ip  frier* 

<)  v»||T  ©  ,m 
c*  in  Ip  P 1  in 
-  p  |v0  m  m 


—  P  |P0 

cc  *ojtv 


©  ler  i 

®>®  i 


o  fr  I© 
0  5*10 
IP  CM1*- 

*o  *oip 


p|9 

P>4 


- 1  m  - 


♦  p  * 

in  ©  |©  < 


©  fr 
to  © 


1  P  "M  • 

PI©  < 


P  PI© 
■r  P|o 

i©  **1© 


i  4*  /*-.  i 

I  ©IP  I 

‘  *i  * 

eelp  i 

pfp  - 
in-*o  i 


rv  *oip  ir  Up  p 


HioaoTcjeHtrsoi _ . 

DrtNMfoifjNininN'Dn) 


©OPOOCOPOPOP 
oopovoapp^pop) 
v  P'  ©  45  p  *0  nn  opn 
cn  h  n  m  (r  c  a  a  co)  in  ml 


•o  p  ,  p  p| 
p  ©  ©  mi 


:  fr'  © 
I  *4'  CO 
i  o:  o 


p  p  i 

cc  i 


•  o  p  o.  o  P 
i  O  *  m  N  m 
I  •**  CO  *0  ©  «- 


1  '  J 

(o  c  m  p-  n  c 

m  in  vO  h 


p  p 
in  4} 


in  co  mi 

m  VD  CO 


co  m 

©  p 

i  m  m 


irt  ©  i 

mi  m  < 


)  P!  © 
»  ♦  • 
I  -4  © 
■  <r  ® 
I  O'  p 


P  K>  . 
TO  P  « 

p,  m  c 


.1 


a>  roi 
p  m<  cm 
•  *i  • 
p  pip 
P  P!  © 
cm  cotm 


HHinH© 

©  4)  ©Jp 


©  p<  -4  ir?' ' 


4*1  * 

PH 

m  45'T 


P'p 
fr  p 
vO-  vO 


iC  h  «r  ol 


P(P 
colrn 
© !  43 
CM  |  *4 


-•■m  i 

Pi<V  < 


»  ip'm 

I  o|* o 


Ip  p 

p  o 


©  vfl 
!»•: 
ip  o 

p  © 
|  CM  P 
W  vO 


pi©  i 

•  i  • 


«  Pi© 

■  w  *n 
»  p/p 
i  4»im 


•  ©|p  i 
»  ©  l  cc  * 

I  vfi.CE  l 

‘  P|^  f 


P  P 
•  • 
p  p 
P  1 


I  4J|-i  v 

i  • 1  • 

■  ©I©  4 
•  mi©  < 

i  f-ip  • 
»  CM  1*0  * 


P  O 

vo  m 

K>  *0 


p  m 
p  © 


p  P.m 


cm  rvl®  m 
*0  010* 


/  m  o  i 
Mini 


•  * 


i  4>  P 
©>•-4 


.  K)  m  CM 
•  p*»o  m| 
I  Pi  4)  mi 

i 

I 

O  CD  O I 
p  *o  m 


I  p.vfl 

>  4)  0 
©  © 


P  '  CM  ©| 
o  ~4  m 
p  *o  o1 


»  p'in  i 
mm  i 


I  T-4-© 

I  ©'*. 

1  in  p 


•  P!© 

i  m;in 


•m!  P  p 
ml  CM  -4 

PIP  pi 


in 'cm  col 

p|©  p, 


pi© 

~«i«* 

t/*im 


p  m 
inp 


PIP  o 
hIP  *0 
WIP  p, 


PI©  CO 
in  1*0 


p  jp  © 
pi©  p 

mip 

m  in | in  in|©  © 


•  r  ip  i 
mi©  i 


»|«*  p 


»  ©IP  4 

I  ©M  • 

01*0  4 
'*■!©» 


'  fr  fr  . 
*0(fV  t 


p  © 

m  © 
I  cm  *n 
®  © 


P  w4 

®  rv 


r  ©  m  p 
r  p|m  m 
•  •!  •  * 


fr  © 
tm  *o 

p  © 

p  © 


p  © 
co  P 


■  i 

1 1 1  - 1 3 


©  K> 
O  P 


CM  O 
©  CJ 


P  P 
P  K> 

*o  *o 


<M  « 
P  © 

P  P 


■O  © 
*0  P 


*o  © 

*-  CM 
©  © 


©  © 

•  • 

O  P 


P  o 

1*0  *o 
jp  o 
©  © 


c  inh 
p  rv  |o 


m  p 

p 

©  © 

r- 

p 

CM  CO 

O'H  H  CO 

© 

© 

CO 

CM 

© 

© 

o 

©ip 

c— 

PlP  CO’P 

© 

O 

O  P 

p 

o  p 

o 

Pi© 

p 

P  O  CO i CM  COI 

O 

OP 

p 

in 

•r 

m 

P  05  P 

P'O  ©<© 

CM 

fr 

© 

O 

m 

p 

co 

CJ 

Hi  O 

© 

© 'p  p  © 

CO 

CO 

*oco 

CM 

CM 

cv 

CO 

CM 

•H 

«-»  ®  ml© 

p 

•  •*4  'P 

■  •-4 

© 

p 

® 

CM 

P 

fr  1© 

p 

o 

o 

fr  'p  ~4-*+ 

o 

CM 

♦-41© 

© 

®  o 

o 

p 

cc 

o 

P  O  43  P 

m 

P 

CM‘P 

o 

© 

©■© 

0? 

* 

o 

POOi© 

© 

CM 

^0 

•-4 

p 

in 

ro 

©  *o  m 

h  m  in  •-* 

co 

*o 

m 

•-4 

COICM 

© 

H 

o  <co  ©  m 

CO 

*-4 

CM 

P 

p 

m  m 

p 

CO 

CO 

— *  tn  cm  -4 

•H 

‘© 

O 

© 

o 

© 

p 

PIO 

p 

p  p  p  >® 

© 

P 

CO 

P 

© 

© 

© 

T- 

c4 

« 

P  ©  CC  ;P 

© 

P 

P 

© 

p 

♦o 

O 

m  co 

CM 

CO 

P  *0  co-p 

m 

P 

CO 

p 

m  picm 

PI®  If) 

ro  m  in  © 

CM 

CO 

CO 

CO 

co 

CO  T-> 

p 

© 

fr 

o 

CM  O  m-«H 

p 

•H  CO 

•CM  -P 

•w* 

CM 

r4 

CM 

in 

O 

P 

p 

1 

P  1-4 

o 

m  -p  tO‘*o 

C-. 

P 

*- 

- 

© 

© 

CO 

o 

CC  lO 

o 

pc  m  ©  © 

o 

©im 

CO 

c»» 

©  1©  p 

H 

o 

in  i*o  cv  © 

T- 

in 

w 

© 

CO 

© 

p 

CO 

o 

CM  'P  P  -P 

o 

p 

p 

P 

*n 

O 

CO 

CO  CO 

*o 

-*  cc  tn  p 

p 

1^4  'P 

-•4 

o 

©  p 

© 

P 

p 

P 

p 

© 

p 

P  ,©  •-*  'P 

*o 

© 

© 

•4 

p 

P 

CO 

P 

©  © 

©  (©  ©  < 

rj 

© 

H 

O 

p 

TO 

© 

CV  o 

© 

o 

o  p  P  ‘CV  P  | 

fr 

pirn 

CO 

•— ( 

p  .©  p 

m  m 

in  ©  ro  © 

p 

CO 

co 

CM 

CM 

CM 

«-4 

•4 

•H 

—t 

—  ©  *0  »*-4 

i© 

•«H 

CT 

CO 

O 

© 

•o 

p  i© 

m 

m  ©  ©  © 

© 

© 

O/CV  <v 

cv  in 

*7 

*?*? 

co  ICO  m  CM 

p 

to  *o»p 

4* 

m 

r-4 

p 

rv 

kF 

cr 

©  4*4  ©  '*0 

«-4 

m 

p 

er 

in 

© 

n 

o 

Cs. 

© 

co 

*o  *p  ©  © 

K. 

p 

<r 

p 

p 

p 

p 

p 

•O  K5 

CO 

*o  to  in  m 

—4 

— *  <T- 

Ts. 

p 

CM 

p 

o 

© 

«■ 

CO 

O 

o 

O  -H  O  W 

CM 

*• 

^|m  in 

OJ 

O  JCM 

* 

o 

CM  P  CJ  '.CC 

in 

P 

©ICO 

«• 

p 

Tw 

o 

© 

E 

o 

»  »  ob 

cc 

©  W 

P 

© 

CO 

CM  1© 

m 

(MAI  CIO 

© 

P 

*- 

»** 

© 

©  p 

*- 

m 

© 

©  P  CO 

m 

—4 

-H 

•rv  ~4 

i  p 

« 

CJl© 

CM 

K. 

CM 

CM 

«--L 

o 

Plea 

fr 

K. 

»o 

*? 

* 

in 

•o  ro  CO  p 

© 

J. 

©IP 

o 

© 

P 

p  p 

O 

p 

DO 

© 

in 

mi© 

© 

© 

P 

©  © 

© 

© 

©  *o  m  /cv  © 

p 

•-* 

© 

© 

ro  co 

<r 

P 

JJ 

p 

©  ro  in  — 

*-< 

•o 

♦ 

P 

CO 

■o 

*> 

CJ 

V 

CO 

CV  P  ©  TO 

*o 

r-  *° 

j  P/ 

p 

© 

O 

p 

CO 

1*0 

p 

© 

«h  ter  p  fp 

o 

"* 

o 

Ip 

m 

m 

CV 

•O  © 

fr 

o 

to  P  ©  TO 

•H 

o 

p 

fr 

fr 

p 

p 

O 

p 

o 

fr  p  p  Ip 

CM 

—4 

p 

*> 

cO 

© 

p 

P 

© 

s 

co  tr  fr  >cr 

T- 

m 

p 

m 

*o 

p 

n 

<=*• 

KC  P  fcv 

P 

» 

*■ 

p 

lO 

CM 

*0 

© 

© 

ro  sk* 

CO 

in  ^  o 

K 


I©  ote 
K  r  ro 
r  ol® 
P  ©kO 


sfc: 


OM  Mh  CO 
P  io  *-U  p 

•  i  •  •  I  •  • 

©  ho  ©|p  p 
©  tm  o  |©  m 
inh)  «h  ee 
fr  p*  M  P  P 
r  *"4 1 

©  ho  |©  *o 

©ICO  ©p>  »0 

•  I  •  •  «  • 
obobo 
CM  »o  CM  Pi  « 
MH  mb  * 

in  im  p  *o  © 


©  w4  ftv  #n  ip  m 
CM  CM  <M  (M  jC*  «M 


I©  p  fcr  P  ©  « 

prso  cm  rv  rv  p  • 


1  c 

«r  to  p  f-* 
p  p  ©U 

p  jp  t-  rv 


p  o  r“  o 


Np'ceom 

f  £ 

©  45  ©  H  © 

P  O  (M  p  O 

•  '  t  •  I  •  • 

enmn  co 
p  m  ©  P  ** 
ir  n  p  h  n 
cm  p  m  p  *o 

r  ts 

i 

P  *o  I*  p  ro 
Mn*>.«* 

•  •  •  I  •  • 
P  O  p  rv  © 

KttQOO 
o  P  -O  I©  © 
TCBOlP* 

F 

7  ftc  2  LI  > 

T-  «C  P  IP  U 

r  t  u  p  o 
r  o 
p  © 


i 


FIGURE  JII-7.  Example  of  Model  Output 


approach  is  deemed  most  appropriate  since  it  utilized  the  best  available 
data  and  provides  the  desired  results. 

The  amount  of  incremental  inflow  between  two  adjacent  stations 
as  calculated  in  the  INFLOW  program  is  simply  the  mean  daily  flow  at  the 
downstream  station  minus  the  mean  daily  flow  at  the  upstream  station. 

The  upstream  flow  is  lagged  prior  to  subtraction;  the  amount  of  lag  de¬ 
pends  on  the  time  of  travel  between  stations.  If  the  wrong  amount  of 
lag  is  applied  to  the  upstream  flow,  negative  inflows  result  between 
stations.  But  the  same  lag  assumption  is  used  to  disaggregate  flows  in 
the  INFLOW  program  as  is  used  to  combine  the  inflows  in  the  NATFLO  pro¬ 
gram,  so  the  negative  inflows  have  no  effect  on  the  natural  flow  simula¬ 
tion.  The  problems  with  negative  inflows  occur  at  nodes  just  downstream 
of  reservoir  sites  when  regulation  is  added  to  the  system.  In  most 
cases,  the  conservation  releases  from  the  reservoirs  are  much  smaller 
than  the  naturally  occurring  discharges.  Thus,  a  small  upstream  dis¬ 
charge  plus  a  negative  incremental  inflow  could  yield  negative  or  ab¬ 
normally  low  discharges  at  nodes  just  downstream  of  reservoir  sites.  If 
abnormally  low  or  negative  flows  were  allowed  to  occur,  the  statistical 
portions  of  this  project  would  become  meaningless. 

It  is  well  known  that  time  of  travel  is  a  function  of  flow.  The 
actual  function,  however,  is  not  well  defined.  Several  attempts  are  made 
to  determine  the  best  lag  time  function  between  each  pair  of  nodes.  As  a 
baseline  case,  the  first  attempt  at  generating  inflow  assumes  zero  lag 
between  all  stations.  This  proves  to  be  wholly  inadequate.  There  are 
often  many  days  in  a  row  which  display  negative  inflows.  The  next 
attempt  uses  a  constant  lag  between  the  various  stations.  While  several 
different  constant  lags  are  tried,  all  prove  to  be  equally  unsatisfactory. 
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The  third  attempt  begins  with  a  derivation  based  on  Manning's 
Equation,  and  results  in  an  equation  of  the  form: 

L  =  aQB  (III-l) 

where 

L  =  flow  through  time  of  a  given  reach, 

Q  =  discharge  in  that  reach  at  the  given  time;  and 
a,£  =  empirically  derived  constants. 

Unfortunately  there  is  insufficient  travel  time  information  to 
properly  determine  the  coefficients  "a"  and  "8."  For  this  reason,  this 
method  also  proves  to  be  unsatisfactory. 

The  final  and  selected  method  for  lag  time  determination  assumes  a 
linear  lag  function  of  the  form, 

L  =  aQ  +  b  ( 1 1 1 -2 ) 

where 

L  =  flow  through  time  of  a  given  reach, 

Q  =  discharge  in  that  reach  at  the  given  time,  and 
a,b  =  empirically  derived  constants. 

The  constants,  a  and  b,  are  initially  determined  from  available  times  of 
travel  and  high  and  low  flow  data.  High  and  low  flow  data  are  readily 
available  for  each  USGS  station.  The  maximum  and  minimum  discharges  which 
are  published  in  the  surface  water  records  are  used  in  determining  the 
constants.  The  initial  values  of  low  and  high  flow  times  of  travel  are 
developed  from  the  "Report  on  the  Comprehensive  Survey  of  the  Water  Resources 
of  the  Delaware  River  Basin. Reasonable  travel  times  are  selected  for  each 
reach.  High  and  low  flow  times  of  travel  are  estimated  using  the  stream  miles 
between  stations  and  average  rates  of  travel  when  not  available  in  the 
report.  The  constant,  b,  is  the  low  flow  time  of  travel  and  represents  the 
maximum  time  between  stations.  The  maximum  travel  time  between  stations  can 


*  "Report  on  the  Comprehensive  Survey  of  the  Water  Resources  of  the  Delaware 
River  Basin,  Appendix  M:  Hydrology"  U.S.  Corps  of  Engineers,  Philadelphia 
District,  1960. 


III-15 


never  be  greater  than  24  hours,  since  the  model  operates  with  a  daily  time 
step.  The  constant,  a,  is  determined  by: 


a  _  Low  Flow  Travel  Time  -  High  Flow  Travel  Time  (III-3) 

Low  Flow  -  High  Flow 

and  is  always  a  negative  value. 

The  times  of  travel  are  further  adjusted  from  the  initial  values 
selected.  The  times  of  travel  are  adjusted  so  that  the  fewest  number  of 
negative  inflows  are  created  when  the  recorded  data  are  lagged  and  subtracted. 
Negative  inflows  are  still  occasionally  produced,  but  these  negatives  pose 
no  problem  to  the  model  except  at  those  nodes  just  downstream  of  an  existing 
or  potential  reservoir  site.  Table  1 1 1-2  shows  the  final  range  of  lag  times 
used  in  the  modelling  effort.  The  constants  "a"  and  "b"  for  the  lag  function 
are  also  shown  in  Table  III-2. 

The  lag  function  described  above  is  applied  in  the  model  simulation 
in  the  following  manner. 

For  each  day  the  total  flow  at,  say,  node  I  is  the  sum  of  the 
incremental  inflow  to  the  node  and  the  routed  flow  from  the  upstream  node, 
i-l.  The  total  routed  flow  from  the  upstream  node  is  calculated  based  on 
routing  a  fraction  of  the  day-before  flow  at  the  upstream  node  and  routing 
a  fraction  of  the  present-day  flow  at  the  upstream  node. 

The  flow  at  node  I  for  any  given  day  is  expressed  as: 

FLOW(I)  =  INFLOW ( I )  +  QOUT(I-l)*  YLA6  +  Q(I  -  1)*(1-XLAG)  (III -4) 

where 

FLOW ( I ) 

INFLOW(I) 

QOUT(I-l) 

YLAG 


=  Present-day  total  flow  at  node  I 

=  Present-day  incremental  inflow  to  node  I 

=  Day-before  flow  at  upstream  node 

=  lag  function  applied  to  day-before  total  flow 
at  upstream  node 
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TABLE  1 1 1 -2 

RANGE  OF  HIGH  AND  LOW  FLOW  LAG  TIMES 
AND  LAG  FUNCTION  CONSTANTS 


Reach 

Lag  Time 
(Hours) 

Low  Flow  High  Flow 

Lag  Function 
a 

(in  hr/cfs  10-3) 

Constants 

b 

(hours) 

Downsville  -  Fishs  Eddy 

16 

0 

-0.6695 

16 

Fishs  Eddy  -  Barryville 

24 

6 

-0.3377 

21 

Stilesville  -  Hale  Eddy 

6 

0 

-0.3429 

6 

Hale  Eddy  -  Barryville 

24 

6 

-0.6429 

24 

Barryville  -  Port  Jervis 

10 

3.5 

-0.0500 

10 

Prompton  -  Lackawaxen  at  Honesdale 

2 

0 

-0.3413 

2 

Dyberry  at  Honesdale  -  Lackawaxen  at  Honesdale  1 

0 

-0.0645 

1 

Lackawaxen  at  Honesdale  -  Hawley 

5 

2 

-0.1613 

5 

Hawley  -  Port  Jervis 

11 

5 

-0.1156 

11 

Wilsonville  -  Port  Jervis 

4 

2 

-0.3106 

4 

Neversink  -  Oakland  Valley 

18 

0 

-0.8072 

18 

Oakland  Valley  -  Montague 

10 

2 

-0.2667 

10 

Montague  -  Belvidere 

24 

0 

-0.0960 

24 

Belvidere  -  Riegelsville 

4 

0 

-0.0062 

4 

White  Haven  -  Walnutport 

16 

5 

-0.2030 

16 

Beltzville  Dam  -  Walnutport 

9 

3 

-1.1320 

9 

Palmerton  -  Walnutport 

3 

1 

-0.1786 

3 

Walnutport  -  Bethlehem 

16 

4 

-0.1542 

16 

Allentown  -  Bethlehem 

4 

1 

-0.1852 

4 

Bethlehem  -  Riegelsville 

12 

3 

-0.0978 

12 

Hackettstown  -  Bloomsbury 

9 

4 

-2.3040 

9 

Bloomsbury  -  Riegelsville 

3 

2 

-0.1437 

3 

Riegelsville  -  Trenton 

14 

0 

-0.0176 

14 

Pipersville  -  Trenton 

8 

5 

-0.1875 

8 

Trenton  -  Torresdale  . 

9 

5 

-0.0122 

9 

Assunpink  Cr.  at  Trenton  -  Torresdale 

9 

5 

-0.7339 

9 

Extonville  -  Torresdale 

10 

6 

-0.7722 

10 

Langhorne  -  Torresdale 

4 

2 

-0.0406 

4 

Pemberton  -  Torresdale 

8 

5 

-0.1734 

8 

Pottsville  -  Landingville 

3 

1.5 

-0.2431 

3 

Cressona  -  Landingville 

2 

0 

-1.8180 

2 

Langingville  -  Berne 

5 

2.5 

-0.2917 

5 

Tamaqua  -  Berne 

7 

4 

-0.3851 

7 

Berne  -  Schuylkill  at  Reading 

9 

5 

-0.0935 

4 

Schuylkill  at  Reading  -  Pottstown 

18 

11 

-0.1827 

18 

Pottstown  -  Philadelphia 

24 

11 

-0.2559 

24 

Graterford  -  Philadelphia 

11 

5 

-0.1504 

11 

Philadelphia  -  Delaware  R. 

5 

2 

-0.0312 

5 

Chester  -  Delaware  R.  at  Chester 

2 

1 

-0.0694 

2 

Newark  -  Delawar  Memorial  Bridge 

8 

4 

-0.2326 

8 

Wooddale  -  Delaware  Memorial  Bridge 

8 

3.5 

-0.2528 

8 

Chadds  Ford  -  Wilmington  (Brandywine  Cr.) 

5 

2.5 

-0.4167 

C 

Wilmington  (Brandywine  Cr)  -  Del.  Mem.  Br. 

3 

1.5 

-0.2366 

3 
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Q ( I  -  1)  =  Present-day  total  flow  of  upstream  node 

(1  -  XLAG)  =  Lag  function  applied  to  present-day  total  flow 
at  upstream  node 

In  the  model,  the  lag  function.  Equation  III-2,  is  transformed  into  two 

expressions.  For  the  day-before  flow  at  the  upstream  node  which  is  routed 

and  accounted  for  in  the  present-day  flow  at  the  downstream  node  the  lag 
function  is: 

YLAG  =  a  *  QQUT( I  -  1)  +  b  (III-5) 

where 

YLAG  =  Fraction  of  the  day-before  flow  at  the  upstream 

node  routed  to  the  downstream  node 

QOUT(I-l)  =  Day-before  flow  at  upstream  node 

a,b  =  Lag  Constants 

For  the  present-day  flow  at  the  upstream  node  which  is  routed  and  accounted 
for  in  the  present-day  flow  at  the  downstream  node  the  lag  function  is: 

XLAG  =  a  *  0(1-1)  +  b  ( 1 1 1-6) 

where 

XLAG  =  Fraction  of  present-day  flow  at  the  upstream  node 

which  is  not  routed  to  downstream  node.  Fraction  of 
present-day  flow  which  routed  to  downstream  node 
is  calculated  in  Equation  1 1 1 - 4  as  (1  -  XLAG). 

Q(I-l)  =  Present-day  flow  at  upstream  node 

a,b  =  Lag  Constants  (identical  to  a,  b  in  III-5) 

Consider  the  following  example.  An  upstream  node  (1-1)  has  a  flow 
of  1500  cfs  for  the  present-day  flow  and  1000  cfs  for  the  da.y-before  flow. 
The  flows  are  to  be  routed  to  the  downstream  node  (I)  with  no  incremental 
runoff  to  the  downstream  node  (INFLOW  (I)  =  0).  The  "a''  and  "b"  of  the 
lag  function  for  the  reach  are  -0.0003  and  24,  respectively.  In  the  model 
the  constants  "a"  (hr/cfs)  and  "b"  (hrs)  are  divided  by  24  to  produce  a 
fraction  of  a  daily  lag  time. 
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YLAG,  the  percent  of  flow  routed,  is  applied  to  the  day-before 
flow  at  the  upstream  node  as  shown  in  Equation  1 1 1 -4  and  is  calculated 
using  Equation  1 1 1 - 5  as: 

YLAG  =  (-0.0003/24)  (1000)  +  24/24  =  0.9875 

where 

1000  =  the  day-before  flow  at  upstream  node 

Ninety-nine  percent  of  the  day-before  flow  at  the  upstream  node  is  routed 
to  the  downstream  node. 

XLA3  is  calculated  using  Equation  1 1 1 -6  as: 

XLAG  =  (-0.0003/24) ( 1500)  +  24/24  =  0.9813 

where 


1500  =  the  present-day  flow  at  upstream  node. 

Thus,  the  amount  of  the  present-day  upstream  node  flow  routed  to  the  down¬ 
stream  node  is  1.87  percent  or  (1  -  XLAG)  as  shown  in  Equation  III-4. 

The  final  calculation  of  the  downstream  node  flow  for  the  present- 
day,  as  computed  by  equation  III -4 ,  with  no  increment  1  runoff  and  a  fraction 
of  the  day-before  flow  and  a  fraction  of  the  present-da^  flow  from  the 
upstream  node  is: 

FL0W(I)  =  0  +  1000(0.9875)  +  1500  (1  -  0.9813)  =  1015.55 

Negative  inflows  are  sometimes  produced  independent  of  the  lag 
functions  applied  during  the  subtraction  of  recorded  flows.  These  negative 
inflows  could  be  the  result  of:  inaccuracies  in  the  recorded  data  (most 
gages  record  to  +  15  percent);  small  unaccountable  diversions;  or  infiltra¬ 
tion  and  evaporation  between  USGS  stations  (particularly  on  parts  of  the 
Schuylkill  and  Lehigh  Rivers  where  the  bedrock  is  limestone). 
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The  models  can  handle  negative  inflows  unless  they  produce  ab¬ 
normally  low  or  negative  flows  when  combined  under  regulated  conditions. 
The  only  nodes  which  cannot  readily  accept  negative  inflows  are  those 
just  downstream  of  an  existing  or  proposed  reservoir  site.  If  such  a 
node  has  a  negative  inflow,  the  negative  inflow  could  produce  a  negative 
or  abnormally  low  flow  when  only  a  conservation  release  is  added  from 
upstream  under  regulated  conditions.  The  inflows  to  the  nodes  downstream 
of  the  existing  and  proposed  reservoir  sites  are  examined  carefully  for 
negative  inflows.  The  nodes  at  which  there  are  problems  include: 

Fishs  Eddy, 

Hale  Eddy, 

Lackawaxen  at  Honesdale 
Oakland  Valley, 

Allentown,  and 
Tulpehocken  at  Reading. 

The  times  during  which  negative  inflows  are  a  recurring  problem 
for  these  gages  are  listed  on  Table  1 1 1 -3 .  During  these  times,  the  in¬ 
flows  are  created  in  a  different  manner.  The  new  methods  of  creating 
inflows  for  each  troublesome  node  are  also  listed  on  Table  1 1 1 - 3 .  For 
anodewhich  has  negative  inflows  during  a  regulated  period  (such  as  Hale 
Eddy),  the  new  inflows  to  the  node  are  calculated  as  a  percentage  of  the 
unregulated  inflow  at  an  upstream  station.  The  percentage  is  estimated 
using  a  drainage  area  ratio.  For  a  node  which  has  negative  inflows 
during  an  unregulated  period  (such  as  Tulpehocken  at  Reading),  new  in¬ 
flows  are  calculated  for  both  the  upstream  and  downstream  nodes  based 
on  the  recorded  flow  at  the  downstream  gage. 

After  the  inflows  to  the  nodes  downstream  of  reservoir  sites  are 
corrected  for  problems  with  recurring  negatives,  the  inflows  to  those 
nodes  downstream  of  reservoir  sites  are  processed  to  remove  any  stray 
negatives.  The  remaining  negative  flows  are  made  equal  to  zero  and 
carried  until  the  following  day’s  positive  inflows  displace  them.  No 
more  than  three  consecutive  days  of  zeros  occur  with  this  process. 
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TABLE  1 1 1-3 

FIXES  FOR  RECURRING  NEGATIVE  INFLOWS 
T  NODES  DOWNSTREAM  OF  RESERVOIR  NODES 
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This  process  is  not  applied  to  the  inflows  to  the  nodes  on  the 
Schuylkill  at  Reading  and  Lehigh  at  Walnutport.  Negative  inflov/s  are 
allowed  at  these  stations,  because  they  are  main  stem  stations  and  there 
is  little  danger  of  producing  unusually  low  flows  with  the  proposed 
relatively  sliqht  regulation. 

After  the  final  adjustments  to  the  lag  functions,  the  negative 
inflows  which  occur  at  the  other  nodes  in  the  model  are  left  intact. 

These  negatives  result  from  lagging  and  subtracting  recorded  data.  By 
leaving  these  negatives  in  the  input  data,  the  recorded  data  are  best 
preserved  and  continuity  is  maintained  for  flow. 

RESULTS  OF  THE  NATURAL  FLOW  MODEL 

As  in  any  modelling  effort,  the  model  which  is  to  be  used 
as  a  predictive  tool  has  the  ability  to  reconstitute  recorded  data.  With 
this  particular  project,  the  shear  volume  of  this  data  makes  this  type  of 
comparison  an  arduous  task.  Since  the  USGS  records  prior  to  1953  are  not 
affected  by  regulation,  the  natural  flow  model  simulation  is  compared  to 
these  records.  Because  of  the  amount  of  output  generated  by  the  model,  this 
comparison  is  done  in  a  two  step  process.  First,  to  assure  that  continuity 
is  preserved,  the  monthly  average  discharge  for  each  of  the  simulated 
stations  is  compared  to  the  data  presented  in  the  various  USGS  Surface 
toaier  Supply  Records  publications.  As  an  example,  recorded  flows  and 
simulated  flows  are  presented  in  Figures  1 1 1 -8  through  1 1 1 - 10 .  Figure 
1II-8  gives  the  monthly  mean  discharge  for  water  years  1951  through  1960 
for  the  East  Branch  Delaware  River  at  Fishs  Eddy  as  presented  in  the  U.S. 
Geological  Survey  Water  Supply  Paper  1722.  Figures  1 1 1 -9  and  1 1 1 - 10  present 
the  model  flow^  at  Fishs  Eddy  for  1953  and  1954,  respectively.  The  monthly 
mean  flows  are  shown  at  the  bottom  of  the  daily  flow  table.  A  comparison 
of  the  monthly  mean  flows  for  1953  and  1954  show  that  the  simulated  flows 
match  the  recorded  flows  when  rounded  to  three  significant  figures.  Note 
that  in  months  subsequent  to  September  1954  the  USGS  records  of  Fishs  Eddy 
do  not  agree  with  the  simulated  records.  This,  of  course,  is  because  the 
simulated  records  have  eliminated  the  effect  of  the  Pepacton  Reservoir 
which  became  operational  at  that  time. 
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DELAWARE  RIVER  BAS XV 


4210.  East  Branch  Delaware  River  at  Flshs  E ddy,  R.  Y. 

Location  .—Lat  41*58' 00*.  Iona  75*10'5Q",  on  left  bank  at  downstream  side  of  highway 
[revised)  bridge  at  Flaha  Eddy,  Delaware  County,  Just  upatreaa  from  Fish  Creek, 

4*  allei  downstream  fro*  Beaver  Kill  and  11  alias  upstreaa  from  confluence  of  East 
and  Vest  Branches  near  Hancock. 

Drainage  area.— 783  aq  ml. 

Records  available.-- October  1912  to  September  1960. 
periods,  published  in  WSP  1302. 


Monthly  discharge  only  for  some 


Qa&f ar~3 ta*e  mortar.  Datv«  of  gage  Is  950.96  ft  above  mean  sea  level,  datum  of 
*929.  Prior  to  Sept.  27,  1928,  staff  gage  at  same  site  and  datum. 

Extremes.— 1912-60:  Maximum  discharge,  53,300  cfs  Aug.  24,  1933  (gage  height,  20.60  ft), 
from  rating  curve  extended  above  22,000  cfs  by  logarithmic  plotting;  mint  rum,  67  cfs 
Aug.  28,  1949;  minimum  gage  height,  1.51  ft  Aug.  4.  5,  1935. 

Plood  of  Oct.  9,  1903,  reached  a  stage  of  23.6  ft,  from  descriptions  obtained  In 
April  1939  from  local  residents  who  had  experienced  the  flood  (discharge,  about  70,000 
cfs,  from  rating  curve  extended  above  22,000  cfs  by  logarithmic  plotting). 

Remarks . —Subsequent  to  September  1954,  entire  flow  from  371  5q  ml  of  dnlnare  area  con¬ 
trolled  by  Pepacton  Reservoir.  Part  of  flow  diverted  for  municipal  supply  of  city  of 
New  York.  Remainder  of  flow  (except  for  conservation  releases  sr.d  spill)  impounded 
for  release  during  periods  of  low  flow  lr.  the  lower  Delaware  River  basin,  as  directed 
by  the  Delaware  River  Master.  Records  of  chemical  analyses  for  the  period  October 
1957  to  September  19S9  are  published  In  reports  of  the  Geological  Survey. 
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FIGURE  1 1 1-8.  Monthly  Mean  Discharges  at  East  Branch 
Delaware  River  at  Fishs  Eddy  as  Reported  by  the  U.S.G.S 
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FIGURE  III-10.  Simulated .Natural  Discharge  1954:  East 
Branch  Delaware  River  at  Fishs  Eddy 


After  demonstrating  that  continuity  is  preserved  on  a  monthly  basis 
for  each  of  the  sixty  stations,  the  simulated  daily  flows  are  compared  to 
assure  their  validity.  Figures  1 1 1 - 1 1  through  1 1 1-17  show  typical  compari¬ 
sons  of  the  simulated  daily  data  to  data  recorded  by  the  US6S.  Figure  III-ll 
presents  recorded  flows  of  various  stations  for  the  first  eight  days  of 
January  1953.  Six  gaging  stations  are  underlined  in  the  figure  for  which 
simulated  flows  are  presented  in  Figures  1 1 1 - 12  through  III-17,  respectively. 
For  example,  the  flows  of  1-8  January  1953  for  the  first  station  underlined 
(01417000)  in  Figure  III-ll  compare  favorably  to  the  first  eight  days  of 
January  for  this  station,  E.  Branch  Delaware  River  at  Downsville,  in 
Figure  1 1 1 - 12 .  The  simulated  daily  data  satisfactori ly  match  the  unregu¬ 
lated  daily  records  of  the  USGS.  Variations  between  the  simulated  and 
recorded  data  is  attributed  to  round-off  applied  during  the  procedures 
used  to  develop  the  total  flow  at  the  gage  locations. 

Duration  and  Frequency  Analysis 

The  natural  flow  model  for  the  Delaware  River  Basin  simulates  the 
natural  daily  flows  for  a  50-year  base  period  from  water  year  1928  through 
water  year  1977,  a  flow  duration  and  low  flow  frequency  analysis  is  per¬ 
formed  on  the  50  years  of  naturalized  flows  for  the  following  44  key 
locations: 

1.  East  Branch  Delaware  River  at  Downsville,  N.Y., 

2.  East  Branch  Delaware  River  at  Fishs  Eddy,  N.Y., 

3.  West  Branch  Delaware  River  at  Stilesville,  N.Y., 

4.  West  Branch  Delaware  River  at  Hale  Eddy,  N.Y., 

5.  Delaware  River  near  Callicoon,  N.Y., 

6.  Delaware  River  near  Barryville,  N.Y., 

7.  West  Branch  Lackawaxen  River  at  Prompton,  PA, 

8.  Dyberry  Creek  near  Honesdale,  PA, 

9.  Lackawaxen  River  at  Honesdale,  PA, 

10.  Lackawaxen  River  at  Hawley,  PA, 

11.  Delaware  River  at  Port  Jervis,  N.Y., 

12.  Neversink  River  at  Neversink,  N.Y., 
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Flows  for  1-8  January  1953 
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13.  Neversink  River  at  Oakland  Valley,  N.Y., 

14.  Delaware  River  at  Montague,  N.J., 

15.  Delaware  River  below  Tocks  Island  Damsite,  PA., 

16.  Delaware  River  at  Belvidere,  N.Y., 

17.  F.E.  Walter  Reservoir  near  White  Haven,  PA., 

18.  Pohopoco  Creek  below  Beltzville  Dam,  PA., 

19.  Aquashicola  Creek  at  Aquashicola  Damsite,  PA., 

20.  Aquashicola  Creek  at  Palmerton,  PA., 

21.  Lehigh  River  at  Walnutport,  PA., 

22.  Jordan  Creek  near  Schnecksville,  PA., 

23.  Jordan  Creek  at  Allentown,  PA., 

24.  Lehigh  River  at  Bethlehem,  PA., 

25.  Lehigh  River  at  Glendon,  PA., 

26.  Musconetcong  River  near  Hackettstown,  N.J., 

27.  Delaware  River  at  Riegelsville,  N.J., 

28.  Tohickon  Creek  near  Pipersville,  PA., 

29.  Delaware  River  at  Trenton,  N.J., 

31.  Schuylkill  River  at  Cressona,  PA., 

32.  Schuylkill  River  at  Landingville,  PA., 

33.  Little  Schuylkill  River  at  Tamaqua,  PA., 

34.  Little  Schuylkill  River  at.  Diehersvil le,  PA., 

35.  Schuylkill  River  at.  Berne,  PA., 

36.  Maiden  Creek  at  Viroinvil le,  PA., 

37.  Tulpehocken  Creek  at  Pine  Marsh  Pansite,  PA., 

38.  Tulpehocken  Creek  near  heading,  PA., 

39.  Schuylkill  River  near  Reading,  PA., 

4D.  Schuylkill  Piver  at  Pottstown,  PA., 

41.  Perkionen  Creek  at  Craterford,  PA., 

42.  Schuylkill  River  at  Philadelphia,  pA. , 

43.  Delaware  River  below  Schuylkill  Confluence,  and 

44.  Delaware  River  at  Delaware  Memorial  Bridgp. 

The  duration  and  low  flow  frequency  tables  and  plots  from  the  natural 
flow  model  are  presented  in  Appendix  A, 


The  flow  duration  analysis  consists  of  the  development  of  a  cumu¬ 
lative  frequency  distribution  which  shows  the  percentage  of  time  the 
indicated  flows  have  been  equaled  or  exceeded.  The  distribution  does 
not  take  into  account  the  chronological  sequence  of  flows  and  therefore 
does  not  indicate  whether  varying  periods  of  low  flow,  for  example, 
occurred  during  one  dry  weather  period  or  were  distributed  over  several 
years. 


The  low  flow  frequency  analysis  shows  the  relationship  between  the 
magnitude  and  frequency  of  the  annual  lowest  mean  flow  for  a  given  number 
of  consecutive  days.  For  each  n-consecutive  days  analyzed,  the  log- 
Pearson  Type  III  method  is  used  to  produce  theoretical  values  which 
correspond  to  non-exceedance  probabil ities  or  recurrence  intervals  in 
years. 


The  duration  and  frequency  results  presented  in  Appendix  A  are  for 
50  years  of  naturalized  flows  at  the  44  key  locations  throughout  the 
Delaware  River  Basin.  These  characteristics  are  the  basis  for  comparison 
to  the  regulated  flow  characteristics  of  the  basin  with  the  operation  of 
the  three  New  York  City  reservoirs  during  the  50-year  base  period. 

The  U.S.  Geological  Survey's  Daily  Values  Statistics  Program  (A969) 
is  used  to  develop  the  duration  and  flow  frequency  tables.  The  Philadelphia 
Corps  of  Engineers  executed  the  program  on  the  Boeing  computer  system  in 
Vienna,  Virginia  from  the  input  data  supplied  by  WRE  to  them  on  magnetic 
tape  and  cards. 

The  A969  program  is  used  to  compute  the  number  of  times  flow  at  a 
node  equaled  or  exceeded  a  certain  magnitude.  The  program  tabulated 
the  number  of  daily  values  in  each  of  35  magnitude  classes  for  each  year 
for  each  key  location.  The  class  limits  are  computed  from  a  relation¬ 
ship  between  the  lowest  positive  non-zero  and  the  record  highest  daily 
value  of  the  period  being  processed.  A  duration  table  is  produced  as 
output.  A  summary  is  produced  with  the  duration  table  which  gives  the 
value  of  each  class  limit,  the  number  of  days  in  the  period  in  each 
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class,  the  total  number  of  days  having  a  value  greater  than  or  equal  to 
each  class  limit,  and  the  percent  of  all  days  in  which  a  class  limit  is 
equaled  or  exceeded.  An  example  of  the  duration  table  and  summary  is 
shown  in  Figure  1 1 1-18. 

The  flow  duration  analysis  results  for  each  of  the  44  key  locations 
are  given  in  Table  A-l  of  Appendix  A.  Nine  percentages  are  chosen  for 
display  in  this  table.  The  flows  for  each  of  these  percentages  are  inter¬ 
polated  from  the  actual  class  limit  values  and  the  corresponding  percent¬ 
ages  given  in  the  duration  tables  developed  by  the  A969  program.  Flow 
duration  curves  for  the  naturalized  daily  flows  at  the  44  key  locations 
are  shown  on  Figures  A-l  to  A-44  in  Appendix  A.  These  figures  also  show 
the  duration  curves  for  the  regulated  flows  which  are  discussed  in 
Chapter  IV. 

The  low  flow  frequency  analysis  is  conducted  using  the  Log-Pearson 
Type  III  distribution  and  is  based  on  the  flow  values  of  the  low  flow 
frequency  table.  The  low  flow  table  contains  the  lowest  "n"-day  mean 
value  for  each  period  of  consecutive  days  in  each  water  year  be¬ 
gins  October  first.  The  table  also  includes  a  ranking  of  the  flows 
according  to  magnitude  in  the  period  of  years  analyzed.  An  example  of 
a  low  flow  table  produced  by  the  A969  program  is  shown  in  Figure  1 1 1- 19. 

The  Log-Pearson  low  flow  frequency  tables  from  the  A969  program 
contain  a  tabulation  of  the  input  data  from  the  low  flow  table,  the 
computed  statistics  of  flow  values  and  of  logarithms,  and  a  list  of 
eleven  non-exceedence  probabilities  with  corresponding  recurrence  inter¬ 
vals  and  the  theoretical  flow  values.  An  example  of  the  Log-Pearson 
low  flow  frequency  table  for  the  one-day  low  value  is  shown  in  Figure 
1 11-20. 


For  each  of  the  44  key  locations  low  flow  frequency  tables  are 
produced  for  periods  of  1,  3,  7,  14,  30,  60,  90,  120,  183  and  365  con¬ 
secutive  days.  The  naturalized  low  flow  frequency  results  are  given 
in  Table  A-2.1  to  A-2.44  of  Appendix  A  for  the  44  locations.  In  each 
table  the  low  flows  for  each  of  the  ten  consecutive  day  periods  are  given 
for  eleven  recurrence  intervals  and  the  corresponding  probabilities. 
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FIGURE  1 1 1 -1 8 .  Sample  Duration  Table  From  Program  A969 
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FIGURE  1 1 1 - 1 9 .  Sample  Low  Flow  Table  From  Program  A969 
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FIGURE  1 1 1-20.  Sample  Loq-Pearson  Low  Flow  Frequency  Table 
Table  From  Program  A969 


Low  flow  frequency  curves  for  the  7-day,  and  120-day  periods  are 
given  on  Figures  A-45  to  A-62  in  Appendix  A  for  18  key  locations.  These 
curves  indicate  for  a  given  magnitude  of  flow  the  recurrence  interval  in 
years  and  the  non-exceedence  probability  in  percent  for  each  of  the  two 
periods  of  n-consecutive  days. 

Duration  and  Frequency  Comparison  of  Observed  and  Simulated  Flows 

Several  key  locations  have  been  chosen  to  compare  the  duration 
and  frequency  curves  of  the  50  years  of  naturalized  flows  developed  by 
the  model  and  the  observed  flows  from  the  unregulated  period  of  record. 
These  locations  are  at  the  following  USGS  gaging  stations: 

East  Branch  Delaware  River  at  Fishs  Eddy,  N.Y., 

West  Branch  Delaware  River  at  Hale  Eddy,  N.Y., 

Neversink  River  at  Oakland  Valley,  N.Y., 

Delaware  River  at  Montague,  N.J.,  and 
Delaware  River  at  Trenton,  N.J. 

The  first  three  stations  are  downstream  from  the  three  New  York  City 
reservoirs  and  the  other  two  are  on  the  mainstem  of  the  Delaware. 

Figure  1 1 1 -21  presents  the  duration  curves  for  Fishs  Eddy,  Hale 
Eddy  and  Oakland  Valley,  and  Figure  1 1 1 -22  presents  the  duration  curves 
for  Montague  and  Trenton.  The  solid  lines  represent  the  curves  for  the 
50  years  of  naturalized  flows  and  the  dashed  lines  represent  the  curves 
for  the  unregulated  period  of  record  of  observed  flows.  Figures  1 1 1 - 23 
through  1 1 1-27  present  for  each  of  the  five  locations  the  7-day,  90-day 
and  365-day  low  flow  frequency  curves.  As  with  the  duration  curves  the 
solid  lines  represent  the  naturalized  flows  for  50  years  and  the  dashed 
lines  represent  the  observed  flows  for  the  unregulated  period  of  record. 

The  duration  curves  for  all  five  stations  show  that  for  the  lower 
flows  the  simulated  values  are  less  than  the  observed  values.  With  a 
few  exceptions  the  frequency  plots  also  give  the  same  results.  In  some 
instances  the  simulated  flows  are  less  than  the  observed  flows  for  the 
entire  range  of  return  intervals. 


These  results  reflect  differences  in  the  data  base  used  to  analyze 
the  frequency  and  duration.  The  unregulated  period  of  observed  flows 
is  less  than  50  years  for  which  naturalized  flows  have  been  developed. 

For  Fishs  Eddy  the  unregulated  period  of  records  is  26  years;  for 
Hale  Eddy,  35  years,  for  Oakland  Valley,  24  years;  for  Montague,  13 
years,  and  for  Trenton,  25  years.  One  of  the  major  factors  in  the  dif¬ 
ference  is  the  fact  that  the  observed  flows  do  not  include  the  drought  of 
the  1960's.  The  naturalized  flow  model  simulates  the  50  years  of  flow 
based  on  the  correlation  and  fill-in  analysis  described  in  Chapter  II. 

The  simulated  flows  include  the  drought  of  the  1960 's  and  several  years 
of  low  flow  periods  are  input  into  the  frequency  and  duration  analysis. 
This  produces  curves  which  when  compared  to  the  observed  curves  have  lower 
flows  for  the  given  return  intervals  on  the  frequency  plots  and  lower 
flows  for  the  percent  the  flow  was  equaled  or  exceeded  on  the  duration 
plots. 
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4-ui  FIGURE  III-23.  COMPARISON  OF  NATURALIZED 
AND  OBSERVED  LOW  FLOW  FREQUENCY 
CURVES  FOR  THE  E.  BR.  DELAWARE 
RIVER  AT  FISHS  EDDY 
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FIGURE  1 1 1-24.  COMPARISON  OF  NATURALIZED 
AND  OBSERVED  LOW  FLOW  FREQUENCY  CURVES 
FOR  THE  W.  BR.  DELAWARE  RIVER  AT  HALE  EDDY 
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IV.  REGULATED  FLOW  MODEL 


PURPOSE 


The  purpose  of  the  regulated  flow  model  is  to  determine  the  effect 
of  New  York  City  reservoir  projects  on  the  low  flow  frequency  and  flow 
duration  analyses  at  key  locations  in  the  Delaware  River  Basin.  Because 
this  study  is  to  analyze  the  low  flow  periods,  sophisticated  routing 
techniques  are  not  required.  The  emphasis  instead  is  on  the  ability  of 
the  model  to  match  the  desired  flow  condition  within  the  constraints  of 
the  reservoir  operating  rules. 

INCREMENTAL  RUNOFF  AND  ROUTING  PROCEDURES 

The  incremental  inflow  which  is  used  in  the  natural  flow  model  is 
also  used  in  the  regulated  flow  model  with  one  major  exception.  The 
natural  inflows,  developed  with  the  correlation  and  extension  procedures, 
for  the  three  New  York  City  reservoirs,  Pepacton,  Cannonsville  and  Never- 
sink,  are  replaced  with  daily  net  inflows  to  the  reservoirs  during  the 
regulated  period  of  record.  The  daily  net  inflows,  obtained  from  the 
New  York  City  Bureau  of  Water  Supply  files  of  Daily  Yields  in  the  Delaware 
Watershed,  take  into  consideration  the  precipitation  and  evaporation  over 
the  reservoir  surface.  The  periods  for  which  the  reservoir  inflow  data 
are  available  are  given  in  Table  IV- 1 . 

The  net  evaporation  for  the  New  York  City  reservoirs  during  the 
unregulated  period  of  record  is  accounted  for  in  the  regulated  flow  model 
by  subtracting  the  net  evaporation  from  the  reservoir  storage.  The  net 
evaporation  which  is  simulated  in  the  daily  model  is  based  on  monthly 
averages  of  evaporation  and  precipitation.  Evaporation  data  are  obtained 


A 


TABLE  IV- 1 

PERIODS  OF  RECORDED  NET  INFLOW  DATA 
FOR  THE  NEW  YORK  CITY  RESERVOIRS* 


Reservoir 

Period  of  Record 

Pepacton 

9/17/54-5/31/55 

7/1/55-9/30/77 

Cannons ville 

10/1/63-9/30/77 

Neversink 

7/1/53-12/31/53 

2/1/54-7/31/59 

1/1/60-9/30/77 

*Periods  are  those  which  were  used  in  the  regulated  flow 
model . 


from  the  Climatic  Atlas  of  the  United  States*.  The  average  annual  lake 
evaporation  for  the  Delaware  River  Basin  in  the  region  of  the  reservoirs 
is  divided  into  two  six-month  periods  with  75  percent  of  the  average  annual 
evaporation  occurring  in  the  months  May  through  October.  Average  monthly 
precipitation  data  are  obtained  from  the  Report  of  the  River  Master  of 
the  Delaware  River^.  The  monthly  averages  are  for  the  Delaware  River 
Basin  above  Montague,  New  Jersey  for  the  period  December  1940  to 
November  1975. 

The  flow  correlation  and  extension  procedures  for  the  unregulated 
period  of  record  account  for  the  runoff  over  the  area  of  the  reservoir 
site.  In  order  to  avoid  double  accounting  of  that  percentage  of  the 
precipitation  which  occurs  as  runoff,  the  net  evaporation  is  calculated 
as  the  evaporation  minus  one-half  precipitation.  This  analysis 
assumes  that  the  direct  runoff  is  50  percent  of  the  precipitation.  Table 
IV-2  gives  the  monthly  average  e vaporat ■ on ,  one-half  the  precipitation 
and  net  evaporation  which  is  used  the  regulated  flow  model  for  the 
unregulated  period  of  record.  The  simulation  of  the  net  evaporation 
is  also  used  for  the  periods  of  missing  daily  net  inflows  to  the  reser¬ 
voirs  during  the  regulated  period  of  record  as  shown  in  Table  IV-1  for 
Pepacton  and  Neversink  reservoirs. 

Identical  routing  procedures  are  used  in  the  regulated  flow  model 
as  are  used  in  the  natural  flow  model. 


*  Climatic  Atlas  of  the  United  States,  U.S.  Department  of  Commerce, 
Environmental  Science  Services  Administration,  June  1968. 

2 

Report  of  the  River  Master  of  the  Delaware  River,  for  the  Period 
December  1,  1975  -  November  30,  1976,  United  States  Department  of 
the  Interior,  Geological  Survey  National  Center,  Reston,  Virginia, 
1977. 
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TABLE  IV-2 


NET  EVAPORATION  FOR  NEW  YORK  CITY  RESERVOIRS 
(INCHES/MONTH) 

Month 

Evaporation* 

2 

One-half  Precipitation 

Net  Evaporation 

Jan 

1.208 

1.400 

-0.192 

Feb 

1.208 

1.375 

-0.167 

Mar 

1.208 

1.620 

-0.412 

Apr 

1.208 

1.800 

-0.592 

May 

3.625 

2.075 

1.550 

Jun 

3.625 

1.995 

1.630 

Jul 

3.625 

2.080 

1.545 

Aug 

3.625 

2.010 

1.615 

Sep 

3.625 

1.800 

1.825 

Oct 

3.625 

1.555 

2.070 

Nov 

1.208 

1.920 

-0.712 

Dec 

1.208 

1.795 

-0.587 

Total s 

28.998 

21.425 

7.573 

Source:  Climatic  Atlas  of  the  United  States,  U.S.  Department  of  Commerce, 
Environmental  Science  Services  Administration,  June  1968  (75%  of 
„  annual  occurs  in  May  through  October). 

c  Source:  Report  of  the  River  Master  of  the  Delaware  River,  For  the  Period 
December  1,  1975  -  November  30,  1976,  United  States  Department 
of  the  Interior,  Geological  Survey  National  Center,  Reston, 
Virginia,  1977. 
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NEW  YORK  CITY  RESERVOIR  OPERATING  PROCEDURE 


The  operation  rules  for  diversion  and  releases  from  Cannonsville, 
Pepacton  and  Neversink  Reservoirs  are  as  follows: 


•  The  total  daily  New  York  City  water  supply  diversion  is  with¬ 
drawn  from  only  one  of  the  three  reservoirs  each  day.  The 
criterion  for  selection  is  based  on  the  reservoir  which  is  the 
fullest  in  percent  of  total  volume. 

•  Basic  conservation  releases  are  made  on  each  day  from  all  three 
reservoirs.  Additional  releases  are  required  when  the  flow  at 
Montague  does  not  equal  or  exceed  a  prescribed  target.  These 
required  releases  are  always  made  from  the  percent  fullest 
reservoir. 

•  The  priority  of  water  use  and  its  effects  are  as  follows.  The 
water  supply  diversion  is  taken  first  from  the  percent  fullest 
reservoir.  After  the  diversion  from  the  reservoir  has  been 
accounted  for,  the  Montague  release  is  made  (if  necessary)  from 
the  percent  fullest  reservoir.  The  reservoir  which  releases  to 
meet  the  flow  objective  at  Montague  can  be  a  different  reser¬ 
voir  than  the  one  from  which  the  water  supply  is  diverted.  Thus 
for  each  day,  accounting  for  diversions  and  releases,  the  reser¬ 
voir  volumes  are  kept  in  balance. 


The  target  flow  at  Montague  and  the  maximum  diversion  allowed  from 
the  New  York  City  reservoirs  is  set  by  the  operating  level  of  the  reser¬ 
voirs.  The  New  York  City  reservoirs  are  either  in  normal,  drought  warning, 
or  drought  condition  depending  on  the  combined  storage  in  the  reservoirs 
each  day  and  the  time  of  year.  Figure  I V- I  shows  the  reservoir  storage 
curves  which  defined  the  reservoir  conditions  and  the  maximum  diversions 
and  Montague  targets  for  each  condition  developed  for  this  study.  The 
storage  curves,  diversions,  and  targets  are  a  modification  of  the  Drought 
Emergency  Operating  Rules  -  Model  #2  for  the  New  York  Delaware  System 
presented  in  the  Task  Group  Report,  DRBC  Docket  No.  D-77-201. 


The  basic  conservation  releases  for  each  reservoir  vary  with  the 
season.  Table  IV-3  presents  the  basic  conservation  releases  for  the  New 
York  City  reservoirs.  The  conservation  releases  are  those  specified  in 
the  DRBC  Docket  D77-20CP . 

*  Task  Group  Report  DRBC  Docket  No.  D-77-20  Appraisal  of  Upper  Basin 
Reservoir  Systems,  Drought  Emergency  Criteria  and  Conservation  Measures, 
Delaware  River  Basin  Commission,  March  1979. 
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FIGURE  IV- 1 .  Modified  Drought  Emergency  Operating  Rules-Model  #2 

New  York  Delaware  System 


The  operating  rules  for  diversion  and  releases  used  in  the  daily 
flow  model  simulate  the  flow  response  in  the  basin  under  the  conditions 
that  all  three  reservoir  volumes  are  kept  in  balance  with  diversions  and 
releases  from  the  percent  fullest  reservoir.  Also,  the  maximum  allowable 
diversions  are  withdrawn  under  each  of  the  three  reservoir  conditions.  The 
model  approach,  which  evaluates  the  total  storage  available  at  the  maximum  di 
version  rates,  is  somewhat  different  than  the  actual  operation  of  the 
New  York  City  reservoirs.  In  the  actual  operation,  the  maximum  diversions 
are  not  always  withdrawn  from  the  three  reservoirs.  The  actual  amount 
depends  on  the  overall  diversion  schedule  from  the  entire  New  York  City 
water  supply  system  of  which  the  three  reservoirs  in  the  Delaware  River 
Basin  are  only  a  part.  Also,  instead  of  withdrawing  the  diversions  from 
the  percent  fullest  reservoir  and  then  releasing  to  meet  the  Montague 
objective  from  the  percent  fullest  reservoir,  the  actual  schedule  is  as 
follows.  The  Montague  target  is  checked  and  conservation  releases  are 
made  from  Pepacton  and  Neversink  reservoirs.  The  deficit  at  Montague 
is  made  up  by  releases  from  Cannons- • >le  Reservoir.  Then  the  New  York 
City  diversions  are  taken  from  all  three  reservoirs  maintaining  them  in 
equal  balance.  The  actual  system  also  has  limitations  on  the  maximum  amount 
of  water  which  can  be  released  from  any  one  reservoir  based  on  the  outlet 
capacities.  This  restriction  is  not  imposed  on  the  model  scheme. 

MODEL  SIMULATION  OF  RESERVOIR  SYSTEM 

The  "egulated  flow  model  is  unable  to  predict  in  advance  whether 
the  flow  at  Montague  will  be  below  the  target  flow.  Therefore  if  the 
flow  ac  Montague  falls  below  the  target  on  any  given  day,  the  model 
"backs  up"  two  days,  releases  the  required  amount  from  a  reservoir,  and 
proceeds  with  the  simulation  from  that  day.  The  flow  at  Montague  in  two 
days  may  not  exactly  equal  the  target  flow.  The  releases  from  the  New 
York  City  reservoirs  must  be  routed  for  approximately  two  days  down  to 
Montague.  The  routing  procedure  with  a  daily  time  step  attenuates  the 
release  over  at  least  three  days.  But  enough  water  is  released  on  one 
day  so  that  the  three-day  average  at  Montague  is  always  nearly  equal  to 
or  greater  than  the  target  flow.  The  attenuation  of  the  release  simulates 
the  attenuation  which  occurs  in  the  real  river  system.  Using  a  smaller 
time-step  would  better  simulate  attenuation  and  the  daily  flows  in  the 
basin,  but  the  magnitude  of  additional  costs  would  not  be  justified  by 
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the  increased  accuracy.  As  the  model  currently  works,  the  three-day 
average  flow  at  Montague  is  simulated  accurately. 


The  backup  procedure  is  based  on  a  simple  concept.  The  model 
simply  backs  up  two  days  and  operates  with  the  values  of  all  variables 
at  that  point  in  time.  The  model  makes  the  required  release  from  one  of 
the  New  York  City  Reservoirs  in  order  to  meet  the  target  flow  at  Montague 
two  days  later. 

Figure  IV-2  shows  the  flow  chart  of  the  reservoir  operating  rules, 
which  are  used  in  the  regulated  flow  model. 

A  review  of  the  model  output  clarifies  the  operating  rules.  Figures 
IV-3,  IV-4  and  I V- 5  show  the  percent  of  remaining  storage  in  each  of  the 
three  New  York  City  reservoirs  in  calendar  year  1953  of  the  simulation. 

The  releases  and  withdrawals  from  the  reservoir  are  modeled  to  maintain  an 
equal  balance  of  the  percent  reservoir  storage  for  all  three  reservoirs. 

The  computer  output  of  the  three  figures  show  that  the  percent  remaining 
storages  are  approximately  equal  for  the  Cannonsvil le,  Pepacton  and 
Neversink  reservoirs.  For  example,  in  January  1953  for  Cannonsvil le 
(Figure  IV-3),  Pepacton  (Figure  IV-4)  and  Neversink  (Figure  IV-5),  the 
minimum  percent  storages  remaining  are  72.80,  71.95  and  71.42;  the  maxi¬ 
mum  percent  storages  remaining  are  87.31,  86.97,  and  86.36;  and  the  mean 
percent  storages  remaining  are  76.89,  76.48  and  75.43. 

Figures  IV-6,  IV- 7  and  IV-8  show  the  total  outflow  from  each  of  the 
three  reservoirs.  Basic  conservation  releases  are  always  made  from  ear.1 
reservoir,  required  releases  to  Montague  are  always  made  from  the  percent 
fullest  reservoir,  and  spills  occur  when  the  reservoirs  are  full.  The  daily 
variations  of  releases,  from  low  to  high  flows,  based  on  the  model  scheme 
does  not  necessarily  reflect  the  actual  reservoir  operations. 

For  example,  conservation  releases  are  made  from  each  of  the  three 
reservoirs  during  the  entire  months  of  January  and  February,  1953  with  no 
spills  occurring.  A  check  against  the  reservoir  storages  during  these  two 
months,  presented  in  Figures  IV-3,  IV-4  and  IV-5,  show  that  the  reservoirs 
are  less  than  full.  These  figures  also  demonstrate  that  releases  are  made 
from  the  percent  fullest  reservoir.  For  31  August  1953  a  release  of 
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ENTER  RESOR 


FIGURE  IV-2  RESERVOIR  OPERATING  RULES  AS 
EXECUTED  IN  THE  MODEL 
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FIGURE  I V- 3 .  Percent  of  Remaining  Available  Storage  in 
Cannonsville  Reservoir  for  1953 


FIGURE  IV-4.  Percent  of  Remaining  Available  Storage 
in  Pepacton  Reservoir  for  1953 


FIGURE  IV-5.  Percent  of  Remaining  Available  Storage 
in  Neversink  Reservoir  for  1953 
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FIGURE  IV-7.  Total  Outflow  from  Cannonsville  Reservoir  in  1953 


i 


425.82  cfs  is  made  from  the  Pepacton  Reservoir  (Figure  I V-7)  while  the 
Cannonsville  Reservoir  (IV-6)  and  the  Neversink  Reservoir  (Figure  IV-8) 
are  making  basic  conservation  releases  of  23.0  cfs  and  15.5  cfs  respective¬ 
ly.  For  30  August  1953,  which  is  checked  to  determine  the  next  day 
release  from  the  percent  fullest  reservoir,  the  percent  remaining  storage 
of  Pepacton  (66.57  shown  in  Figure  IV-4)  is  greater  than  the  percent 
remaining  storage  of  Cannonsville  (65.60  shown  in  Figure  IV-3)  and  Never¬ 
sink  (65.46  shown  in  Figure  IV-5). 

For  certain  dates  these  figures  show  that  a  release  from  Neversink 
occurs  on  the  same  day  as  a  release  from  one  of  the  other  two  reservoirs. 
The  model  is  programmed  to  release  the  Neversink  flow  augmentation  for 
Montague  one  day  later  than  the  day  for  which  the  release  is  calculated 
to  occur.  Because  the  route  time  from  Neversink  to  Montague  is  approx¬ 
imately  one  day  less  than  the  route  time  of  the  other  two  reservoirs,  this 
adjustment  is  made  to  make  certain  that  only  one  release  flow  from  any  of 
the  three  reservoirs  arrives  at  Montague  on  any  one  day. 

Figure  IV-9  shows  the  amount  of  water  diverted  to  New  York  City. 

Note  that  in  1953  the  model  is  in  a  normal  condition  for  most  of  the  year 
during  which  800  mgd  is  diverted  to  New  York  City.  During  the  period 
October  24  through  December  7  the  model  is  in  a  drought  warning  condition 
and  600  mgd  is  diverted  to  New  York  City. 

The  final  figures  in  this  series  show  the  flow  at  Montague  for  the 
natural  model  (Figure  IV-10)  and  the  regulated  model  (Figure  IV-11). 

Figure  IV-10  shows  that  for  many  days  especially  during  the  months  of  July 
through  October  the  natural  flow  at  Montague  is  less  than  the  Montague 

target  used  in  the  regulated  model  simulation.  The  regulated  flow  model 
operates  to  maintain  the  Montague  target.  However,  Figure  IV-11  shows 
several  days  which  fall  slightly  below  the  target  of  1750  cfs.  This 
variation  is  due  to  the  attenuation  of  the  releases  over  three  or  more 
days  which  is  built  into  the  model.  Releases  are  attenuated  as  they 
travel  from  the  reservoirs  to  Montague.  The  amount  of  attenuation  of  the 
releases  is  different  for  each  reservoir  depending  on  the  total  time  of 
travel  from  the  reservoir  to  Montague.  If  releases  are  made  for  several 
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FIGURE  IV-9.  Average  Daily  Diversion  (mgd)  to  New  York 
City  Simulated  for  1953 
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FIGURE  IV-10.  Natural  Flow  at  Montague  Simulated 

for  1953 
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FIGURE  IV- 1 1 .  Regulated  Flow  at  Montague 
Simulated  for  1953 


days  in  a  row  (as  in  August  1953),  and  are  made  from  one  reservoir  on 
one  day  and  another  reservoir  on  the  next,  the  discharges  at  Montague 
fluctuate  about  the  target  flow  depending  on  how  the  routed  releases  have 
combined.  However,  the  three  day  average  of  the  discharges  at  Montague 
always  nearly  equals  or  exceeds  the  target  flow  as  discussed  previously 
in  this  chapter. 

RESULTS 


The  question  addressed  in  this  aspect  of  the  modelling  project  is, 
"Can  the  upper  basin  reservoirs  provide  enough  water  to  meet  the  demands 
both  at  New  York  City  and  Montague?"  As  shown  in  the  natural  model  results 
in  Chapter  III,  the  periods  of  concern  are  the  droughts  of  the  thirties  and 
the  sixties  or  more  specifically  the  most  critical  years  1931  and  1965. 

The  regulated  flow  model  further  substantiates  the  fact  that  the  sixties 
drought  was  the  most  severe. 

The  Drought  of  the  1930's 

To  analyze  the  drought  of  the  thirties,  two  operating  rules  are 
investigated.  The  first  rule  is  the  modified  Drought  Emergency  Operating 
Rules-Model  #2  which  is  discussed  earlier  in  this  report.  Figures  IV-12, 

I V- 1 3  and  I V- 14  show  the  remaining  useable  storage  in  each  of  the  three 
reservoirs  in  1931  using  this  operating  rule.  The  minimum  daily  percent 
remaining  storage  for  Cannonsville  (Figure  IV-12)  is  23.98;  for  Pepacton 
(Figure  I V - 1 3 ) ,  23.99;  and  for  Neversink  (Figure  IV-14),  23.01  all  of 
which  occur  during  the  month  of  February.  Therefore,  there  is  enough 
reservoir  storage  in  1931  to  maintain  the  flow  target  at  Montague. 

Figures  I V- 1 5  and  I V- 16  list  for  1930  and  1931,  respectively,  the  average 
daily  diversion  to  New  York  City.  Figure  I V- 1 5  shows  that  on  November  27, 
1930  the  reservoir  system  is  shifted  into  the  Drought  Warning  Condition. 

At  this  time  the  diversion  to  New  York  City  is  reduced  from  800  mgd  to 
600  mgd  and  the  required  flow  at  Montague  remains  at  1,750  cfs.  Figure 
I V- 16  shows  that  on  January  3,  1931  the  system  goes  into  the  Drought 
Condition.  This  condition  provides  430  mgd  to  New  York  City  and  maintains 
1,525  cfs  at  Montague. 


IV-21 
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FIGURE  IV-12.  Percent  of  Remaining  Usable  Storage  in  Cannonsville 
Reservoir  for  1931  under  modified  Operating  Rules  #2 
(95.7  BG  Total  Storage) 
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FIGURE  I V- 13 .  Percent  of  Remaining  Usable  Storage  in  Pepacton 

Reservoir  for  1931  under  modified  Operating  Rule  #2 
(140.2  BG  Total  Storage) 
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FIGURE  IV-15.  Average  Daily  Diversion  to  New  York  City 
for  1930  under  modified  Operating  Rules  #2 
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A  second  set  of  operating  rules  is  also  tested  for  the  thirties. 

This  second  case  always  provides  800  mgd  to  New  York  City  and  also  main¬ 
tains  1,750  cfs  at  Montague  at  all  times.  These  conditions  are  those 
established  in  the  1954  Supreme  Court  Decree^. 

It  should  be  noted  that  the  model  is  designed  to  guarantee  the 
required  flow  to  both  places.  Thus,  to  preserve  continuity  the  storages 
of  the  reservoirs  are  allowed  to  become  negative  if  not  enough  water  is 
available  for  both  places.  The  advantage  of  this  approach  is  that  the 
user  can  determine  the  amount  of  additional  storage  required  to  overcome 
any  deficiency. 

The  1954  Supreme  Court  Decree  conditions  are  tested  with  the  equal 
drawdown  of  the  New  York  City  reservoirs.  That  is,  the  diversions  and 
releases  are  taken  from  the  percent  fullest  reservoir.  When  the  reser¬ 
voirs  are  drawndown  evenly,  the  minimum  reservoir  storage  during  the  1930's 
drought  is  about  26  billion  gallons  (9.5  percent  of  total  useable  storage). 
There  is  enough  storage  remaining  after  meeting  the  1954  Supreme  Court  Decree 
conditions  to  provide  an  additional  diversion  of  39  mgd  to  New  York  City 
during  the  1930‘s  drought.  Therefore,  with  equal  drawdown  of  the  New  York 
City  reservoirs,  the  maximum  diversion  rate  of  the  1930‘s  drought  is  839  mgd 
with  a  1,750  cfs  flow  objective  at  Montague.  This  result  compares  favorably 
with  the  results  of  the  River  Master3  and  the  DRBC  (855  mgd  and  848  mgd, 
respectively)  presented  in  the  DRBC  Task  Group  Report  DRBC  Docket  No. 

3 

D-77-20  .  The  maximum  diversion  rate  of  839  is  estimated  with  a  daily 
flow  model,  and  with  the  attenuation  of  reservoir  releases  to  Montague. 

The  River  Master  and  DRBC’s  results  are  estimated  with  no  attenuation  of 
releases  to  Montague,  and  with  a  monthly  analysis  for  some  periods.  The 
River  Master  and  DRBC  also  estimated  flows  at  Montague  and  inflows  to 
Pepacton,  Cannonsville,  and  Neversink  Reservoirs  with  correlations  to 


United  States  Supreme  Court  Decree,  New  Jersey  v.  New  York,  347  U.S.  995 
(1954),  Approved  by  the  United  States  Supreme  Court  June  7,  1954. 

O 

Office  of  the  Delaware  River  Master,  U.S.  Geological  Survey  National 
Center,  Reston,  Virginia. 

3  Task  Group  Report  DRBC  Docket  No.  D-77-20,  Appraisal  of  Upper  Basin 
Reservoir  Systems,  Drought  Emergency  Criteria  and  Conservation  Measures, 
Delaware  River  Basin  Commission,  March  1979. 
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different  USGS  stations  than  in  the  regulated  flow  model.  The  variation 
between  results  is  no  surprise  when  all  of  the  above  factors  are  considered. 

The  Drought  of  the  1960 ' s 

The  regulated  model  operates  to  provide  the  water  required  for 
the  New  York  City  diversion  and  for  the  Montague  target  allowing  the  reser¬ 
voir  storages  to  go  negative.  During  November  of  1965,  the  worst  year 
of  the  1960's  drought,  the  modeled  reservoir  storages  are  barely  negative. 

Thus,  the  actual  reservoir  storages  are  almost  able  to  provide  enough  water  for 
the  New  York  City  diversion  and  the  target  flow  at  Montaaue  with  the  modified 
Drought  Emergency  Operating  Rules-Model  #2  and  equal  drawdown.  Figures 
IV-17,  IV- 18  and  IV- 19  present  the  percent  storage  remaining  in  1965  for 
Cannonsville,  Pepacton  and  Neversink  Reservoirs,  respectively.  For  the 
month  of  November,  the  minimum  daily  percent  storage  remaining  for  Cannons¬ 
ville  is  -1.68  (Figure  IV-17);  for  Pepacton,  -1.78  (Figure  IV- 18) ;  and 
for  Neversink,  -2.74  (Figure  IV-19).  November  1965  is  the  only  month 
in  the  sixties  drought  at  which  time  the  New  York  City  diversion  and  the 
Montague  flow  target  cannot  be  met.  During  this  month  and  the  other 
months  of  1965  the  reservoir  system  is  operating  under  the  drought  condi¬ 
tions  of  the  modified  Drought  Emergency  Operating  Rules-Model  #2.  The 
Montague  target  under  the  drought  condition  is  1525  cfs.  The  diversion 
from  the  New  York  City  reservoirs  is  430  mgd  which  is  shown  in  Figure 
IV-20,  Average  Daily  Diversions  to  New  York  City  for  1965  under  Modified 
Operating  Rules  #2.  These  and  other  outputs  indicate  that  the  reservoir 
system  had  been  in  a  drought  condition  for  over  one  year  before  the  reser¬ 
voir  emptied. 

The  1954  Supreme  Court  Decree  Rules  are  also  tested  with  even 
drawdown  during  the  1960's  drought,  and  as  expected,  the  800  mgd 
diversion  and  1,750  cfs  Montague  flow  objective  cannot  be  maintained 
with  the  current  New  York  City  reservoir  storage.  The  critical  period 
of  the  1960 ' s  drought  was  June  1964  through  November  1965.  If  the 
reservoirs  were  full  on  June  1964,  the  maximum  rate  of  diversion  which 
could  be  sustained  during  the  critical  period  is  466  mgd.  The  River 
Master  and  the  City  of  New  York  computed  maximum  rates  of  diversion  for 
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FIGURE  IV-18.  Remaining  Usable  Storage  in  Pepacton 

Reservoir  for  1965  under  modified  Operating  Rules  #2 
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FIGURE  IV-20.  Average  Daily  Diversion  to  New  York  City 

for  1965  under  Modified  Operating  Rules  #2 


this  period  of  482  and  482.3  mgd,  respecti vely1 .  Both  the  River  Master 
and  New  York  City  assumed  no  attenuation  of  the  releases  made  to  Montague, 
which  accounts  for  the  slight  discrepancy  in  results. 

DURATION  AND  FREQUENCY  ANALYSIS 

Flow  duration  and  low  flow  frequency  analyses  are  performed  on  the 
simulated  50  years  of  regulated  daily  flows  for  44  key  locations  which 
were  also  analyzed  for  the  naturalized  flows.  The  U.S.  Geological  Survey's 
Daily  Values  Statistics  Program  (A969)  as  described  in  Chapter  III,  is 
used  to  perform  the  analysis.  For  the  regulated  flow  simulation  the  A969 
program  was  run  by  the  Philadelphia  Corps  of  Engineers  at  Boeing  Computer 
Service  in  Vienna,  Virginia.  The  tables  and  frequency  plots  of  the  re¬ 
gulated  daily  flows  are  given  in  Appendix  B. 

The  duration  tables  for  each  location  are  given  on  Table  B-l  of 
Appendix  B.  The  flow  values  equaled  or  exceeded  for  nine  percentages  are 
displayed  in  the  duration  table.  Reservoir  regulation  of  the  flows 
causes  a  variation  of  the  frequency  curves  from  those  which  were  produced 
by  the  naturalized  flows.  The  duration  curves  are  presented  in  Figures  A-l 
to  A-44  in  Appendix  A.  The  curves  are  plotted  on  the  same  graph  as  the 
naturalized  flow  duration  curves. 

On  the  mainstem  of  the  Delaware  River  below  the  reservoirs,  the 
regulation  of  flow  causes  the  higher  flows  of  naturalized  conditions  to 
be  dampened  out  and  reduced.  It  also  increases  the  lower  naturalized 
flow  as  flow  requirements  at  Montague  are  being  met  by  releases  from  the 
reservoirs.  Figures  A-5,  A-14,  A-29,  A-44  of  Appendix  A  demonstrate  this 
fact  for  the  locations  at  Callicoon,  Montague,  Trenton,  and  the  Delaware 
Memorial  Bridge  at  Wilmington,  Delaware,  respecti vely .  The  square  symbols 
on  the  Figures  represent  the  natural  duration  curves  for  these  locations 
and  the  circular  symbols  represent  the  regulated  duration  curves. 


‘Task  Group  Report  DRBC  Docket  No.  D-77-20,  Appraisal  of  Upper  Basin 
Reservoir  Systems,  Drought  Emergency  Criteria  and  Conservation  Measures, 
Delaware  River  Basin  Commission,  March  1979. 
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For  each  of  the  ^4  key  location  Loq-Pearson  low  flow  frequency 
tables  are  produced  for  periods  of  1,  3.  7,  14,  30,  60,  99,  129,  183 
and  365  consecutive  days.  The  regulated  low  flow  freauency  results 
are  given  in  Tables  B-2.1  through  B-2.44  in  Appendix  B.  In  each 
table  the  low  flov/s  for  each  of  the  tpn  consecutive  hay  periods  are 
given  for  eleven  recurrence  intervals  and  the  corresponding 
probahi 1 ities. 

The  l.og-Pearson  distribution  is  a  cumulative  probability 
distribution  determined  from  the  base-10  logarithms  of  the  sample 
observations.  These  '-'nservations  are  the  low  1  day,  7  day,  etc.  for 
each  of  the  years  analyzed.  Pecause  of  the  regulation  of  flows  to  meet 
the  flow  requirements  at  Montague  the  lowest  flows  for  the  69-year 
period  are  similar  in  magnitude  for  the  60,  90,  129  and  133  consecutive 
day  periods  for  all  the  mainstem  locations  below  the  reservoirs.  ror 
example  the  lowest  90-day  flow  for  the  60-year  period  at  Montague  is 
1520  cfs,  the  lowest  120-day  flow  is  1530  cfs,  and  the  lowest  183-day 
flow  is  1530  cfs. 


For  several  locations  on  the  mainstem  below  the  reservoirs  the 
calculated  low  flow  frequency  values  at  the  lower  probabilities  do  not 
consistently  increase  with  the  increasinu  number  of  consecutive  days  as 
should  be  expected.  For  example,  in  some  casps  the  129-riay  flow  for  the 
109  year  recurrence  interval  is  calculated  to  be  less  t^an  the  °0-day 
flow  for  the  190  year  recurrence  interval.  This  is  a  result  of  the  I  no- 
Pearson  distribution  fit  through  the  sample  data  in  whirh  the  lowest 
flows  are  similar  for  each  of  these  two  consecutive  day  periods  but  the 
higher  flows  of  the  129  day  period  are  greater  than  the  high  Hows  for 
the  °0  day  period.  This  causes  the  frequency  curvp  for  the  129  day 
period  to  actually  fall  below  the  frequency  curve  for  the  90  day  period. 
Whenever  these  inconsistencies  occurred,  the  flov/s  are  plotted  against 
duration  on  log-log  paper  for  the  particular  probability.  Adjustments 
are  made  so  that  low  flow  always  increases  with  duration. 
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Regulated  low  flow  frequency  curves  for  the  7-day  and  120-day 
periods  are  shown  in  Figures  B-l  through  B-18  in  Appendix  B  for  the  18 
key  locations  given  in  Table  IV-4.  The  square  symbols  on  the  figures 
represent  the  7-day  frequency  curves  and  the  circular  symbols  represent 
the  120-day  frequency  curves.  The  frequency  curves  of  the  natural  flows, 
given  in  Appendix  A  show  a  wider  variation  of  the  flows  for  given  recurrence 
intervals  than  do  the  frequency  curves  at  locations  where  regulation  occurs. 
The  flows  for  the  7-day  and  120-day  periods  for  the  regulated  flows  are 
similar  in  magnitude  during  low  flow  conditions  due  to  the  flow  regulation 
by  the  reservoirs  to  maintain  higher  flows  in  the  Delaware  River  during 
dry  weather  periods. 
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TABLE  IV- 4 


KEY  LOCATIONS  FOR  LOW  FLOW 
FREQUENCY  CURVES 


Location 


East  Branch  Delaware  River 
at  Downsville,  NY 

West  Branci  Delaware  River 
at  Stilesvi lie,  NY 

West  Branch  Delaware  River 
at  Hale  Eddy,  NY 

Delaware  River  near 
Cal li coon,  NY 

Delaware  River  near 
Barryville,  NY 

Lackawaxen  River  at 
Hawley.  PA 

Delaware  River  at 
Port  Jervis,  NY 

Neversink  River  at 
Neversink,  NY 

Delaware  River  at 
Montague,  NJ 

Pohopoco  Creek  at 
Beltzville  Damsite,  PA. 

Lehigh  River  at 
Bethlehem,  PA 

Tohickon  Creek  at 
Pipersville,  PA 

Delaware  River  at 
Trenton,  NJ 

Tulpehocken  Creek  at 
Blue  Marsh  Damsite,  PA. 

Schuylkill  River  at 
Reading,  PA 

Schuylkill  River  at 
Philadelphia,  PA 

Delaware  River  below 
Mouth  of  Schuylkill 

Delaware  River  at 
Delaware  Memorial  Bridge 


USGS  Station 
01417000 

01425000 

01426500 

01427405 

01428500 

01431500 

01434000 

01436000 

01438500 

01449800 

01453000 

01459500 

01463500 

01470960 

01471500 

91474500 

None 

None 
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I.  INTRODUCTION 


The  purpose  of  Phase  II  is  to  analyze  various  reservoir  combinations 
in  the  Delaware  River  Basin,  used  specifically  to  supplement  the  flow  of  the 
Delaware  River  at  Trenton,  New  Jersey  during  low  flow  conditions.  Combina¬ 
tions  include  both  proposed  and  existing  reservoirs,  with  modifications  to 
some  of  the  existing  reservoirs'  operating  rules.  The  following  reservoirs 
are  systematically  analyzed  for  augmenting  low  flow  at  Trenton:  Beltzville, 

FE  Walter,  Prompton.Nockamixon,  and  Hackettstown.  Two  other  reservoirs, 
Merrill  Creek  and  Cannonsville  with  a  modified  storage,  are  also  analyzed 
even  though  they  do  not  release  directly  to  maintain  the  Trenton  flow.  Table 
1-1  lists  all  the  reservoirs  involved  in  the  analysis  along  with  pertinent 
information  of  each:  location,  USGS  gaging  station,  drainage  area  and  usable 
storage.  All  other  reservoirs  in  the  basin  are  operated  as  they  have  been 
in  the  past.  The  New  York  City  reservoirs  of  Pepacton,  Cannonsville  and 
Neversink  are  operated  according  to  specified  operating  rules  discussed  in 
Chapter  II. 

With  reservoirs  on  line,  additional  storage  and  modifications  to 
operating  rules  become  an  increased  regulated  supply  of  water  which  will  un¬ 
doubtedly  improve  the  low  flow  characteristics  of  the  basin.  The  question 
addressed  here  is  which  reservoir  combinations  are  more  beneficial  than 
others.  Table  1-2  gives  the  17  reservoir  combinations  analyzed  in  Phase  II. 
The  analysis  tests  the  ability  of  each  combination  to  sustain  a  maximum  low 
flow  objective  at  Trenton  for  a  series  of  historical  hydrologic  conditions 
such  as  the  worst  event  in  the  past  50  years  and  an  average  event  in  the  same 
50  years.  Note  that  any  combination,  which  includes  Prompton,  routes  Prompton' 
releases  around  Montague,  using  a  methodology  which  at  the  same  time  gives 
Montague  credit  for  Prompton's  natural  inflow.  A  detailed  discussion  is 
given  in  Chapter  III. 

The  maximum  low  flow  objective  is  achieved  when  the  total  usable 
storage  of  the  reservoirs  as  given  in  Table  1-1  is  depleted  with  no  considera 
tion  of  reservoir  or  release  water  quality.  Each  reservoir  does  maintain  a 
minimum  pool  level  that  is  not  considered  a  part  of  the  usable  storage. 
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TABLE  1-1 


RESERVOIR  INFORMATION 


Reservoir 

Location 

USGS 

Gaging  Station 

Drainaqe 
Area 
(sq.mi . ) 

Osable 

Storage 

(BP) 

'Cannonsville 

(modified) 

West  Branch 

Delaware  River 
at  Stilesville, 

N.Y. 

01425000 

456 

108.8 

Prompton 

Lackawaxen  River 
at  Prompton,  PA 

01429000 

59.7 

10.1 

Merrill  Creek 

Off-Channel  Delaware 
River  near  Belvidere 

none 

none 

15.0 

FE  Walter 

Lehigh  River  near 
White  Haven,  PA 

01447800 

290 

22.8 

Beltzville 

Pohopoco  Creek 
near  Parryville,PA 
draining  into  Lehigh 
River 

01449800 

109 

13.0 

Hackettstown 

Musconetcong  River 
at  Hackettstown,  N.d. 

01456000 

70.0 

9.9 

Nockamixon 

Tohickon  Creek  at 
Pipersville,  PA 
draining  into  the 
Delaware  River 

01459500 

97.4 

13.0 

1,2Blue  Marsh 

Tulpehocken  Creek 
at  Blue  Marsh 

Damsite,  PA 

01470960 

175 

7.5  3 

(4/1-9/30) 

(10/1-3/31) 

*Not  used  to  maintain  specific  Trenton  flow 
2Used  only  in  the  Total  Basin  Model 
includes  dead  storage 


TABLE  1-2 

RESERVOIR  COMBINATIONS 


COMBINATION 

RESERVOIRS 

COMBINATION 

RESERVOIRS 

1 

Beltzville  (13. 0)1 

Cannonsville  existing 

11 

Beltzville  (26.0) 
Cannonsville  modified 
Merrill  Creek 

2 

Beltzville  (13.0) 

Cannonsville  modified 

12 

Nockamixon 

Beltzville  (36.1) 

3 

Beltzville  (35.8) 

Cannonsville  modified 

FE  Walter  modified 

Cannonsville  modified 
Prompton  modified2 
Nockamixon 

4 

Beltzville  (13.0) 
Cannonsville  modified 
Merrill  Creek 

13 

Beltzville  (45.9) 

Cannonsville  modified 
FE  Walter  modified 
Merrill  Creek 

5 

Beltzville  (23.1) 

Cannonsville  modified 

Prompton  modified7 

Prompton  modified2 

14 

Beltzville  (48.8) 

Cannonsville  modified 

6 

Beltzville  (26.0) 
Cannonsville  modified 
Nockamixon 

FE  Walter  modified 
Merrill  Creek 
Nockamixon 

7 

Beltzville  (35.8) 

Cannonsville  modified 

FE  Walter  modified 

Merrill  Creek 

15 

Beltzville  (58.9) 

Cannonsville  modified 
FE  Walter  modified 
Merrill  Creek 

Prompton  modified2 

8 

Beltzville  (45.9) 

Cannonsville  modified 

Nockamixon 

FE  Walter  modified 
Prompton  modified2 

16 

Beltzville  (58.7) 
Cannonsville  modified 
FE  Walter  modified 

9 

Beltzville  (48.8) 

Cannonsville  modified 

FE  Walter  modified 
Nockamixon 

17 

Merrill  Creek 

Nockami xon 
Hackettstown 

Beltzville  (68.8) 

10 

Beltzville  (23.1) 
Cannonsville  modified 
Merrill  Creek 

Prompton  modified2 

Cannonsville  modified 
FE  Walter  modified 
Merrill  Creek 

Prompton  modified2 

Nockamixon 

Hackettstown 

Storage  contributing  to  maintain  Trenton  flow  in  BG  does  not 
include  Cannonsville  or  Merrill  Creek 

2Uses  a  routing  scheme  around  Montague 
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The  Phase  II  analysis  is  separated  into  three  major  categories: 
base  run,  reservoir  combinations  and  total  basin  modeling.  Each  preceding 
category  is  an  essential  step  to  the  next  category. 

The  first  step  consists  of  a  base  run  which  sets  the  general  basin 
conditions  for  the  reservoir  combinations  analysis.  This  base  run  emphasizes 
the  differences  be*"  :en  Phase  I  and  Phase  II.  It  encompasses  future  condi¬ 
tions  such  as  consumptive  uses  based  on  a  year  2000  projection,  plus  the 
augmented  conservation  releases  from  the  New  York  City  reservoirs  of  Pepacton, 
Cannonsville,  and  Neversink.  The  Delaware  and  Raritan  (D&R)  Canal  is  specifi¬ 
cally  operated  rather  than  using  historic  diversions.  Phase  II  uses  a  diver¬ 
sion  scheme  for  the  D&R  Canal  based  on  the  same  system  conditions  as  used  for 
the  New  York  City  reservoirs. 

The  second  step  is  the  use  of  the  Lower  Basin  Model  for  the  reservoir 
analysis  as  described  in  Chapter  III.  Unlike  the  Total  Basin  Model,  this 
model  has  a  reduced  study  area.  Montague's  50-year  flow,  already  prescribed 
by  the  previous  Base  Run,  used  in  order  to  keep  the  New  York  City  reservoirs' 
operations  constant,  creates  the  "headwaters"  to  the  Lower  Basin  Model  at 
Montague.  This  also  eliminates  the  backup  procedure  used  to  operate  the 
New  York  City  reservoirs.  The  study  area  is  further  reduced  by  running  the 
model  only  as  far  as  Trenton,  the  node  whose  flow  governs  the  operation  of 
the  lower  basin  reservoirs.  Figure  1-1  presents  pictorial ly  the  difference 
of  the  basin  study  area  between  the  two  models.  The  Lower  Basin  Model  con¬ 
siders  the  middle  portion  of  the  basin  with  Montague  acting  as  its  headwaters 
and  Trenton  as  the  end  of  the  system.  The  Total  Basin  Model  uses  the  basin 
in  each  simulation  from  the  headwaters  to  the  Delaware  Memorial  Bridge  at 
Wilmington,  Delaware. 

Several  of  the  combinations  treat  Prompton  as  a  source  for  meeting  the 
flow  objectives  at  Trenton.  Because  Prompton  is  upstream  of  Montague,  a  spe¬ 
cial  consideration  of  the  Montague  "headwaters"  is  made.  Briefly,  the  Prompton 
inflows  are  subtracted  from  the  Montague  flow  after  the  New  York  City  reser¬ 
voirs  have  completed  their  operations.  Then  the  Lower  Basin  Model  treats  the 
Prompton  Reservoir  as  another  node  flowing  into  the  Delaware  River  at 
Montague.  For  more  detail,  see  the  discussion  in  Chapter  III. 
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The  Total  Basin  Modelling,  or  third  step,  ties  the  reservoir  analy¬ 
sis  done  with  the  Lower  Basin  Model  to  the  total  basin  study  area  (Chapter  IV).  The 
worst  event's  maximum  sustainable  Trenton  objective  is  used  in  all  of  the 
total  basin  simulations.  Using  the  Total  Basin  Model,  two  separate  model¬ 
ling  analyses  are  performed. 

The  first  total  basin  modelling  analysis  deals  directly  with  the 
Prompton  Reservoir  releases  and  its  affects  on  both  the  New  York  City 
reservoir  releases  and  the  lower  basin  reservoir  releases.  For  this  analy¬ 
sis,  Beltzville  and  FE  Walter  are  on  line  to  augment  the  Trenton  flow  as 
well  as  Prompton.  This  is  Combination  8  from  Table  1-1.  Three  cases  are 
tested  in  this  category.  The  first  alternative  keeps  the  New  York  City 
reservoirs  releasing  the  same  amount  as  the  Base  Run  by  not  allowing  the 
Prompton  release  to  augment  the  Montague  flow.  This  is  done  by  changing 
the  Montague  target  whenever  the  Prompton  outflow  is  different  from  the 
natural  inflow.  This  is  referred  to  as  "Prompton  Around  Montague"  and  is  dis¬ 
cussed  more  fully  in  Chapter  IV.  The  second  alternative  allows  the  Prompton 
releases  to  help  augment  the  Montague  flow  while  being  operated  by  the 
Trenton  flow.  This  is  called  "Prompton  Through  Montague."  The  New  York  City 
reservoirs'  operating  rules  are  identical  for  each  of  these  two  simulations. 

The  third  and  last  alternative  has  a  different  set  of  operating  rules  for 
the  New  York  City  reservoirs  to  compare  with  the  second  alternative.  In 
all  three  cases,  the  Trenton  flow  objective  is  held  constant  to  the  worst 
event  maximum  objective  obtained  from  the  lower  basin  reservoir  analysis. 

The  second  analysis  compares  the  results  from  the  total  basin  runs 
for  Beltzville  only  (Combination  One)  and  all  reservoirs  on  line  (Combina¬ 
tion  17)  to  the  results  from  the  Phase  II  Base  Run.  Duration  and  low  flow 
frequency  analyses  for  these  three  cases  are  presented  in  Chapter  IV. 
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II.  BASE  RUN  PHASE  II 


The  base  run  conditions  for  Phase  II  differ  from  the  Phase  I 
base  run.  The  differences  are  due  to  changes  in  the  Phase  II  base 
run  to  reflect  the  addition  of  augmented  conservation  releases  re¬ 
quired  by  the  Delaware  River  Basin  Commission  for  Pepacton,  Cannonsville 
and  Neversink;  projected  consumptive  uses  for  the  year  2000;  and  speci¬ 
fied  Delaware  and  Raritan  (D  &  R)  Canal  diversions. 


The  New  York  City  reservoir  operating  rules  remained  the  same  in 
Phase  II  as  Phase  I.  Table  1 1 - 1  gives  the  operation  matrices  used  for  both 
Phase  I  and  II.  The  Montague  objectives  and  New  York  City  diversions  are 
identical.  The  conservation  releases  are  the  same  for  each  reservoir  for 
the  system  conditions  of  Drought  Warning  and  Drought.  However,  in  the 
Normal  Condition,  augmented  releases  are  made  from  each  reservoir  in  Phase 
II  while  in  Phase  I  it  is  again  the  basic  conservation  release.  In  Phase 
II,  future  projections  of  the  D&R  Canal  diversions  are  made  whereas  in 
Phase  I,  the  historical  diversions  were  used.  The  consumptive  use  increases 
in  the  Basin  from  1975  to  2000  are  incorporated  in  Phase  II.  Phase  I  has  the 
consumptive  uses  already  accounted  for  in  the  incremental  inflow. 

CONSERVATION  RELEASES 

The  basic  conservation  releases  used  in  Phase  I  have  been  updated 
to  include  augmented  releases  during  Normal  conditions.  These  additional 
releases  help  to  enhance  the  flow  conditions  of  the  rivers  downstream  of 
the  reservoirs.  They  are  required  by  the  New  Y  irk  State  Department  of 
Environmental  Conservation  and  have  been  agreed  upon  by  the  Delaware  River 
Basin  Commission  (DRBC).  Table  1 1 - 1  presents  the  difference  between  the 
Phase  I  and  II  conservation  release  schedule.  Table  II-2  gives  the  re¬ 
leases  specified  by  the  DRBC  Docket  No.  D- 77-20  for  the  New  York  City 
reservoirs  .  TI  , 


TABLE  II-l 


OPERATIONS  MATRICES 


PHASE 

I 

Normal 

Drought 

Montague  (cfs) 

1750 

1750 

1525 

NYC  Diversion 
(mgd) 

800 

600 

430 

D&R  Canal  (cfs) 

Actual  historic  record 

Conservation  Release* 

Basic 

Basic 

Basic 

Consumptive 

Uses 

Actual  historic  accounting 

PHASE 

II 

Normal 

Drought 

Montague  (cfs) 

1750 

1750 

1525 

NYC  Diversion 
(mgd) 

800 

600 

430 

D&R  Canal  (cfs) 

155 

124 

93 

Conservation  Release* 

Augmented 

Basic 

Basic 

Consumptive 

Uses 

Incremental 

increase  2000-1975 

-  See  Table  1 1 -4 

*See  Table  1 1 -3 

Source:  Task  Group  Report  DRBC  Docket  No.  D-77-20,  Appraisal  of  Upper 
Basin  Reservoir  Systems,  Drought  Fnorar>ncv  Criteria  and  Conservation  Measures, 
Delaware  River  Basin  Commission,  March  1979. 
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TABLE  I I -2 

NYC  RESERVOIRS  RELEASE  SCHEDULE  AS  SPECIFIED  IN  DRBC 


DOCKET  D-77-20* 


Reservoir 

Operation  Dates 

Basic  Conservation 

Neversink 

4/8 

-  10/31 

15.5  cfs 

11/1 

-  4/7 

4.6 

Pepacton 

4/3 

-  10/31 

18.5 

11/1 

-  4/7 

6.2 

Cannonsville 

4/16 

-  10/31 

23.2 

12/1 

-  3/15 

7.7 

Augmented  Release 

Neversink 

4/1 

-  10/31 

45  cfs 

11/1 

-  3/31 

25 

Pepacton 

4/1 

-  10/31 

70 

11/1 

-  3/31 

50 

Cannons vi lie 

4/1 

-  6/14 

45 

6/15 

-  8/15 

325 

8/16 

-  10/31 

45 

11/1 

-  3/31 

33 
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Because  of  programming  constraints,  these  releases  are  input  into 
the  model  on  a  monthly  basis  as  seen  in  Table  I I -3  and  not  on  a  daily  basis 
(Table  1 1-2 )  which  would  be  required  to  simulate  the  exact  changes  which 
occur  within  a  month.  For  the  basic  conservation  releases,  the  monthly  in¬ 
put  uses  the  maximum  release  required  for  that  particular  month  for  Neversink 
and  Pepacton.  Cannonsville  uses  the  lower  release  to  offset  the  higher  re¬ 
leases  of  the  other  two  reservoirs.  For  the  augmented  releases,  Neversink 
and  Pepacton  already  follow  a  monthly  schedule  in  Table  1 1-2 .  Cannonsville  does 
not.  For  the  months  of  June  and  August,  the  larger  of  the  two  releases  for 
Cannonsville  are  made  for  the  entire  month.  These  two  simplifications  are 
analyzed  to  determine  the  extent  of  the  effects  of  modelling  the  monthly 
releases  rather  than  the  daily  releases. 

An  analysis  is  performed  outside  the  model  to  show  the  differences 
between  the  two  release  schedules,  daily  and  modified  monthly.  For  the 
Drought  condition  releases,  the  month  in  question  is  April.  The  other  months 
are  identical  to  the  DRBC  specified  release  schedules.  The  model  releases 
41.8  cfs  from  all  three  reservoirs.  From  Table  1 1-2 ,  the  week  of  April  1-7 
should  release  18.5  cfs.  This  total  difference  is  23.3  cfs.  For  the  week 
of  April  8-15,  the  DRBC  Docket  states  that  the  total  should  be  41.8  cfs  and 
for  the  last  two  weeks,  the  total  release  should  be  57.3  cfs,  15.5  cfs  more 
than  in  the  model.  For  the  entire  month  of  April,  a  total  of  69.4  cfs-days 
will  not  be  released  in  the  model,  or  an  average  of  2.3  cfs  per  day.  By 
the  time  this  gets  routed  down  to  Montague,  the  difference  is  negligible. 

There  is,  however,  a  rather  large  flow  difference  for  the  augmented 
Cannonsville  release:  325  cfs  in  the  model  versus  45  cfs  in  the  Docket  at 
the  beginning  of  June  and  the  end  of  August.  This  increased  release  will 
drop  Cannonsville  into  the  Drought  Warning  Condition  from  the  Normal 
Condition  two  weeks  faster  than  normally  expected.  This  would  then  reduce 
the  conservation  release  from  the  augmented  release  of  325  cfs  to  the  basic 
release  of  only  23.2  cfs  instead  of  releasing  45  cfs  for  the  next  two  weeks. 

It  should  be  noted  that  these  releases  only  affect  the  immediate  downstream 
reach.  However,  the  operation  of  Cannonsville  in  conjunction  with  the 
other  two  reservoirs  adequately  models  the  flows  at  Montague. 
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Source:  Recommendation  of  New  York  State  Department  of  Environmental  Conservation 


CONSUMPTIVE  USES 


Because  Phase  II  is  dealing  with  proposed  conservation  releases, 
modification  to  existing  reservoirs  and  the  construction  of  new  reservoirs, 
it  is  necessary  to  include  all  other  future  considerations,  particularly 
consumptive  uses.  Through  the  technique  of  subtracting  upstream  flow 
from  downstream  flow  to  achieve  the  incremental  inflows,  the  historical 
consumptive  uses  have  already  been  accounted  for.  Therefore,  only  the  in¬ 
creased  consumptive  use  is  considered.  The  year  2000  is  chosen  for  the 
baseline  future  consumptive  uses  and  1975  as  the  end  of  the  historical 
accounting.  The  appropriate  consumptive  use  increases  are  supplied  by  the 
Delaware  River  Basin  Commission  (DRBC)  at  specified  nodes  in  the  model  in 
cfs.  Table  I 1-4  shows  the  consumptive  uses,  broken  down  by  node  and  whether 
it  is  a  power  plant  consumptive  use  or  some  other  type  of  consumptive  use. 

The  power  plant  separation  was  done  in  order  to  isolate  the  contribu¬ 
tion  of  Merrill  Creek  Reservoir  to  the  Delaware's  flow.  The  design  of 
Merrill  Creek  is  such  that  when  it  is  in  operation,  specific  power  plants' 
uses  will  be  replaced  by  the  Merrill  Creek  Reservoir  releases.  These  power 
plants  are  scattered  throughout  the  basin  and  are  listed  in  Table  1 1 -4 .  A 
more  detailed  discussion  is  presented  in  Chapter  III. 

D&R  CANAL  DIVERSIONS 

The  final  change  from  Phase  I  to  Phase  II  is  the  operation  of  the 
Delaware  und  Raritan  Canal.  In  the  Phase  I  modelling  the  historical  diver¬ 
sions  are  used.  Phase  II  uses  the  proposed  diversions  of  the  canal.  These 
diversions  are  specified  by  DRBC.  The  proposed  diversions  are  a  function 
of  the  system  conditions  based  on  the  New  York  City  reservoir  storage  levels. 
When  the  reservoirs  are  in  a  Normal  condition,  the  diversion  from  the 
Delaware  River  into  the  canal  is  155  cfs.  The  diversions  for  the  Drought 
Warning  and  Drought  conditions  are  124  cfs  and  93  cfs  respectively  as  pre¬ 
sented  in  Table  I I - 1 . 
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TABLE  1 1-4 

CONSUMPTIVE  USES1 
NET  INCREASE 
(CFS) 
1975-2000 


Node 

Power 

Plant 

All 

Other 

Total 

Port  Jervis 

— 

1.95 

1.95 

Riegelsville 

9.022 

62.56 

71.58 

Pt.  Pleasant 

70. 53 

— 

70.50 

Trenton 

— 

30.39 

30.39 

Chester 

26. 94 

— 

26.905 

Mouth  of  Schuylkill 

— 

187.3 

187. 30 5 

Delaware 

Memorial  Bridge 

— 

71.98 

71.98s 

TOTAL 

106.42 

354.18 

460.60 

Delaware  River  Basin  Comriss  ion's  Level  B  Study 

2Electric  Generating  Station  Martins  Creek  #3  and  #4 

"  3Electric  Generating  Stations  Gilbert  #8  and  #9  and 

Limerick  #1  and  #2 

4Electric  Generating  Stations  Eddystone  #3  and  #4  and 
**  Chester  #10  and  #11 

‘  5A1 ready  adjusted  for  replacement  factor  for  estuarine  use 

4 

1 

1 
1 
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III.  LOWER  BASIN  MODEL-RESERVOIR  COMBINATION  ANALYSIS 


The  Lower  Basin  Model  is  a  condensed  version  of  the  Total  Basin 
Model.  Beginning  with  the  flow  at  Montague  as  the  headwaters,  the  study 
area  extends  down  to  Trenton,  as  seen  in  Figure  1-1.  It  is  here  at  Trenton 
that  flow  objectives  are  tested  to  find  the  maximum  sustainable  low  flow 
objective  for  various  reservoir  combinations.  Using  the  17  different 
reservoir  combinations  described  in  Table  1-1,  the  objectives  are  analyzed 
for  nine  historical  events,  ranging  from  the  worst  event  in  the  past  50 
years  to  an  average  event. 

LOWER  BASIN  MODEL 

Fixed  Montague 

By  limiting  the  study  area  to  begin  with  the  flow  at  Montague, 
the  back-up  procedure  used  to  operate  the  New  York  City  reservoirs  is  elim¬ 
inated  from  the  Lower  Basin  Model.  Primarily,  this  fixed  Montague  serves 
the  purpose  of  maintaining  constant  New  York  City  reservoir  operations  through¬ 
out  the  lower  basin  reservoir  analysis. 


The  fixed  Montague  flow  is  determined  from  the  Base  Run  with  no 
lower  basin  reservoirs  on  line.  The  flow  at  Montague  is  then  written  on  a 
separate  tape  to  be  used  as  input  to  the  Lower  Basin  Model.  Because  the 
proposed  diversions  of  the  D&R  Canal  are  based  on  the  New  York  City  reservoir 
storage,  they  are  also  written  on  this  separate  tape.  This  new  tape  is 
used  along  side  the  oriainal  incremental  inflow  tape  to  create  the  driving 
force  of  the  Lower  Basin  Model. 


III-l 


From  Table  1-1,  there  are  two  different  New  York  City  reservoir  input 
conditions:  1)  Cannonsville  existing  storage  at  95,7000  mg  and  2)  Cannonsville 
modified  storage  at  108,750  mg.  This  therefore  requires  two  separate  tapes 
to  be  generated  as  input  into  the  Lower  Basin  Model,  each  being  used  for  the 
appropriate  reservoir  combination. 

Prompton  Reservoir  Scheme 

All  the  reservoirs  analyzed  with  the  Lower  Basin  Model  are  located 
between  Montague  and  Trenton,  save  one.  The  one  important  exception  is 
Prompton.  It  is  located  in  the  upper  basin,  above  Montague.  Therefore,  a 
fixed  flow  condition  at  Montague  cannot  accurately  describe  the  upper  basin 
when  Prompton  is  on  line.  A  fixed  Montague  would  include  the  natural  flow 
from  Prompton,  not  as  it  is  operated  as  a  reservoir.  For  this  case,  another 
Montague  scheme  is  used  because  the  lower  basin  runs  are  designed  to  keep  the 
New  York  City  reservoirs'  releases  and  diversions  the  same  as  the  Base  Run 
or  the  increased  Cannonsville  storage  run. 

In  order  to  maintain  a  constant  operation  of  the  New  York  City 
reservoirs  regardless  of  whether  Prompton  is  on  line,  the  following  scheme 
is  devised  to  allow  the  New  York  City  reservoirs  to  receive  "credit"  for 
the  natural  inflow  at  Prompton  and  to  allow  Pennsylvania  to  use  Prompton 's 
additional  releases  to  maintain  Trenton's  flow.  Figure  III-l  shows  schema¬ 
tically  this  procedure.  The  previous  total  basin  runs  (with  Cannonsville 
existing  and  modified)  routed  Prompton' s  natural  inflow  down  to  Montague 
with  a  real  lag  time  of  about  18  hours.  This  then  fixes  the  New  York  City  re¬ 
servoir's  operations.  An  estimate  of  the  Montague  flow  without  the  Prompton 
natural  inflow  is  made.  This  is  done  so  the  lower  basin  can  operate 
Prompton  as  a  reservoir  and  use  its  natural  inflow  as  input  to  its  storage. 

When  necessary,  releases  can  then  be  made  to  maintain  the  specified  Trenton 
objective.  This  estimate  took  the  actual  flow  at  Montague  and  subtracted 
out  the  previous  day's  natural  inflow  to  Prompton,  a  safe  estimate  of 
Prompton' s  contribution.  Then  the  Lower  Basin  Model  adds  together  this 
ficticious  Montague  flow  and  Prompton's  reservoir  operations  with  Prompton’s 
releases  (including  spills,  conservation  and  Trenton  releases)  lagged  one 
day.  This  one  day  lag  is  done  to  be  consistent  with  the  subtraction  of  the 
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Montague  Real  with  Prompton  Natural 
(NYC  Reservoirs  operation  "keyed"  to 
this  node) 


I  Subtract  Yesterday's  Prompton  Inflow 
1 


FIGURE  III-l 

Prompton  Reservoir  Schematic 
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previous  day's  natural  Prompton  inflow  from  the  Montague  flow.  With  this 
scheme,  the  New  York  City  reservoirs'  operations  are  constant  throughout 
the  entire  reservoir  analysis  because  they  were  operated  on  Montague's 
flows  accounting  for  Prompton 's  natural  flow  before  it  was  subtracted  out. 

Trenton  Release  Scheme 

Unlike  the  Montague  release  scheme  which  backed  the  model  up  three 
days,  recalibrates  and  reroutes  the  New  York  City  releases,  the  Trenton 
scheme  is  far  simpler.  A  next-day  release  is  used  such  that  yesterday's 
Trenton  flow  operates  today's  releases  from  the  lower  basin  reservoirs. 

This  will  naturally  lead  to  some  "missed  days".  A  missed  day  is  one  on 
which  the  Trenton  flow  is  lower  than  the  objective.  With  no  alteration 
to  this  next  day  release  scheme,  a  severe  oscillation  is  set  up  due  to 
the  "miss  one  day  -  release  next  day"  type  of  schedule.  This  is  compounded 
by  the  fact  that  most  of  the  releases  are  not  seen  until  three  days  after 
the  point  of  low  flow.  This  three-day  lag  occurs  because  the  routing  time 
can  be  as  much  as  two  days  from  some  of  the  reservoi rs .  This  oscillation  is 
more  apparent  as  the  objective  goes  higher,  causing  the  actual  monthly 
average  flow  at  Trenton  to  be  as  much  as  150  cfs  lower  than  the  target. 

To  overcome  this  missed  target,  a  strategy  is  set  up  such  that  today's 
release  is  based  on  the  release  made  three  days  ago  by  subtracting  (that  re¬ 
lease  of  three  days  ago)  from  today's  Trenton  flow  to  better  estimate  the 
natural  flow  at  Trenton.  This  strategy  takes  into  consideration  both  the 
next  day  release  scheme  plus  the  two-day  maximum  route  time.  With  this 
method,  there  will  always  be  that  first  "missed  day".  However,  any  oscilla¬ 
tion  will  be  much  less  noticeable  and  monthly  averaged  missed  targets  are 
rare.  When  they  do  occur,  they  are  on  the  order  of  5-10  cfs  lower  than  the 
target,  an  acceptable  deficit. 
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Merrill  Creek 


The  Merrill  Creek  reservoir  is  used  specifically  to  replace  the 
water  to  be  used  by  the  power  plants  within  the  Delaware  River  Basin  (see 
Table  II-4).  That  consumptive  use  is  126.9  cfs.  This  includes  a  20.48 
cfs  power  plant  use  below  the  Delaware  Memorial  Bridge.  Merrill  Creek  is 
required  to  replace  these  consumptive  uses  when  the  Trenton  flow  is  below 
3000  cfs,  releasing  into  the  Delaware  River  just  above  Riegelsville.  The 
water  for  its  storage  is  skimmed  at  a  maximum  rate  of  145  cfs  from  the 
Delaware  River  at  Belvidere  when  the  Trenton  flow  is  greater  than  3145  cfs. 
These  release/skim  targets  are  operated  using  the  next  day  scheme  as  is 
the  Trenton  objective.  Because  of  this  and  the  fact  that  not  all  of  the 
126.9  cfs  release  will  be  seen  immediately  at  Trenton  due  to  routing,  the 
release/skim  targets  are  changed  to  3100  and  3200  cfs  respectively.  In 
the  simulation,  Merrill  Creek  acts  as  a  reservoir  but  is  not  operated  to 
maintain  the  Trenton  objective.  It  merely  skims  and  releases,  based  on 
the  previous  day's  flow  at  Trenton.  Because  Merrill  Creek  is  not  used  to 
maintain  any  particular  flow  objective,  its  evaporation  and  conservation 
release  are  set  equal  to  zero. 

This  release/skim  scheme  is  not  applicable  when  the  target  approaches 
and  exceeds  the  3100  cfs  mark.  Consider  an  objective  of  3300  cfs:  all 
the  other  reservoirs  are  releasing  to  maintain  3300  cfs  while  Merrill  Creek 
is  skimming  to  replenish  its  storage.  Once  Merrill  Creek  is  full,  it  would 
never  lose  any  volume  through  evaporation  or  conservation  release  nor  would 
it  release  anything  because  Trenton's  flow  would  always  be  above  3100  cfs. 
Merrill  Creek  would  therefore  quit  skimming  and  would  no  longer  be  counter¬ 
acting  the  other  reservoirs'  releases.  However,  the  initial  reason  for 
including  Merrill  Creek  to  replace  power  plant  uses  would  not  be  accomplished. 
For  this  reason,  the  release/skim  scheme  is  altered  when  the  objective  exceeds 
3100  cfs.  For  those  targets,  Merrill  Creek  would  release  at  the  target  and 
skim  at  100  cfs  above  the  target.  While  this  scheme  does  not  use  the  full 
potential  of  Merrill  Creek's  storage  capacity,  it  would  at  least  contribute 
to  the  basins  operations  by  making  some  releases  rather  than  none  at  all. 


Evaporation  Calculations  and  Conservation  Releases 


Information  received  from  the  Philadelphia  District  Corps  of  Engineers 
allowed  updating  and  improvement  in  the  evaporation  calculations.  The  New 
York  City  reservoirs'  evaporation  values  are  not  changed  so  that  the  total 
basin  runs  of  Phase  II  can  be  compared  with  the  final  regulated  run  of 
Phase  I.  The  lower  basin  reservoirs  use  the  new  information.  Table  1 1 1 - 1 
gives  the  values  used  for  the  lower  basin  reservoirs.  Appendix  M1  assigns  a 
percentage  of  the  yearly  lake  evaporation  for  each  reservoir  to  the  months  of 
May  through  September.  The  average  monthly  rainfall  for  each  reservoir  was 
taken  from  plate  8  of  Appendix  M.  The  net  evaporation,  as  discussed  fully 
in  the  Phase  I  report,  is  calculated  as  the  evaporation  minus  one-half  the 
precipitation  on  a  monthly  basis.  Table  1 1 1-2  gives  the  evaporation  cal¬ 
culated  for  each  of  the  lower  basin  reservoirs  including  Blue  Marsh. 

Also  specified  in  Table  1 1 1 - 2  are  the  conservation  releases  for  each 
of  the  lower  basin  reservoirs.  The  conservation  releases  are  set  equal  to 
the  7-day  low  flow  which  has  a  recurrence  interval  of  10  years.  The  Q7-10 
flows  are  selected  from  the  frequency  analysis  of  the  natural  simulation 
performed  in  Phase  I. 

RESERVOIR  COMBINATIONS  ANALYSIS 

The  Lower  Basin  Model  as  described  above  'is  used  to  analyze  the  17 
combinations  in  Table  1-2.  The  purpose  is  to  define  the  maximum  maintainable 
flow  objective  at  Trenton  for  each  of  the  reservoir  combinations  for  selected 
low  flow  annual  events. 

The  model  defines  a  near-perfect  operation  of  the  reservoirs.  It 
should  be  emphasized  that  the  resultant  maintainable  flows  are  best  used  for 
the  relative  comparison  of  reservoir  combinations.  It  is  likely  that  the 
modeled  flow  objectives  are  higher  than  could  actually  be  obtained  by  the 
physical  operation  of  the  existing  and  proposed  reservoirs. 

Report  on  the  Comprehensive  Survey  of  the  Water  Resources  of  the  Delaware 
River  Basin,  Appendix  M. 
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TABLE  III-l 

LAKE  EVAPORATION  FOR  LOWER  BASIN  RESERVOIRS 


Reservoir 

Annual 

Lake  Evaporation 
(inches) 

Percent 

Occurring 

May  -  Sep 

F.E.  Walter 

30.0 

71 

Prompton 

30.0 

71 

Hackettstown 

31.5 

70 

Nockamixon 

33.5 

68 

Blue  Marsh 

33.5 

68 

Beltzville 

29.0 

66 

Source:  Report  on 

the  Comprehensive  Survey  of  the  Water 

Resources 

oT  the  Delaware  River  Basin,  Appendix M 

U.S.  Army  Engineer  District,  Philadelphia,  April, 

1960. 
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TABLE  1 1 1-2 

OVitP  BASIN  RESERVOIRS  EVAPORATION  AND  CONSERVATION  RELEASES 


Definition  and  Calculation  of  Events 

Nine  events  are  chosen  to  analyze  the  different  reservoir  combina¬ 
tions.  These  are  the  worst  event,  2nd  worst  annual  event,  3rd  worst 
annual  event,  4th,  5th,  10th,  15th,  20th  and  25th  worst  annual  events.  The 
model  simulation  uses  only  these  nine  selected  events,  not  a  full  50-year 
simulation.  Therefore,  these  events  are  isolated  from  the  total  50  years 
of  information.  In  order  to  choose  these  events,  test  runs  are  made  on 
combinations  One  and  Seventeen.  Both  of  these  combinations  are  run  for 
two  Trenton  objectives,  2500  and  3200  cfs.  An  analysis  of  the  drawdowns 
of  the  reservoirs  during  these  four  50-year  simulations  is  carried  out  to 
aid  in  the  prediction  of  the  worst  events. 

The  definition  of  an  event  is  the  period  of  time  between  two  full 
reservoir  states,  i.e.,  from  full  reservoir,  drawdown,  recovery,  to  full. 

The  New  York  City  Water  Year,  June  1  to  May  31  is  chosen  as  the  period  to 
describe  annual  events  because  the  lower  basin  reservoirs  fill  up  almost 
entirely  every  June  or  July.  There  is  one  exception:  the  drought  from 
1964-1966  is  considered  a  single  event  because  the  reservoirs  never  achieve 
a  full  storage  level  durinn  those  two  years. 

The  severity  of  each  event  is  determined  by  the  maximum  drawdown 
of  the  reservoir  system  during  the  New  York  City  Water  Year.  This  drawdown  is 
the  absolute  change  in  volume  of  the  reservoir  storage  level.  This  defini¬ 
tion  will  eliminate  any  possible  overlapping  of  drought  severity.  The  two 
targets  are  chosen  to  give  a  wide  range  of  responses  by  the  reservoirs. 

With  only  Beltzville  on  line  (Combination  One)  the  3200  cfs  objective  produces 
many  consecutive  years  of  negative  storage.  An  absolute  difference  is  still 
calculated.  The  2500  cfs  objective  produces  far  fewer  negative  years.  For 
the  case  of  all  the  reservoirs  on  line  (Combination  17)  and  a  2500  cfs  ob¬ 
jective  there  are  no  negative  years  at  all,  however  the  absolute  change  in 
volume  is  still  used  to  predict  the  worst  events.  The  following  is  an  example 
of  only  Beltzville  trying  to  maintain  a  large  target,  resulting  in  consecu¬ 
tive  years  of  negative  storage. 
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A  beginning  point  for  an  event  could  be  negative  storage,  let  us  say 
-5  BG.  If  the  storage  then  dropped  to  -20  BG,  the  severity  is  only  -15  BG, 
not  -20  BG.  This  method  will  also  credit  a  year  with  recovery.  For  instance, 
reverse  the  example  above.  Begin  at  -20  BG  and  end  at  -5  BG.  The  "severity" 
is  +15  BG,  thus,  not  even  accountable  in  the  worst  event  list.  Only  draw¬ 
downs  in  the  storage  level  are  considered  deficits. 

The  maximum  drawdown  is  calculated  on  every  year  of  each  of  the  four 
runs  described  above.  Then,  each  year  within  each  run  is  ranked,  with  the 
highest  deficit  being  the  worst.  These  four  rankings  are  combined  to  produce 
the  final  rankings,  from  which  the  nine  events  are  chosen.  Table  1 1 1 -3  gives 
the  final  results.  As  expected,  1964-1966  was  the  worst  event  in  the  50  year 
history.  The  25th  worst  annual  event,  June  1,  1935  to  May  31,  1936  is  an 
average  year. 

The  nine  events  are  then  run  back  to  back  in  the  simulation.  Because 
they  are  totally  exclusive  events,  the  reservoir  storages  are  artificially  set 
to  their  maximum  storage  on  June  1,  the  beginning  of  each  annual  event.  The 
worst  event  which  covers  the  period  from  June  1,  1964  to  May  31,  1966  is  of 
course,  allowed  to  run  naturally  past  June  1,  1965.  The  3rd,  4th,  and  5th 
worst  events  have  very  similar  inflow  characteristics.  Because  of  this, 
different  reservoir  combinations  bring  about  a  different  order  of  worst  events 
This  is  the  case  for  the  events  of  1931-32  and  1957-58.  For  the  most  part, 
1957-58  is  the  third  worst  event  and  1931-32  is  the  fourth.  However,  in  the 
combinations  where  FE  Walter  or  Nockamixon  are  not  used  as  reservoirs  (Combina 
tions  1,  2,  4  and  5),  1931-32  is  the  third  worst  event.  Therefore,  for  these 
four  combinations,  the  final  ranking  in  Table  1 1 1 -4  uses  1931-32  for  the  third 
worst  event  instead  of  1957-58  that  all  the  other  combinations  use.  Examina¬ 
tion  of  the  Natural  Flow  Simulation  of  Phase  I  shows  that  the  inflows  to  FE 
Walter  and  Nockamixon  are  much  higher  for  1931-32  than  for  1957-58.  This 
allows  for  more  available  storage  to  maintain  a  higher  Trenton  objective. 

Procedure  of  the  Reservoir  Combination  Analysis 

The  definition  of  the  maximum  sustainable  flow  objective  is  that 
objective  which  draws  the  reservoirs  exactly  to  their  zero  usable  storage 
level.  Daily  releases  needed  to  maintain  the  Trenton  flow  are  made 


TABLE  I I 1-3 

EVENT  DEFINITION 
June  1  -  May  31 


EVENT 


Worst 


Simulation 
Group  1 


Simulation 
Group  2 


*0nly  2-year  event 


III-ll 


1964-66* 

1930-31 


1957-58 


1931-32 

1966-67 

1962-63 

1936-37 

1968-69 

1935-36 


1 

entirely  from  the  percent  fullest  reservoir  on  that  particular  day.  In 
order  to  exactly  predict  the  objective  for  a  particular  event,  an  inde¬ 
finite  number  of  objecitves  would  have  to  be  tested.  Therefore,  the 
number  of  objectives  is  limited  to  a  range  which  will  produce  positive  and 
negative  reservoir  storages.  From  these  objectives  and  storages,  a  zero 
storage  level  with  corresponding  objective  is  determined. 

The  nine  events  in  Table  1 1 1 -3  are  split  into  two  groups  for  simu¬ 
lation  purposes.  The  first  group  contains  the  first  five  events,  the  second 
group  contains  the  last  four  events.  Each  group  is  then  run  separately, 
testing  various  objectives  for  different  conditions.  This  is  done  because 
the  objectives  for  the  last  four  events  would  not  be  able  to  reliably  pre¬ 
dict  a  zero  storage  for  the  first  five  events  and  vice  versa.  The  first 
five  use  a  range  of  300-400  cfs  at  increments  of  100  cfs.  Therefore,  no 
more  than  four  objectives  are  tested  for  the  first  five  events.  The  last 
four  events  use  a  much  wider  range,  generally  on  the  order  of  900-1200  cfs, 
at  increments  of  300-400  cfs.  A  maximum  of  three  objectives  test  the  last 
four  events. 

The  methodology  of  estimating  the  flow  objectives  is  done  mathe¬ 
matically  and  graphically.  Given  two  runs  for  one  reservoir  combination 
where  one  produces  a  positive  storage  due  to  an  objective  too  small  and  the 
other  a  negative  storage  due  to  an  objective  too  large  for  a  particular 
event,  the  drawdown  periods  in  days  are  estimated  from  a  full  storage  level 
to  the  minimum  storage  level.  In  the  case  of  the  objective  being  too  small, 
the  remaining  storage  is  divided  by  this  number  of  days  and  an  approximate 
increase  in  the  target  for  each  day  is  found.  The  same  is  done  for  the  case 
of  the  negative  storage:  this  case  determines  the  decrease  in  the  target 
for  each  day.  These  two  estimates  are  averaged  and  become  the  mathematical 
average.  To  check  this,  the  specified  objective  is  plotted  aqainst  the 
resulting  storage,  one  negative  and  one  positive,  and  the  zero  deficit  point 
is  found  on  the  graph  by  connecting  the  two  points.  The  mathematical  average 
and  the  graphical  zero  point  are  used  to  determine  the  final  low  flow  ob¬ 
jective.  Using  both  of  these  methods,  the  separate  objective  estimates  are 
quite  close,  most  times  within  20  cfs.  The  final  results  are  given  in 
Table  1 1 1-4. 
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Merrill  Creek,  Pror'pton  (r),  Noeka^i *on, 

Hackettstown  ?970  3060  3170  3310  3380  3610  3920  4260  4840 


Results  of  the  Reservoir  Combination  Analysis 


The  results  show  that  as  more  reservoirs  are  put  on  line,  the 
maximum  sustainable  low  flow  objective  increases.  As  the  events  become 
less  severe,  the  objective  also  increases.  The  lowest  worst  event  target 
is  2270  cfs  for  Combination  One,  and  the  highest  is  2970  cfs  for  Combina¬ 
tion  17.  For  an  average  event  (25th  worst  event)  this  range  is  from  3330 
cfs  (Combination  One)  to  4840  cfs  (Combination  17).  There  are,  however, 
anomalies  to  the  general  trend  of  an  increased  target  with  increased 
available  storage. 

Prompton  has  less  storage  available  than  does  Nockamixon.  This 
would  lead  one  to  surmise  that  a  combination  which  replaces  Prompton  with 
Nockamixon  would  be  able  to  sustain  higher  maximum  flow  objectives.  The 
results  show  this  is  not  the  case  in  almost  half  of  similar  events.  Com¬ 
bination  5  (with  Prompton)  and  6  (with  Nockamixon)  are  qood  examples  of 
the  anomaly.  The  worst  event  has  a  higher  objective  using  Prompton  than 
Nockamixon,  2430  cfs  and  2400  cfs,  respecti vely . 

The  events  oscillate,  first  having  Prompton  better  for  flow  con¬ 
ditions,  then  Nockamixon.  The  years  where  Prompton  consistently  produces 
a  higher  objective  are  1964-1966  (2-year  worst  event),  1931-1932  (3rd 
worst  event)  and  1936-1937  (15th  worst  event).  On  a  closer  look  at  the 
incremental  inflows,  on  an  average  monthly  basis  Prompton  indeed  has 
much  more  flow  than  Nockamixon  for  these  years  in  question.  The  combina¬ 
tion  pairs  of  8  and  9,  10  and  11,  13  and  14  (with  Prompton  and  with 
Nockamixon  respectively),  all  show  similar  patterns  as  5  and  6.  The  basin 
obviously  reacts  differently  under  stress  conditions.  For  an  average  year 
(25th  worst  event),  Prompton  does  not  come  close  to  maintaining  the  objective 
that  Nockamixon  can.  This  is  just  the  case  for  Combinations  5  (with  Prompton) 
and  6  (with  Nockamixon)  where  they  maintain  an  objective  of  3650  and  3730 
cfs  respectively. 
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A  second  slight  irregularity  in  Table  III -4  is  Combinations  1  and  2, 
Cannonsville  existing  and  modified  storage,  respectively.  With  Cannons vi lie 
releasing  to  meet  the  Montague  objective,  there  is  very  little  difference  in 
the  maintainable  Trenton  objective,  particularly  in  the  first  five  events. 
The  same  is  true  for  the  last  four  events  with  only  an  exchange  of  ranks. 
Again,  the  difference  is  very  slight  and  can  be  considered  almost  identical 
objectives. 
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IV.  TOTAL  BASIN  MODELING 


The  purpose  of  the  Lower  Basin  Model  is  to  predict  a  flow  objective 
for  a  particular  reservoir  combination  for  a  given  event  as  described  in 
Chapter  III.  It  is  to  do  this  calculation  as  efficiently  as  possible, 
looking  at  a  very  reduced  study  area.  Once  the  maximum  sustainable  low 
flow  objective  has  been  determined  for  the  worst  event,  the  next  step  is 
to  model  and  analyze  the  effects  of  using  that  objective  with  the  Total 
Basin  Model  for  various  New  York  City  reservoir  operating  schemes  and 
different  techniques  used  to  credit  Prompton’s  releases  to  Montague  or 
Trenton.  The  statistical  flow  characteristics  will  also  be  analyzed  for 
the  Base  Run,  Combination  One  (Beltzville  only)  and  Combination  Seventeen 
(all  reservoirs  operating). 

For  the  total  basin  modeling  analysis  of  the  Prompton  and  New  York 
City  reservoir  operating  schemes,  a  single  reservoir  Combination  from  the 
lower  basin  analysis  has  been  chosen.  This  is  Combination  Eight  in  Table 
1-2  which  includes  Beltzville,  Cannonsville  modified,  FE  Walter  and 
Prompton  Reservoirs.  This  combination  is  suitable  for  examining  the  effects 
of  Prompton's  different  crediting  schemes  and  for  the  analysis  of  the 
special  set  of  New  York  City  operating  rules.  The  Trenton  flow  objective 
is  the  result  of  the  reservoir  analysis  in  Chapter  III.  It  is  the  worst 
event's  maximum  sustainable  flow  for  the  specified  reservoir  combination 
of  2770  cfs.  This  will  be  held  constant  for  this  portion  of  the  Total 
Basin  Modeling. 

Because  the  study  area  with  the  Total  Basin  Model  encompasses  the 
basin  down  to  the  Delaware  Memorial  Bridge,  another  basin  operation  is 
included.  Blue  Marsh  Reservoir  is  located  in  the  Schuylkill  River  Basin. 
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The  Schuylkill  River  makes  its  confluence  with  the  Delaware  River  between 
Trenton  and  the  Delaware  Memorial  Bridge.  Blue  Marsh  does  not  release  to 
maintain  any  specific  flow  objective.  It  is  however,  modeled  as  a  reservoir 
with  a  conservation  release,  evaporation  factors  and  supports  a  water  supply 
withdrawal  for  the  Schuylkill  River  basin.  This  reservoir  operation  is 
essential  to  proper  modeling  of  the  basin  around  the  area  of  Tulpehocken 
Creek  and  Schuylkill  River.  Blue  Marsh  has  a  conservation  release  of  41  cfs 
and  a  water  supply  withdrawal  of  9  cfs.  In  the  summer  months  between 
May  1  and  September  30,  Blue  Marsh  is  used  as  a  recreational  area.  For  this 
reason,  its  maximum  storage  is  raised  to  22,900  acre-feet  from  17,600  acre- 
feet  which  it  maintains  for  the  rest  of  the  year. 

PROMPTON  RESERVOIR 

Prompton  Reservoir  is  located  in  the  upper  basin  above  Montague. 

Its  releases  are  naturally  routed  to  Montague  via  the  Lackawaxen  River. 
Because  it  is  releasing  to  maintain  a  Trenton  flow  objective  and  at  the 
same  time  must  pass  by  Montague,  its  releases  could  help  augment  either 
flow  objective.  Two  different  techniques  of  crediting  Prompton's  augmen¬ 
tation  are  analyzed.  The  first  is  identical  to  the  Lower  Basin  Model  in 
that  the  New  York  City  reservoirs  are  not  allowed  to  change  their  releases 
from  the  base  run.  This  run  is  referred  to  as  Prompton  Around.  Montague. 
Prompton's  releases  are  not  used  to  augment  Montague's  flow. 

The  second  run  does  just  the  opposite:  Prompton's  releases,  even 
though  they  are  determined  according  to  the  Trenton  flow  objective,  are 
used  to  augment  Montague's  flow  to  help  meet  its  objective.  This  is 
called  Prompton  Through  Montague  and  results  in  the  New  York  City  reser¬ 
voirs  operating  differently. 

Prompton  Around  Montague 

In  this  case,  Prompton's  reservoir  releases  are  modeled  to  maintain 
the  flow  objective  at  Trenton  without  crediting  the  releases  to  Montague's 


flow.  The  releases  are  actually  routed  to  Montague.  A  special  technique 
is  used  to  adjust  the  Montague  objective  to  take  into  account  the  additional 
flow  created  by  Prompton's  releases  and  to  keep  the  operations  of  the  New 
York  City  reservoirs  identical  to  the  Base  Run  made  with  Cannonsville 
modified  storage.  This  is  accomplished  in  the  Total  Basin  Model  by  using  a 
simple  mass  balance  to  raise  or  lower  the  Montague  objective: 

Change  in  Objective  =  (Prompton  Q710  Conservation  Releases  + 
Prompton  Spills  +  Prompton  Trenton  Releases)  -  Prompton 
Natural  Inflow 

Montague  Objective  =  Original  Objective  +  Change  in  Objective. 

The  change  in  objective  is  positive  (raising  the  Montague  objective) 
when  Prompton  reservoir  releases  are  more  than  its  natural  inflow.  This 
occurs  when  Trenton  is  in  a  low  flow  condition.  The  additional  Prompton 
releases  naturally  raise  the  Montague  flow.  If  the  objective  were  not 
raised  also,  the  New  York  City  reservoirs  would  release  less  than  in  the 
base  run.  The  change  in  objective  is  negative  (reducing  the  Montague  ob¬ 
jective)  when  the  natural  inflow  is  greater  than  a  release.  This  occurs 
usually  when  Prompton  is  recovering  from  a  drawdown  period,  thus  holding 
back  most  of  the  natural  inflow.  The  decrease  in  natural  inflow  from 
Prompton  produces  a  lower  flow  at  Montague  which  would  cause  the  New  York 
City  reservoirs  to  release  more  if  the  objective  were  not  also  reduced. 

Exhibits  I V- 1  through  IV-4  give  samples  of  the  output  for  Prompton 
Around  Montague.  Exhibits  I V- 1  and  IV-2  are  the  actual  monthly  average 
flows  at  Montague  and  Trenton  respectively.  Exhibits  IV-3  and  IV-4  are 
the  average  monthly  total  storages  for  the  New  York  City  reservoirs  and 
the  lower  basin  reservoirs,  respectively.  These  four  basin  results  are 
compared  to  the  Prompton  Through  Montague  simulation  and  are  discussed 
in  the  following  paragraphs. 
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Prompton  Through  Montague 


For  this  case,  Prompton 's  releases  are  governed  by  the  Trenton 
objective,  yet  allowed  to  help  augment  the  Montague  flow.  Therefore,  no 
changes  are  made  in  the  Montague  objective.  Using  the  same  Trenton  ob¬ 
jective  of  2770  cfs  as  Prompton  Around  Montague,  the  results  are  shown 
in  Exhibits  IV-5  through  IV-8.  Exhibits  IV-5  and  IV-6  are  the  monthly 
average  flows  at  Montague  and  Trenton  respectively.  Exhibits  IV-7  and 
IV-8  are  the  average  monthly  total  storages  for  the  New  York  City  reser¬ 
voirs  and  the  lower  basin  reservoirs,  respectively. 

Prompton 's  releases  which  in  this  case  are  now  credited  to  the 
Montague  flow  reduce  the  required  releases  from  the  New  York  City  reser¬ 
voirs.  Exhibit  IV-7,  which  is  the  total  New  York  City  reservoir  storage 
shows  a  surplus  of  12.1  bg  through  the  60 ' s  drought.  Exhibit  IV-3,  total 
New  York  City  reservoir  storage  for  the  Prompton  Around  Montague  case, 
shows  a  deficit  of  6.2  bg  through  the  60 ’ s  drought.  Therefore,  a  net  gain 
of  18.3  bg  is  seen  by  the  New  York  City  reservoirs,  when  Prompton1 s  releases 
help  augment  Montague's  flow. 


In  the  case  of  Prompton  Through  Montague,  the  flows  at  Montague 
seen  in  Exhibit  IV-5  will  be  less  than  the  other  case  of  PromDton  Around 
Montague  given  in  Exhibit  IV-1.  With  Prompton  credited  to  the  Montague 
flow,  the  objective  never  changes  and  therefore  will  result  in  lower  average 
flows.  The  periods  to  examine  are  August  -  November  of  1964  and  June  - 
September  of  1965.  Prompton  Around  Montague  produces  an  average  flow  for 
each  period  of  1760  cfs  and  1660  cfs  respectively.  Prompton  Through  Montague 
produces  1660  and  1550  cfs  for  each  respectively,  100  cfs  lower  than  Prompton 
.Around  Montague. 


This  reduction  of  the  Montague  flow  makes  less  water  available  to 
the  lower  basin.  With  the  Trenton  objective  constant  for  the  two  cases,  the 
lower  basin  reservoirs  must  therefore  release  more  to  maintain  the  same 
Trenton  objective.  The  flows  at  Trenton  for  each  case  are  nearly  the  same. 
The  Trenton  flows  for  those  same  periods  average  2790  cfs  and  2810  cfs 
respectively  for  the  Prompton  Around  Montague  case  as  seen  in  Exhibit  1 1/-2 
and  2810  cfs  and  2830  cfs  for  Prompton  Through  Montague  given  in  Exhibit 
IV-6. 
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The  storage  of  the  lower  basin  reservoirs  also  reflect  the  additional 
releases  required  to  maintain  the  Trenton  objective.  During  the  60 ' s  drought, 
when  Prompton's  releases  are  credited  to  Montague  in  the  Prompton  Throupl. 
Montague  case,  five  consecutive  months  of  negative  storage  occur  in  the  lower 
basin  reservoirs  as  seen  in  Exhibit  IV-8  with  a  maximum  deficit  of  13.7  bg. 
When  Prompton's  releases  are  not  credited  to  Montague  in  the  Prompton  Aj-c.eS 
Montague  case,  the  total  lower  basin  reservoir  storage  show  a  slight  surplus 
of  4.1  bg  as  seen  in  Exhibit  IV-4.  The  total  loss  bv  the  lower  basin  rpser- 
voirs  is  17.8  bg,  almost  exactly  equal  to  the  gain  by  the  New  York  City 
reservoirs . 

In  conclusion,  when  Prompton's  releases  are  credited  to  the  Montague 
flow,  the  New  York  City  reservoirs  benefit  by  making  less  releases  while 
the  Montague  flow  is  reduced  in  magnitude  and  the  lower  basin  reservoirs  re¬ 
lease  much  more  water  thus  reducing  their  available  storage. 

NEW  YORK  CITY  RESERVOIR  OPERATIONS 

For  the  Prompton  Time  Montague  case  described  above,  the  increase 
in  the  New  York  City  reservoir  storages  would  permit  greater  diversions.  A 
new  set  of  New  York  City  operating  rules  specified  by  the  New  York  State 
Department  of  Environmental  Conservation  are  analyzed.  In  this  section, 
the  results  from  Prompton  \  Montague  are  compared  to  the  results 
from  the  special  operating  rules. 

Table  IV-1  gives  both  the  special  operating  rules  and  those  used 
in  the  rest  of  the  reservoir  analysis,  referred  to  here  as  the  Cvini*ut~ 
operating  rules.  The  Normal  System  Condition  remains  the  same.  The 
Special  Drought  Warning  Condition  reduces  the  Montague  objective  by  100  cfs 
and  raises  the  diversions  by  50  mgd.  The  Special  Drought  Condition  raises 
the  diversions  by  120  mgd  and  keeps  the  Montague  objective  the  same.  Both 
of  these  cases  use  the  same  flow  objective  at  Trenton  of  2770  cfs  and  use 
the  Prompton  releases  to  credit  Montague's  flow.  Exhibits  IV-9  through 
I V- 1 2  of  the  Special  rules  are  compared  with  Exhibits  I V- 5  through  IV-8  of 
the  Prompton  Tim-:  ah  Montague  case.  Exhibits  IV-9  and  IV-10  are  the  average 
monthly  flows  at  Montague  and  Trenton  respectively.  Exhibits  IV-11  and 
I V- 12  are  the  average  total  storage  of  the  New  York  City  reservoirs  and  the 
lower  basin  reservoirs,  respectively. 
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TABLE  IV-1 


NEW  AND  EXISTING  NEW  YORK  CITY  RESERVOIR 
OPERATING  RULES 


CONDITION 


Operation 

Normal 

Drought 

Warning 

Drought 

Montague 

Objective 

Special* 

1750 

1650 

1525 

(cfs) 

Original 

1750 

1750 

1525 

New  York 

City 

Special* 

800 

650 

550 

Diversions 

(mgd) 

Original 

800 

600 

430 

^Specified  by  New 

York  State 

Department 

of  Environmental 

Conservation 

The  most  striking  comparison  is  the  total  New  York  storage  avail¬ 
able  during  the  60's  drought  shown  in  Exhibit  IV-11  for  the  Special  rules 
and  in  Exhibit  IV-7  for  the  Original  rules.  The  Original  rules  never 
create  a  deficit  throughout  the  60's  and  actually  end  up  with  a  surplus  of 
12.1  bg.  The  Special  rules  drawdown  the  reservoirs  such  that  six  months 
of  consecutive  deficits  occur  from  September  1965  to  February  1966  with  a 
maximum  monthly  average  deficit  of  -34.5  bg.  Seven  months  later  in 
October  1966  (water  year  1967),  a  four-month  deficit  begins  with  a  maximum 
monthly  average  deficit  of  -26.2  bg. 

Because  both  of  these  cases  use  the  Prompton  Through  Montague  con¬ 
dition  for  crediting  the  Prompton  releases,  there  is  not  too  much  difference 
in  the  average  monthly  flows  at  Montague.  Comparing  Exhibit  IV-9  from 
the  Special  rules  and  Exhibit  IV-5  from  the  Original  rules,  the  only  dif¬ 
ferences  occur  when  one  considers  that  the  Drought  Warning  Montague  objective 
is  reduced  from  1750  cfs  to  1650  cfs  in  the  Special  rules.  This  does, 
however,  affect  the  operations  of  the  lower  basin  reservoirs.  With  a 
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reduction  of  100  cfs  in  the  Montague  flow,  less  water  is  routed  down  to 
Trenton.  Thus,  the  lower  basin  reservoirs  must  release  more  to  maintain 
the  same  Trenton  objective  as  before.  This  is  only  compounded  by  the 
treatment  of  Prompton's  releases,  allowing  them  to  help  augment  the 
Montague  flow  and  therefore  are  not  directly  augmenting  the  Trenton  flow. 
Exhibit  IV-8  gives  the  total  storage  for  the  lower  basin  reservoirs  for 
the  fvipinal  rules.  Exhibit  IV-12  shows  the  same  for  the  Special  rules. 

During  the  30's  drought,  the  number  of  consecutive  average  monthly  negative 
storages  increases  with  the  Special  rules  as  well  as  the  severity  from 
-4.0  bg  to  -7.4  bg.  The  60's  drought  produces  similar  results.  The  maximum 
deficit  increases  from  -13.7  bg  using  the  Base  Run  rules  to  -18.1  bg  using 
the  special  rules  and  extends  the  consecutive  negative  months  by  one. 

The  change  in  the  New  York  City  reservoir  operating  rules  do  not  have  as 
strong  an  effect  on  the  lower  basin  reservoirs  as  for  the  New  York  City 
reservoirs.  Nonetheless,  any  operating  rule  change  for  the  upper  basin  will 
for  this  case  draw  down  the  lower  basin  reservoirs  even  further  by  requiring 
increased  releases. 

The  increased  releases  from  the  lower  basin  reservoirs  maintain 
approximately  the  same  Trenton  flow  with  either  operating  rules.  This  is 
shown  in  Exhibit  I V- 10  which  gives  the  average  monthly  flow  at  Trenton 
for  the  Special  rules  and  Exhibit  IV-6  gives  the  results  from  the  Original 
rules. 

LAKE  WALLENPAUPACK 

Lake  Wal lenpaupack  has  been  regulated  by  a  power  plant  since  1925. 

Its  regulated  releases  have  been  designated  as  natural  observed  outflow 
at  that  node  and  used  as  such  in  all  analyses  in  both  Phase  I  and  II. 

However,  because  of  its  regulation,  it  is  conceivable  that  Lake  Wal lenpaupack 
could  have  supplemented  the  flow  at  Montague  during  the  60' s  drought  even  more 
than  the  historical  contribution  already  made.  With  this  in  mind,  a  release 
scheme  for  this  study  only  was  devised  in  coordination  with  the  Pennsylvania 
Department  of  Environmental  Resources  such  that  maximum  releases  are  made 
without  considering  optimal  generation  of  hydroelectric  power.  Any  releases 
from  Wallenpaupack  will  help  augment  Montague's  flow  to  meet  the  1750  cfs 
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objective,  thus  reducing  the  releases  from  the  New  York  City  reservoirs. 
This  section  compares  the  results  of  using  these  special  Wallenpaupack 
releases  to  the  results  of  using  the  historical  outflows  for  two  reser¬ 
voir  combinations,  2  (Beltzville  and  Cannonsville  modified)  and  17  (all 
reservoirs  on  line). 

The  scheme  in  use  begins  in  June  1962  when  the  New  York  City 
reservoirs  first  fall  into  Drought  Warning  Condition  and  ends  in  May  1967 
when  the  New  York  City  reservoirs  are  again  in  the  Normal  Condition. 

Lake  Wallenpaupack's  maximum  releases  for  each  day  of  the  month  are  de¬ 
termined  from  a  minimum  elevation  schedule.  Lake  Wallenpaupack  is  normally 
full  at  elevation  1187  feet  (not  including  extra  flood  storage  capacity). 
The  minimum  elevation  of  the  reservoir  is  1165  feet.  A  strict  schedule 
is  set  up  from  June  1  to  December  1  of  minimum  allowable  elevations  as 
given  in  Table  IV-2.  The  releases  will  drawdown  Lake  Wallenpaupack  in 
order  to  arrive  at  these  storage  elevations. 

TABLE  IV-2 

MINIUM  ELEVATIONS  OF  LAKE  WALLENPAUPACK 


June  1 
July  1 
August  1 
September  1 
October  1 
November  1 
December  1 


1187  feet 
1184  feet 
1180  feet 
1177  feet 
1172.5  feet 
1168  feet 
1165  feet 
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The  maximum  releases  for  each  month  are  determined  by  first 
calculating  net  inflows  to  the  reservoir  from  actual  monthly  changes  in 
storage,  and  actual  monthly  releases.  Natural  inflows  to  Lake  Wallen- 
paupack  are  not  available  as  such  and  cannot  be  determined  from  a  corre¬ 
lation  and  extension  procedure;  the  USGS  station  at  Wal lenpaupack  has 
always  recorded  regulated  releases  from  the  water  power  reservoir.  From 
the  elevation  schedule,  the  change  in  storage  is  calculated.  Accounting 
for  the  net  inflow,  the  maximum  release  possible  is  determined.  For  ex¬ 
ample,  the  June  releases  will  draw  down  the  storage  elevation  by  three 
feet  to  obtain  a  July  1  elevation  of  1184  feet.  Adding  the  inflow  to 
this  amount,  the  maximum  releases  for  each  day  are  calculated.  The 
releases,  being  a  function  of  variable  inflow  change  from  month  to  month, 
and  each  month  from  year  to  year.  For  example,  although  the  month  of 
June  releases  are  set  as  high  as  possible  as  are  all  the  months,  each 
June  will  not  release  the  same  amount. 

From  December  1  to  June  1,  Lake  Wallenpaupack  refills  using  the 
calculated  net  inflows  such  that  the  maximum  possible  storage  is  avail¬ 
able  on  June  1.  This  means  that  although  a  conservation  release  of 
70  cfs  will  try  to  be  maintained,  it  is  not  always  possible  to  release 
70  cfs  in  order  to  have  a  full  reservoir  on  June  1. 

The  results  of  the  two  combinations  show  that  Lake  Wallenpaupack 
is  indeed  helpful  in  maintaining  the  Montague  objective.  Fewer  releases 
are  made  by  the  New  York  City  reservoirs  thus  allowing  more  water  to  be 
available  to  raise  the  minimum  storage  remaining  in  the  New  York  City 
reservoirs  at  the  critical  low  point  of  the  simulation.  Both  Combinations 
2  and  17,  when  comparing  the  special  releases  versus  historical  releases 
from  Lake  Wallenpaupack  result  in  an  increased  remaining  storage  capacity  of 
19  bg  for  the  New  York  City  reservoirs  when  the  special  releases  are  made. 
The  lower  basin  reservoirs  which  are  still  operating  to  maintain  a  set 
Trenton  objective,  are  hardly  affected.  Slight  increases  in  Montague's 
flow  require  less  releases  to  be  made  from  the  lower  basin  reservoirs,  pro¬ 
ducing  an  increase  in  the  remaining  total  storage  of  the  lower  basin  reser¬ 
voirs  of  approximately  1  bg  for  each  combination. 
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STATISTICAL  ANALYSIS 


As  part  of  the  total  basin  modeling  process,  two  additional  combina¬ 
tions  were  run  for  the  entire  50  years.  The  two  extreme  conditions  were 
chosen:  Combination  One  (Beltzville  only)  and  Combination  17  (all 
reservoirs  on  line) .  Statistics  were  generated  for  each  of  these  runs  as 
well  as  the  Base  Run.  Appendix  C  contains  the  results  from  the  Base  Run 
described  in  Chapter  II.  Appendix  0  and  E  contain  Combination  One  and  17, 
respectively.  For  each  of  the  44  key  locations,  duration  tables  and  curves 
are  produced.  Table  C-l,  D-l,  and  E-l  contain  these  tables  for  the  Base 
Run,  Combination  One  and  Combination  17,  respectively.  Figures  C-l  through 
C-44,  D-l  through  D-44  and  E-l  through  E-44  are  the  duration  curves  for 
the  respective  runs. 

For  each  of  the  44  key  locations  Log-Pearson  low  flow  frequency 
tables  are  produced  for  periods  of  1 ,  3,  7,  14,  30,  60,  90,  120,  183  and 
365  consecutive  days.  Tables  C-2.1  through  C-2.44  give  the  results  of  the 
Base  Run  in  Appendix  C.  Tables  D-2.1  through  D-2.44  show  the  results  of 
Combination  One.  Finally,  Tables  E-2.1  through  E-2.44  give  the  corresponding 
results  of  Combination  17.  Eighteen  stations  listed  in  Table  IV-3,  chosen 
from  the  44  key  locations,  have  low  flow  frequency  curves  for  the  7-day  and 
120-day  consecutive  day  portion  of  the  tables.  They  are  Figures  C-45 
through  C-62  for  the  Base  Run,  Figures  D-45  through  D-62  for  Combination  One 
and  Figures  E-45  through  E-62  for  Combination  17. 

Many  of  the  stations  that  are  downstream  from  a  reservoir  have 
irregular  statistics.  The  Log-Pearson  method  used  cannot  accurately  account 
for  basic  conservation  releases.  It  tries  to  fit  a  consistent  statistical 
line  through  the  data.  Because  of  this,  the  minimum  statistical  flow  is 
always  less  than  the  basic  conservation  release.  This  is  true  for  Neversink, 
Pepacton  and  Cannonsville  as  well  as  the  lower  basin  reservoirs.  The  low 
flow  frequency  tables  have  been  adjusted  to  show  the  basic  conservation  re¬ 
lease,  not  the  calculated  flow  value.  Also,  as  the  number  of  consecutive 
days  increases,  so  do  the  maximum  flow  values.  However,  directly  downstream 
of  a  reservoir,  the  minimum  flow  will  not  change  as  rapidly,  if  at  all. 


TABLE  I V- 3 


KEY  LOCATIONS  FOR  LOW  FLOW 
FREQUENCY  CURVES 


Location 


East  Branch  Delaware  River 
at  Downsville,  NY 

West  Branch  Delaware  River 
at  Stilesville,  NY 

West  Branch  Delaware  River 
at  Hale  Eddy,  NY 

Delaware  River  near 
Cal li coon,  NY 

Delaware  River  near 
Barryville,  NY 

Lackawaxen  River  at 
Hawley,  PA 

Delaware  River  at 
Port  Jervis,  NY 

Neversink  River  at 
Neversink,  NY 

Delaware  River  at 
Montague,  NJ 

Pohopoco  Creek  at 
Beltzville  Damsite,  PA. 

Lehigh  River  at 
Bethlehem,  PA 

Tohickon  Creek  at 
Pipersville,  PA 

Delaware  River  at 
Trenton,  NJ 

Tulpehocken  Creek  at 
Blue  Marsh  Damsite,  PA. 

Schuylkill  River  at 
Reading,  PA 

Schuylkill  River  at 
Philadelphia 

Delaware  River  below 
Mouth  of  Schuylkill 

Delaware  River  at 
Delaware  Memorial  Bridge 


USGS  Station 
01417000 

01425000 

01426500 

01427405 

01428500 

01431500 

01434000 

01436000 

01438500 

01449800 

01453000 

0145«500 

01463500 

01470960 

01741500 

01474500 

None 

None 


IV-23 


This  causes  the  frequency  curve  for,  say,  the  120  day  period  to  actually  fall 
below  or  cross  the  frequency  curve  for  the  90  day  period.  Whenever  these  in¬ 
consistencies  occur,  the  flows  are  plotted  against  duration  on  log-log  paper 
for  the  particular  probability.  Adjustments  are  made  so  that  low  flow  always 
increases  with  duration.  Therefore,  as  in  Phase  I,  these  tables  were  adjusted 
to  smooth  out  the  curves.  In  such  cases,  the  adjusted  values  have  been  de¬ 
signated  in  the  appendices  by  asterisks  and  called  "Recalculated". 

A  few  of  the  stations  at  the  smaller  streams'  headwaters  actually 
have  a  few  values  of  zero  flow.  Because  the  statistical  analysis  uses  the 
Log- Pearson  method,  these  zero  values  cannot  be  used  to  determine  the  pro¬ 
babilities  and  corresponding  flows.  Instead,  the  A969  program  uses  all  non¬ 
negative  flow  values  to  determine  the  flow  values  and  then  assigns  an  adjusted 
probability  which  assumes  there  were  50  values  used  instead  of  just  the  non¬ 
negative  flow  values.  To  arrive  at  the  required  probabilities,  the  flows  are 
plotted  versus  the  adjusted  probabilities.  The  required  probabilities  are 
then  picked  off  the  graph.  These  cases  are  designated  in  the  Appendices  as 
being  "Calculated  from  the  adjusted  probability". 

The  final  statistical  irregularity  involves  Riegelsville.  Because 
of  some  negative  inflow  values,  Riegelsville  has  flow  values  which  are  less 
than  its  upstream  node,  Belvidere.  The  Log-Pearson  statistical  method  calcu¬ 
lated  frequencies  that  do  not  consistently  increase  from  Belvidere  to  Riegelsvi  11 
to  Trenton.  The  frequencies  for  Riegelsville  have  therefore  been  adjusted 
based  on  a  drainage  area  ratio  and  the  flows  at  Belvidere  and  Trenton, whose 
frequencies  are  not  adjusted. 

The  differences  in  the  durations  and  low  flow  frequencies  between 
the  three  runs  arise  from  the  various  reservoirs  being  put  on  line  to  help 
supplement  Trenton's  flow.  The  Base  Run  is  identical  to  Combination  One 
except  where  Beltzville  is  introduced.  Comparing  the  duration  tables  from 
Appendix  C  and  D,  Base  Run  and  Combination  One  respectively,  Pohopoco  Creek 
below  the  Beltzville  damsite  shows  that  the  low  flows  have  been  increased 
from  the  Base  Run  by  the  conservation  release  in  Combination  One.  The  high 
flows  in  the  duration  curve  have  been  reduced  marginally  In  Combination  One, 
accounting  for  the  refilling  of  the  reservoir.  This  change  in  the  flow 
characteristics  then  affects  everything  further  downstream  on  the  Lehigh 
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River  as  well  as  the  mainstem  of  the  Delaware  River.  Of  particular  concern 
is  Trenton.  While  the  higher  flows  of  the  duration  curve  are  calculated  to 
be  the  same  as  the  basin  run,  the  lower  flows  have  been  augmented  by  100 
cfs.  This  augmentation  is  due  to  the  additional  reservoir  releases  needed 
to  maintain  a  flow  objective  at  Trenton.  A  Trenton  flow  objective  was  not 
used  in  the  Base  Run.  Similar  characteristics  to  the  duration  curves  can 
be  found  in  the  low  flow  frequency  tables  of  Trenton  and  Beltzville  for  the 
two  runs. 


Combination  17  with  all  the  reservoirs  on  line  affects  the  down¬ 
stream  characteristics  differently  for  each  of  the  reservoirs.  Using  the 
duration  curves  of  Combination  17  with  the  Base  Run  results.  Appendix  E  and 
Appendix  C  respectively,  comparisions  of  the  reservoirs  affects  can  be  made. 
Beltzville  and  Nockamixon  react  the  sameway:  increasing  low  flows  and  de¬ 
creasing  high  flows.  The  other  three  reservoirs,  Prompton,  FE  Walter  and 
Hackettstown,  react  just  the  opposite:  high  flows  increase  while  low  flows 
remain  approximately  the  same.  The  larger  high  flows  in  the  duration  tables 
are  due  to  the  releases  being  made  for  the  Trenton  objective  by  a  single  re¬ 
servoir  each  day,  raising  the  number  of  high  flows  considerably. 

Because  Prompton  Reservoir  is  located  in  the  Upper  Basin,  Montague's 
flow  will  also  be  altered.  First  of  all,  the  duration  curve  below  Prompton 
is  changed  by  the  reservoir's  operations.  The  low  flows  are  approximately 
equal  in  each  case.  The  high  flows  have  been  increased  by  the  releases 
necessary  to  maintain  Trenton's  flow.  Moving  down  the  basin  to  Montague, 
the  low  flow  frequencies  have  reflected  the  new  characteristics  due  to  the 
reservoir  operation.  All  of  the  flows  increase  from  the  Base  Run  to  Combina¬ 
tion  17.  However,  because  of  the  usage  of  flow  classes  rather  than  indivi¬ 
dual  flow  values  for  the  calculation  of  the  duration  curve,  the  number  of 
flows  in  each  class  remained  nearly  the  same.  Therefore  the  resulting 
duration  curves  are  identical. 

Trenton's  flows  are  also  of  particular  interest.  The  duration  curve 
of  Combination  17  shows  the  higher  flows  have  decreased  and  the  lower  flows 
have  increased  from  the  Base  Run.  Because  of  the  maintenance  of  the  Trenton 
flow  objective,  the  lower  flows  seen  in  the  Base  Run  are  supplemented  by 
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releases  from  the  reservoirs.  The  higher  flows  in  the  Base  Run  are  reduced 
as  more  and  more  reservoirs  are  put  on  line.  The  natural  high  flows  are 
retained  in  storage  by  the  reservoirs  to  replace  the  water  used  to  augment 
Trenton’s  flow.  The  low  flow  frequencies  of  each  run  have  similar  charac¬ 
teristics  as  the  duration  curves:  the  lower  flows  are  increased  and  the 
higher  flows  are  decreased  in  Combination  17  compared  to  the  results  of  the 
Base  Run. 


